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ADAPTATION IN FEEDBACK CONTROL SYSTEMS 


BY 
IAN McCAUSLAND'! 


ABSTRACT 


The scope of the word ‘‘adaptive”’ in the literature of feedback control systems is 
discussed. Approaching the subject from biological adaptation, the author shows 
that any control system which makes effective use of feedback is adaptive, and there- 
fore recommends that ‘adaptive’ should not be used to denote a separate class of 
feedback control systems. 

INTRODUCTION 


The word ‘‘adaptive” has recently become prom ‘ent in the litera- 
ture of feedback control systems. It is generally supposed that, 
through the study of adaptive control systems, ideas from adaptive 
processes in biological systems can be applied to the solution of prob- 
lems in engineering control systems. Unfortunately, however, papers 
on the subject of adaptive control systems have not given a clear and 
unequivocal definition of what adaptation in a control system really is. 
The purpose of this paper is to discuss the meaning of the term, with a 
view to clarifying its proper usage in the literature. 


ADAPTATION IN BIOLOGY 


Since it is generally assumed that the term ‘‘adaptive’’ is used in 
a sense similar to that in which it is used in biology (1, 2),? it is logical 
to begin by studying the literature of biology in order to discover the 
essential characteristics of adaptive behavior. One definition has been 
given by Ashby (3), who states that ‘‘a form of behaviour is adaptive if 
it maintains the essential variables within physiological limits."’ It 
seems obvious from this definition that any control system which makes 
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effective use of feedback can be described as being adaptive. The 
definition adopted, however, should not depend on one opinion alone, 
and so a typical adaptive process in biology will be considered in order 
to determine what are the essential characteristics of adaptation in such 
a situation. 

Suppose, for example, the process of temperature control of warm- 
blooded animals is considered. One particular aspect of this control 
process serves as an illustration of what is meant by adaptation, and 
this will now be discussed. A fur-bearing animal can control the rate 
of heat dissipation from its body by erecting or relaxing the hairs of its 
fur; this process is describable in terms of a simple feedback control 
system. Furthermore, it is known that animals which live in arctic 
climates grow heavier and thicker coats of fur than those of similar 
species which live in more temperate climates ; this is certainly an adap- 
tation, where “adaptation” is defined according to Webster (4) as 
“modification of an animal or plant (or of its parts or organs) fitting it 
more perfectly for existence under the conditions of its environment.”’ 

It is the author’s contention that these two processes, the simple 
adjustment of existing fur and the adaptive growth of a more suitable 
coat of fur, are essentially the same from the control engineer’s point 
of view. The control engineer is not concerned with whether or not 
any chanves occur in the animal’s weight or appearance, but only with 
whether or noi i hese changes affect its body temperature ; he will there- 
fore view each of .:ese processes as a change in the thermal insulation 
between the animal’s body and its environment, the fact that one 
process is more permanent than the other being merely a difference in 
time scale. 

The study of biological adaptation therefore indicates that any feed- 
back control system may be described as being adaptive. It is, how- 
ever, desirable to examine the field of control systems in a general way 
in order to determine whether the term “‘adaptive”’ is relevant or not, 
without any preconceived ideas as to what is meant by adaptation in 
any other discipline. 


ADAPTATION IN CONTROL SYSTEMS 


It is obvious that adaptation must involve change, but it is equally 
obvious that control of any kind must involve change. The basic 
question therefore is: what kind of change is adaptive, and what kind 
is not? Authors of previous papers on the subject have answered this 
question by saying that a parameter change is adaptive, whereas a 
change in signal is not. There appear to be two basic objections to this 
somewhat arbitrary separation of types of change. In the first place, 
if a quantity in a control system is automatically and continuously 
adjustable, it is highly doubtful that it can in any real sense retain the 
name of parameter. In the second place, the variation of parameter 


| 
| 
a 
: 


Sept., 1959.] 


ADAPTATION IN FEEDBACK CONTROL 145 


can only be effected by the application of a suitable signal at the ap- 
propriate point in the system. 

An “‘adaptable” transducer is generally represented in a diagram as 
a block with two inputs and one output (2). One of the inputs is the 
conventional input signal, whereas the other input is a signal, derived 
from some suitable source, which varies the transfer characteristics of 
the transducer. If an external observer studied the behavior of that 
part of the system by observing the conventional input variable and the 
output he would find the behavior indeterminate, and if the inde- 
terminacy were beneficial he might say that the behavior was adaptive. 
On the other hand, if the observer studied its operation as a two-input 
one-output device he would find it quite determinate. Furthermore, 
the observer would notice that if either input was kept fixed at a par- 
ticular value there would be a definite transfer relationship between the 
other input and the output, and if he had no further knowledge he 
could not say which was the parameter-varying input. The system 
can be thought of as a one-input device with an adaptive feature or as 
a two-input device with no adaptation. If the two-input case is chosen, 
the source from which the ‘‘parameter-varying”’ signal comes is merely 
a compensator. That is to say, whether the control is called adaptive 
or not depends on the system chosen, where system is defined according 
to Ashby (3) as ‘“‘any arbitrarily chosen set of variables.”’ 

In short, the only noteworthy feature about the control systems 
which have been called adaptive seems to be that they automatically 
vary some arbitrary quantities which happen to remain fixed in the 
conventional feedback control system. It seems very doubtful that 
this is sufficient ground for glorifying the new class of control system 
with an adjective which implies a higher level of control. In any case, 
the control problem is rarely amenable to these methods. The usual 
situation is that one has a transducer or process which has an unalter- 
able transfer function, the designer’s problem being to synthesize a 
compensating device which will allow the desired over-all performance 
to be achieved. Therefore any apparent variations of parameters willi 
be in the controller or compensator ; a controller which behaves in such 
a way is described by Johnson (5) as an adaptive controller. However, 
the controller is merely a device which operates on certain input vari- 
ables in such a way as to produce a desired set of output variables, 
although its operation may be described by saying that its characteristics 
are alterable in accordance with certain conditions. The alterable 
characteristics must be adjusted, directly or indirectly, by the inputs 
to the controller, so that the output is a function of the present and past 
values of the input variables. The controller may operate as a sequen- 
tial transducer, but it is doubtful whether this feature makes the con- 
troller adaptive. 

In June, 1957, the author contributed to the discussion of a paper on 
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‘‘Adaptive Servomechanisms” (1), criticizing the definition of ‘‘adap- 
tive’ for reasons similar to those outlined in this paper; that is, that any 
feedback control system may be said to be adaptive. In April, 1958, 
the author was present at a panel discussion on the subject of ‘Control 
Systems which Automatically Optimize Performance,’’ which de- 
generated into a discussion of the meaning of ‘‘adaptive,’’ and which 
showed that the increase in the number of papers on the subject in the 
intervening months had tended to increase the confusion in the terminol- 
ogy. Inthe report (6) on the conference at which this panel discussion 
took place, no mention is made of the subjects discussed by the panel, 
though the report does quote a delegate from the U.S.S.R. as saying in 
an interview that ‘‘In the Soviet Union, as here, there is the argument 
as to what is ‘self-adaptive, self-adjusting.’ ’’ It seems surprising that 
so little notice appears to be taken of the importance of formulating a 
satisfactory definition of adaptive behavior in control systems. 

Aseltine et al. (2) have discussed the usage of ‘‘adaptive’’ by listing 
the various types of adaptation which have been studied. Their first 
category is Passive Adaptation, and the opening sentences of their 
discussion of this category are as follows: 

“The classical example of passive adaptation is feedback itself. 
The reduction of effects of element changes through the use of feedback 
is well known. There are, however, special configurations which have 
inherent adaptive behavior beyond that found in conventional feedback 


systems. 

“The degree of adaptive behavior in a linear feedback system can be 
increased by prefiltering the input signal appropriately.”’ 

In other words, the authors are saying, the so-called adaptive sys- 
tems differ from the conventional feedback control system in degree 
but not in kind; this means that if the class of adaptive control systems 
is treated as a sub-class of all feedback control systems the dividing 
line can only be chosen arbitrarily. 

Drenick and Shahbender (1), in their rebuttal to the author's dis- 
cussion of their definition of ‘“‘adaptive,’’ have given the following list 
of terms which have been used to describe systems similar to those 
described by them: “input-controlled,”’ “‘sympathetic,’’ ‘“‘learning,”’ 
“self-adjusting,”’ “‘self-designing,’’ and ‘‘adaptive.”’ It is not the pur- 
pose of this paper to discuss all the possible terms which might be ap- 
plied to such systems, except to suggest that the use of six different 
terms to mean the same thing is not conducive to the clarity of thought 
and precision of expression which are essential to the dissemination and 
furthering of scientific knowledge. 


CONCLUSION 


Because of the foregoing arguments, it is the conclusion of the 
author that ‘‘adaptive’’ should not be used as a technical term to specify 


3 


Sept., 1959.] ADAPTATION IN FEEDBACK CONTROL 147 


a special class of feedback control systems. However, the word may 
be quite appropriate as a descriptive term applicable to all feedback 
control systems, possibly a better word than ‘‘feedback”’ itself. The 
only thing the word ‘‘feedback”’ tells one about a system is that it has a 
certain kind of topographical representation; on the other hand, 
“adaptive” tells one something about its function, which is what 
primarily interests the control engineer. In other words, the sig- 
nificant thing about the feedback system is not that it happens to 
possess a feedback loop, but that it possesses the ability to adapt to 
disturbances and to variations in the parameters of the forward loop. 
All control systems which make effective use of feedback are adaptive ; 
therefore the word ‘‘adaptive’”’ should not be used to describe an arbi- 
trarily-selected class of such systems. 

It is hoped that this discussion of the subject will help to clarify 
the proper usage of the word “‘adaptive,’’ and that some attempt will be 
made to establish a consistent and logical terminology which will 
minimize ambiguity and confusion. 
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MODES AND NATURAL FREQUENCIES OF SUSPENSION 
BRIDGE OSCILLATIONS 


BY 
D. B. STEINMAN! 


SYNOPSIS 

The natural frequencies of potential modes of suspension bridge oscillations con- 
stitute basic data required in studies for aerodynamic stability. 

The formulas presented herein were derived by the author in 1941-1943 and have 
since been tested, for simplicity and practical usefulness, by fifteen or more years of 
applicacion in the author's office practice. They have also been independently 
checked, for validity and accuracy, by others using more difficult formulas involving 
transcendental (complex or hyperbolic) functions. The author’s formulas appear 
to be the simplest and most practical thus far developed for the calculation of the 
natural frequencies of all possible modes of suspension bridge oscillations. 

For convenience of reference and application, these formulas are here assembled, 
derived and presented in definitive form, together with numerical examples to illustrate 
their practical application. The numerical examples are based on computations for 


the Mackinac Bridge.” 


DEFINITIONS AND BASIC RELATIONS 


If N is the natural frequency (cycles per second) of any harmonic 
oscillation, and w is the circular frequency (radians per second), then 


o=2¢N = (1) 
m 


where m = w/g is the oscillating mass (in the case of a span, m is the 
mass per unit span-length) and K is the ‘‘spring constant,”’ or ratio of 
the elastic force of restitution (Aw), per unit span-length, to the dis- 
placement » at any point. Hence, by definition, 


(2) 


K is also termed the ‘Coefficient of Rigidity’’ of the oscillating system, 
a significant criterion for aerodynamic stability. 

For an oscillation of any span (/) in the form of a simple sine-curve 
(see Figs. 1 and 2), with m segments (half-waves) in the span, the am- 
plitude » at any point x is given by: 


(3) 


. 
n = asin 7 x 


1 Consulting Engineer, New York, N. Y. 
? Epitor’s Note: See this JouRNAL, Vol. 262, p. 453 (1956). 
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where a is the amplitude at mid-segment. The displacement », at any 


instant ¢ is 
9 sin wt (4) 


for n, = 0, att = 0. 


SIDE SPAN MAIN SPAN SIDE SPAN, 


2 n=4 
Fic. 1. Antisymmetrical modes. 


In the derivations that follow, useful integrals are: 


ndx = 0 (for even) 


ndx = (for odd) 
0 


rn 


f = 
0 


Useful differentials are : 


é 

| | | | 

} 
| | 

( a ) | 
| 
n=O 
| 

| 
(6) 
Fy =wn (at = 0) (8) 
7 
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d’n 
dx* 


d‘n 
dx* 


SIDE SPAN__ MAIN SPAN SIDE SPAN 
SIDE SPAN __MAIN SPAN _ 


Fic. 2. Symmetrical modes. 


The cable curve (with sag f) is assumed parabolic : 


y = 4x(] — x) 


J. 
(10) 

4 
(11) 

| 
(a) 

| 

| n,=| n=1,3 n=! | 
“ing 
‘ 
(c) 

met | n=3,! 
| 
(d) 
| 
| n=l n=3,5 | 

| 
| 

kee 
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yielding 
13) 
dx? 
By Eq. 9, the force equation at any instant on any unit length of 
span is 
— w'nm + Kn, = 0 (14) 
yielding 
K = mw. (15) 
This is the derivation of Eq. 1. 
BASIC FORMULA FOR K 


For any element (dx) of the span-length, the potential energy at full 
displacement (n, = ») will be 


W, = 3nAwdx. (16) 
By Eqs. (8) and (15), the kinetic energy at zero displacement (n, = 0) 
will be 


W, = = }Kn'dx. (17) 


Equating the energy transformation (W, = W,) and applying the 
summation to cover all spans and all component sine-curves, the basic 


formula for A is obtained : 
l 
> f nAwdx 
0 


If a span has a coupled synchronous longitudinal oscillation of am- 
plitude e, then W,, given by integration of Eq. 17, is augmented by 


K 


AW, = = = 1 (19) 


and Eq. 18 becomes 


l 
p f nAwdx 


K = 
f ndx + ¥ el 


GENERAL FORMULA FOR K 


The fundamental equilibrium equation for a suspension bridge is 


Aw = = - + yAH — 


I5! 
| 

tee 

‘ (20) 
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where // is the horizontal cable tension, A// is the increment produced 
by displacement, J is the moment of inertia of the section of the stiffen- 
ing truss, and E is the modulus of elasticity for the truss. 

By Eqs. 10, 11 and 13, Eq. 21 yields: 


Aw = = Hn +" EIn + AH. 


For any case in which AH = 0, Eq. 2 or 18 then yields: 


(23) 


Equation 23 is the simplest and most useful formula for K. It is 
exact for all even values of n (n = 2, 4, 6, ---), the so-called ‘‘antz- 
symmetric modes” (Fig. 1), since, for these modes, AH = 0. 

The value of K for n = 2 is the ‘‘Coefficient of Rigidity’’ of the sus- 
pension bridge, a useful criterion for aerodynamic stability. 

In Eq. 23, the H-term and the EJ-term represent the respective 
contributions of the cable and the stiffening truss to the rigidity of the 
suspension bridge. The ratio of the truss term to the total is desig- 
nated by R and is significant for aerodynamic stability. 

From Eqs. 20, 22, 6 and 7, 


EI + 32 (2 AH 


(wal) + (we’l) 


Wo 


This is the most general form of the equation for K. It covers all 
possible combinations of spans, modes, and superimposed harmonics. 

In Eq. 24 and the following equations, and in all similar formulas 
involving summations of functions of the amplitudes, it is not necessary 
to know the absolute magnitudes of the several concurrent amplitudes 
(a and e), but merely the relative magnitudes. 

Equation 24 is completely general for any number of spans of any 
respective lengths (/), simple or continuous, with any respective cable 
sags (f), including straight or unloaded backstays, for any combina- 
tions of values of m and a in simple or superimposed harmonic oscilla- 
tions in the different spans, with any concomitant longitudinal oscilla- 
tions (e) in any spans, with different values of w, H/, EJ, in different 
spans, and with any value of AH producing harmonic cable stretch. 

In the denominator of Eq. 24, wo is the value of w for the main span. 


(22) 
ay 
l 2 l 4 
GE 
K 
2 
: 
: 
; 
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Example 1 
Antisymmetric Modes—K by Eq. 23 
(Fig. 1) 
For the Mackinac Bridge: 
= 3800 ft. 1; = 1800 ft. 
w = 9400 lb. ft. 
= 48,478,000 lb. (for two cables) 
= 3,768 X 10° lb. ft.? (for two trusses) 
(K in Ib. /ft.?) 
By Eq. 23 By Eq. 1 
= 132.5+ 28.2 = 160.7 = 7.09/min. 
530.1 + 450.6 980.7 N = 17.50/min. 
147.7 + 35.0 182.7 = 7.55/min. 
= 590.7 + 559.4 = 1150.1 N = 18.95/min. 


EVALUATION OF AH 


The virtual work represented by the vertical displacements 7 is 


Sf! 
Wi= f ndx. 


The virtual work represented by the cable stretch due to A// is 


AH ['ds 
dx 
Hence, for W,; = W2, with Eq. 6, 


af 


nl 


_ “f 
AH = > 


tds3 
0 dx? 


sp 


In these equations, and those that follow, E. is the elastic modulus 
and A, the section of the cable, and y is the inclination of the cable- 
chord. Also, © denotes summations covering all similar expressions 


in all spans, and +’ denotes a summation covering only the odd values 
of n in the several spans. 


COMPOSITE MODES 


A simple or “primary mode’ is one represented by a single sine- 
curve—hence, single values of m and of a, as in Eq. 23. The term 


153 
(25) 
= ) 
where 
> 
(28) 
| 
: 
: 
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“composite mode’ is used herein to denote oscillations comprising a 
plurality (superimposed or coupled, or both) of simple modes, as in 
Eq. 24 (see Fig. 2). 

From Eq. 24, using the expressions of Eqs. 23 and 27, the general 
formula for K for any single component mode (any simple sine-curve) is 


na nl 


where 
(30) 


a dimensional constant for the structure, and A, is the value of K given 
by Eq. 23. For ©’ = 0, Eq. 29 reduces to Eq. 23. 

For symmetric modes (necessarily composite), as illustrated in 
Fig. 2, Eq. 29 is valid for each component mode; and, with the con- 
sistent (normal) relative amplitudes substituted, it must yield the same 
value of AK (and, hence, the same frequency N) for each component 
mode. To satisfy this condition, Eq. 29 yields two simple governing 
relations between AK and the respective values of m, a, and K,, of the 
component modes : 


and 


na « (32) 


Equation 31 yields the desired value of A for the normal (balanced) 
composite mode, and Eq. 32 gives the relative amplitudes of the com- 
ponent modes. 

In practical application, the largest number of component modes 
necessary to be considered and substituted in Eq. 31 is three, and that 
is the case of the lowest symmetric mode (m = 1,3; and”, = 1, Fig. 2d). 
In the higher symmetric modes, the superimposed harmonics prove 
negligible, and only two primary A,-values need to be substituted, one 
for each span (Example: » = 3 and m, = 1, Fig. 2d). Accordingly, 
Eq. 31 is in reality a quadratic, or cubic, equation for AK. Each real 
root of this equation will represent a normal composite mode with a 
mathematically and physically consistent set of relative amplitudes. 

In practical application, Eq. 31 lends itself to a simplified numerical 
solution by an expedited method of trial substitution. It is the sum 
of two or three fractions equated to unity, with the sole unknown, K, 
in each denominator. If the mode represented by the first fraction is 
dominant, substitute an approximate trial value of A in the other term 


is 
£ 
| 
: fx SENS. 
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or terms, and solve the first fraction for K. The first or second trial 

will generally suffice. Similarly, the one or two other roots are obtained 

by treating the other fraction or fractions, respectively, as dominant. 
The subscripts in @,, 7, fi, etc., designate values for the side spans. 


Example 2 
Symmetric Composite Modes—A by Eq. 29 
Mode: = 1, 3; = 1 
(Figs. 2b, 2c) 

Components: 

Amplitudes : 

K, (by Eq. 23): 34.9 440.8 
(Compare Example 1.) 


For the Mackinac Bridge: 
744 sq. in. L, = 9140 ft. 
C = 2866 |b/ft? 


= 264.0 


By Eq. 31: 


264.0 264.09 2(125.0) 
K — 34.9 K — 440.8 © K — 182.7 


Three solutions: 

Eq. 32 
a;/a = — 0.063 
a,/a = — 0.776 


Dominant 


n 


=> > 0.285 
a,/a3; > 0.474 


Example 3 
Mode: m = 3, 5; m, = 1 (Fig. 2e) 
Components: 
Amplitudes: 
K, (by Eq. 23): 440.8 1928.6 


Eq. 31: 264.025, -2(125.0) _ 


ny = 1 
ay 
182.7 
f 5 
C= = 125.0 
l 
n=3 K = 405.6 
fs (N = 11.26) 
nm, = 1 K = 696.7 aa; = 0.777 
es (N = 14.75) a3/a, = 0.670 on 
ay 
182.7 
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Three solutions: 
K 


Dominant 


= 3 363.2 a, = — 1.290 
(N = 10.65) as ‘a3 = 0.030 

n= 5 1941.2 a;/a, = 0.014 
(N = 24.62) a;/as = 0.036 

ny = 1 537.6 a3 ‘ay = 1.222 
(N = 12.96) as/a, = — 0.051 


For n = 5, dominant, the small value of a;/a; indicates that the super- 
posed harmonic m = 3 is negligible. Similarly, for m = 3, dominant, 


the superposed harmonic = 5 is negligible. 


Example 4 


Mode: m = 3, m, = 1 (Fig. 2d) 


Amplitudes: as ay 

K, (by Eq. 23): 440.8 182.7 

Eq. 31: 264.0/9 2(125.0) _ 
K — 440.8 © K — 182.7 


Two solutions: 


Dominant K Eq. 32 

n= 3 363.8 a;/a3 = — 1.276 
(N = 10.66) 

ny = 1 539.1 a;/a, = 1.209 


(N = 12.98) 


Note that these results are practically identical with those obtained 
in Example 3. The inclusion of the third component (m = 5) in Ex- 
ample 3 was unnecessary ; the effect is almost imperceptible. 


Example 5 


Mode: nm = 5, m, = 1 (Fig. 2f) 


Amplitudes: ds ay 

K, (by Eq. 23): 1928.6 182.7 

264/25 2(125.0) 

K — 1928.6 K — 1827! 


Two solutions: 


aq. 32 
0.035 


K 
1940.9 
(N = 24.62) 

430.9 a;/a; = — 0.033 
(N = 11.60) 


Dominant 


a 
Eq. 32 
4 
m=d 
1 = 1 
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These results for nm = 5 are practically identical with those obtained 
in Example 3. 

For n = 3,5, ---, with m, = 1, the inclusion of a third component 
in applying Eqs. 29-32 is unnecessary. 


STRAIGHT OR UNLOADED BACKSTAYS 


For the case of straight backstays, assuming zero sag (f: = 0), all 
of the foregoing equations, including Eqs. 29-32, are applicable without 
any modification. The sole effect is in the reduced numerical value of 
L, given by Eq. 28, reflected in the increased numerical value of C, 
given by Eq. 30. 

For application to unloaded backstays with the small cable-sag (/1) 
considered, Eqs. 29, 31 and 32 require a simple modification. Instead 
of the usual relation of equal unit dead loads in main and side-spans 
(w, = w), the weight of the cable becomes the sole dead load (w,) in 
the side spans. Consequently, in reducing Eq. 24 to Eq. 29 for a side- 
span mode, the multiplying factor w/w, must be retained. The C 
term is unchanged, since C, (Eq. 30), contains the factor f//? which is 
proportional to w;/wo. Accordingly the sole change in Eqs. 29, 31 and 


32 is to replace K,, for any side-span component by a K,,. This yields: 


UW) n,Q, nl 


nN 
na 


Example 6 
Straight Backstays 


For the hypothetical case of the Mackinac Bridge with straight back- 
stays, assume L, = 7540, C = 3474, C f/l = 320.0. 
For the lowest symmetric mode (m = 1, 3), Eq. 31 becomes: 


320.0/9 
K — 440.8 


320.0. 
K — 34.9 


1 


| 
— C= 
vi 4 
Wo 
ra — 
w 
= ny 
Wy 
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Two solutions: 


Eq. 32 
a;/a = — 0.581 


+ 0.582 


Dominant K 
1 292.9 


3 538.4 a/a3 


n 


ll 


n 


Compare Example 2. The straight backstays raise the correspond- 
ing values of K from 99.5 to 292.9, and from 405.6 to 538.4. 


Example 7 


Unloaded Backstays (with Cable Sag Considered) 


Wi) 2670 
= = (0.284. Ass = 1000 ft. 
Wo 9400 4. Assume /; ft. Hence fi//; 


= 0.0069, C f/l = 320.0 and C f,/1; = 24.0. 


For this case, 


For the mode: » = 1, 3; 2, = 1, (Compare Example 2.) 


320 3209 2(24.0) 


Three solutions: 


Dominant K Eq. 35 
n= 1 273.3 a;/a = - 0.474 
} = 9,24) = — 0.644 


= 0.364 
— 1.147 
0.152 
0.117 


Compare Examples 2 and 6. Unloaded backstays with cable sag 
considered yield A-values somewhat lower than the assumption of 
“straight backstays,’’ but substantially higher than the A-values with 


loaded side spans. 


COMPACT FORM OF GENERAL EQUATION FOR K 


The general formula for A, Eq. 24, may be rewritten in a simple 
compact expression of equal generality : 


> (Kal) 
») (36) 
(wal) + — ¥ (wel) 
Wo We 


in which A for each span or for each component mode has been com- 
puted separately by Eq. 29 or any other pertinent formula. 

Equation 36 may also be written directly from the energy relation, 
Eq. 20, with the aid_of Eqs._2 and _7. 


: 
= 
mn = 3 496.8 a a3 
(N = 12.46) a; 
ny 1 752.4 a ay co 
(N = 15.33) a3 ay sp 
4 
ties 
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With w equal in all spans and with e = 0, Eq. 36 reduces to the 
simpler form: 


=F (37) 


which may also be written directly from the energy relation, Eq. 18. 

Equation 37 is an exceedingly convenient and useful formula for 
practical application. It may be used for beams and other structural 
combinations, as well as for suspension bridges. 


OSCILLATIONS WITH LONGITUDINAL MOTION OF SPAN 


If mn = 2, 6, 10, ---, without participation of the side spans, the 
longitudinal motion A(//2) of the midpoint of the cable may be im- 
parted to the suspended span through any physical coupling at midspan, 
as when center ties are provided. For this case, by Eqs. 6, 25 and 27, 


16 f 


38 
a nl (38) 


and Eq. 36 yields: 


in which Ko is given by Eq. 23. 

If the two side spans oscillate in opposite phase, with the main span 
free of vertical oscillation (a = 0), Fig. la, the parallel and synchronous 
movement A/, of the two tower tops produces longitudinal motion of 
the main span equal to A/;. For this case, 


(40) 


and Eq. 36 yields: 
] 


K (41) 


d, 


in which A, is given by Eq. 23 for n, = 1, 3, 

Equations 39 and 41 express an apparent reduction in K as the 
equivalent of an increase in the effective value of the oscillating mass m 
in Eq. 1. 

In any segment (/) of a parabolic cable, by Eqs. 25 and 6, 


Sf 
I? Jo nl 


Ke 
(39) 
142 (+ ry’ = 
r nl 
3 
16 fy 
: 
a, il 1 
i 
: 
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e=A (5) and Eq. 42 yields Eq. 38. 


For /,, m1, fs, e = Al,, and Eq. 42 yields Eq. 40. 


CABLE ANCHORED AT MID-SPAN 


If the cable is anchored at midspan, as by center ties, and longi- 
tudinal oscillation of the span is prevented (as by friction brakes, or 
automatically as in the case of torsional oscillations), the antisymmetric 
modes n = 2, 6, 10, ---, will induce side-span participation, thereby 
raising the value of K. For this case, Eq. 24 or Eq. 29 is applied over 
the half-span and side span, or from the point of fixation to either 
anchorage, using the obviously modified values of m, f, and 1. (Use 
1/2 for 1, n/2 for n, f/4 for f, }L, for L,, 2C for C, and 2C f,/l; for C fi/li; 
a, K,, and C f/l are unchanged.) 


Example 8 
Antisymmetric Mode with Cable Anchored at Mid-Span 
Mode: m = 2, m,; = 1 (Fig. 10) 
K, = 160.7, 182.7 (Example 1.) 
Data: C f/l = 264.0; 2C fi/l, = 249.9 
Applied over the half-span and one side span, Eq. (31) for this case 
becomes: 
264.0 249.9 


K — 100.7 + K— 182.77! 


Two solutions: 
Dominant K N 
n=2 171.7 
nm, = 1 685.5 14.63 


Compare Example 1. The anchorage of the cable at midspan 
augments K and N. For n = 2, by inducing side-span participation, 
K is increased from 160.7 to 171.7 (increasing N from 7.09 to 7.32 per 
minute). For m,; = 1, by inducing main-span participation, K is 
increased from 182.7 to 685.5, augmenting N from 7.55 to 14.63 per 
minute. 

APPROXIMATE METHODS 

The last term in the numerator of Eq. 24 represents the effect of AH. 
It is also represented by the C-term in Eq. 29. AH is ordinarily either 
zero (in antisymmetric modes) or of relatively minor magnitude (except 
form = lorn = 1,3). For AH = Oor neglected, K is given by Eq. 23 
for a simple mode or by Eq. 36 or 37 for a composite mode. 


160 (J. FL 
Ln f 
For 2’ 4’ 
3 
a;/a: = — 1.000 
a2/a, = 0.958 
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There will be no cable stretch between anchorages (AH = 0) if 


_ 
2. = 0 (43) 


in each span—also, of course, for antisymmetric modes. 

A simple short-cut for estimating the value of K is to assume that 
AH = 0. Consistent relative amplitudes are then given by Eq. 43 or 
44, and their substitution in Eq. 36 or 37 yields K with sufficient 
accuracy for all practical purposes. 

Forn = 3,5,7, ---, the governing requirement AH = 0, represented 
by Eq. 44, may be satisfied, without involving the adjoining spans, by 
superimposing a single-segment oscillation of relative amplitude —a/n, 
yielding, by Eq. 3, 


over all spans, or if 


(45) 


Substituting K, and the relative amplitudes in Eq. 37, 


\n? + 1 

Equation 46 may be generalized for any case of superimposed undu- 
lations of an even number of different odd values of n. If the same 
relation is maintained—namely, that the respective amplitudes a of the 
component harmonics are proportional to their values of n—the general- 


ized form, Eq. 37, reduces to 
(47) 


For the three spans coupled (to oscillate in unison) as by tower-top 
movement, with AH = 0, Eq. 37 may be written in the form: 


Ko+ 

K= 48 

i+g¢ 

in which Ky and K, are the K,-values by Eq. 23 for the respective 
spans, and 


(49) 


te 
: 
= asin — sin—-. 
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The main and side spans may oscillate in unison, with a single 
segment or any odd number of segments (m and m,, respectively) in each 
span (Fig. 2). If the bridge is symmetrical, the relative amplitudes to 
satisfy Eq. 43, substituted in Eq. 49, yield 


n,*l 
2n*l, 


(50) 


K is then given by Eq. 48. 

If 1,/n, = 1/n, the main and side spans are in resonance and 
K Ko K,. 

If m = 2, 6, 10, ---, then } A(//n) is zero for the full span but is 
not zero for the half-span. If the midpoint of the cable is not free to 
take the longitudinal motion A(//2), as when effectively anchored by 
center ties, the lengthening and shortening of the respective half-spans 
will tend to be taken up by side-span participation to produce Al, 
= A(//2) at each tower top. For this case, Eq. 43 is applied over //2 
and /,. The resulting values of the relative amplitudes, substituted in 
Eq. 49, yield 
= 


q (51) 


K is then given by Eq. 48. 


Example 9 
Symmetric Mode—Approximate Method by Eqs. 45-47 

For the symmetric composite mode n = 3, 1, Eq. 46 yields: 

K = 271.7 + 128.5 = 400.2 (N = 11.18). For this case (but with 
side-span participation, Fig. 2c), the more exact method (Example 2) 
yields K = 405.6 (N = 11.26). The error of the approximate method 
is less than 1.5 per cent in the value of K and only 0.7 per cent in the 
value of N. 


Example 10 
Antisymmetric Mode—Approximate Method by Eqs. 48-51 


For the antisymmetric composite mode n = 2, m, = 1 (Fig. 16), with 
center ties forcing side-span participation, Eq. 51 yields g = 1.056. 
With A, = 106.7 and K, = 182.7 (from Example 1), Eq. 48 then yields 


K = 172.0 (N = 7.33) 
These results are almost identical with K = 171.7, N = 7.32, given 
by the exact method Eq. 31, in Example 8. 
SHORT-CUTS FOR ESTIMATING K AND N 


Using the cable term alone, Eqs. 23 and 1 yield the convenient ap- 
proximate relations: 


: 
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K = 1.234n? 


(53) 


in which f is simply the sag of the cable. (Note: In Eq. 53, P is in 
seconds and f is in feet.) These simplified formulas, Eqs. 52 and 53, 
may be used for quickly calculating or estimating the coefficient of 
rigidity K and the natural period P or frequency N of an unstiffened 
suspension bridge or of one with negligible stiffening system, for all 
values of m > 1. 

Similarly, using the girder term alone, Eqs. 23 and 1 yield the con- 
venient approximate relations: 


(54) 


(55) 


in which », is the central deflection theoretically producible (with E 
constant) by the dead load w acting on the span as a simple beam. 
(Note: In Eq. 55, P is in seconds and »,, is in feet.) Equations 54 and 55 
may be used for simple girder spans and for self-anchored suspension 
spans, for all values of n. 

For a stiffened suspension bridge (except the self-anchored type), 
the composite rigidity AK is the sum of the respective contributions 
(Eqs. 52 and 54) of cable and stiffening girder. Hence, by Eqs. 53 
and 55, 

¥ Nw 

Equations 53 and 55 are convenient formulas to remember for 
quick estimates. 


(56) 


Example 11 
K and N by Short-Cut Method 


For the Mackinac Bridge: 
w = 9400 Ib. /ft. f = 350 ft. Nw = 6773 ft. 
For n = 2, Eq. 52 yields K 
Eq. 54 yields K = 


Total K 


163 

w 
2 (52) 

: P=—=-Vf 
N n 
: 

K = 1.27n‘ — 
Nw 
and 

2 
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This result is almost identical with the exact value given by Eq. 23. 
(See Example 1.) 
Also, for n = 2, Eq. 53 yields N = 6.41 
Eq. 55 yields N = 2.92 
Hence, Eq. 56 yields N = 7.05/min. 


This result is also almost identical with the exact value given by Eqs. 
23 and 1. (See Example 1.) 


Similar excellent checks of Eqs. 52—56 are obtained for all even values 
of n. 


TORSIONAL OSCILLATIONS 
Torsional oscillations may be treated simply as opposed-phase 
vertical oscillation of the two cables. The angular amplitude at any 
point of the span is given by 


(57) 


in which d is the width of the bridge measured center to center of girders 
or trusses. 

Ordinarily, the torsional stiffness of the suspended span is relatively 
negligible. Hence, unless substantial torsional stiffness is provided by 
special construction (as with a double system of lateral bracing), the 
values of K and N determined for vertical oscillations, multiplied by 
simple dimensional constants, yield with sufficient accuracy the corre- 
sponding values for torsional oscillations. The spring constant of the 
span in torsion, or the coefficient of torsional rigidity, will be 


( 5) K (58) 


and the natural frequency will be 


(59) 


(60) 


in which r is the polar radius of gyration of the mass of the cross section 
of the span. It is important to note that, for torsional resistance of a 
suspension bridge, the polar moment of inertia (r*m) pertains, not to 
a localized built-up cross section at the panel points, but to the average 
cross section throughout, including the long ranges between panel 
points. For the usual proportions, b/(2r) = 1.3 to 1.5. 
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The polar radius of gyration r is given by 
& (mr) (61) 


m 
in which the component r for the cable masses is 6/2, but for all elements 
of the suspended section it is the polar radius to each element from the 
center of oscillation of the section. 


TORSIONAL RESISTANCE OF WIND TRUSSES 


In the derivation of the foregoing formulas for K, the effect of web 
shear in the stiffening truss or girder has been neglected. The effect 
of the web shear contribution is found to produce a comparatively 
infinitesimal or negligible reduction in K. The normal omission of the 
web shear contribution from the K formulas is completely justified. 

If, however, wind trusses are provided in the horizontal planes of the 
top and bottom chords (or flanges) of the stiffening trusses (or girders), 
new relations between the bending and shear components come into 
play in torsional displacements. ‘The interaction of the horizontal and 
vertical trusses (or girders), with common chords, imposes a relative 
reversal of sign in the bending components of the horizontal trusses, 
with a consequent enforced increase in all of the shear components in 
the deflection and energy equations. In such case the shear terms are 
no longer negligible but contribute a valuable increase in K. 

Accordingly, a highly effective method of augmenting resistance to 
torsional oscillations is by providing two planes of lateral bracing, at or 
near the top and bottom chords (or flanges) of the stiffening trusses (or 
girders), respectively, so as to secure the integral effect of a hollow rec- 
tangular section in torsion. The cross section is assumed rigidly con- 
structed or braced so as to preserve its rectangular form. The resultant 
transverse shearing forces in the four planes of the tubular member, 
when it is subjected to torsion, will then contribute very material 
increases in K, and R. 

Let » and z represent the vertical and lateral displacements, respec- 
tively, at any section of the chord (or flange). Each is the algebraic 
sum of bending and shear contributions: 


2= Z = Sm + (62) 


Between the corresponding elements in the vertical and horizontal 
planes (distinguished by subscripts where necessary), the five governing 
static conditions of equal rotation, equal chord stress, identical chords, 
shear equilibrium, and ratio of shear deflections, in the order named, 
yieid the following five simple basic relations: 

b 


Me 
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in which d is the vertical truss depth and 


(64) 


(non-dimensional). A, is the vertical web section of the girder, multi- 
plied by E,/E; or, in the case of a truss, A, is the equivalent solid web 


section and is given by 


A, = Aga (cos’ ¢ sin ¢) (65) 


in which A, is the sectional area of the diagonal member or members in 
a truss panel, ¢ is the inclination of the diagonal from the vertical, and 
E, is the shear modulus of elasticity. (For steel, E,/E may be taken 
as 0.385.) A, is similarly defined in each horizontal plane of bracing. 
In Eqs. 63, note the enforced reversal of sign for the z,, ratio; this is 
significant. 

From Eqs. 62 and 63, the component deflections 9», 7., 2m, Zand z 
are easily written as simple ratios of », containing only }, d, and k. 
Also, substituting , in the shear geometry of a girder web yields 


2k A dn (66) 


In the case of a truss, the geometry of a truss panel, with the aid of 
Eq. 65, yields the identical formula, Eq. 66. 

Substitution of the foregoing relations, Eqs. 62-66, in the energy 
equation, as in Eqs. 16-22, yields the desired K-formula: 


in which the three terms represent the respective contributions of cable 
displacement, chord (or flange) stresses, and diagonal (or shear) stresses. 
This last term is new, and constitutes a very significant addition to K 
and an even more significant addition to R. 

In Eq. 67, the moment contributions of the vertical and the hori- 
zontal trusses are equal; but the corresponding coefficient 2 is omitted 
from the second term of Eq. 67, because the equality in this case 
represents a duplication of identical chords counted in both the hori- 
zontal and vertical trusses. 

Without the two planes of horizontal trussing, A, = 0, k = 0, 
and Eq. 67 reduces to Eq. 23. 
pe | Equation !67 replaces K, in previous formulas. If AH is not zero, 
add_the AH/ term as in{Eq. 29. For simple antisymmetric modes ( = 2, 
4, 6, ---), Eq. 67 is complete. 
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K should preferably be written to represent the complete cross 
section of the bridge, rather than the value for one side or one cable. 
Accordingly, in Eq. 67, H refers to the two cables, J and A, refer to the 
two vertical trusses or girders, and A, refers to the two horizontal trusses. 

For the usual case of two cables and stiffening trusses, the fota/ 
value of K multiplied by 62/4 gives K, (in torsion), as in Eq. 58. 

Equation 67 assumes the two vertical stiffening trusses (and cables) 
oscillating in opposite phase. Consequently Eq. 67 does not apply to 
ordinary vertical oscillations, and does not represent vertical rigidity of 
the structure. 

By Eq. 67, the addition of wind trusses (in two horizontal planes) 
is a highly economical and effective method of providing desired re- 
sistance to torsional oscillations. It is vastly surpassed, however, in 
economy and efficiency by the simple device of installing transverse 
diagonal stays, between the cables and the trusses, as described and 
successfully applied by the author in 1943. 


Example 12 
Antisymmetric Mode ( = 2) in Torsional Oscillations 
For n = 2, Example 1, by Eq. 23, gives 
K = 160.7 (N = 7.09) 
Hence, by Eqs. 57-61, with 6 = 68 ft., r = 30.4 ft. 
K, = 185,800 (N, = 7.93/min.) 
With center ties, causing side-span participation (Fig. 1b), Example 8, 
by Eq. 31 applied over the half-span and one side span, gives 
K = 171.7 (N = 7.32). 
Hence, by Eqs. 57-61, 
K, = 198,500 (N, = 8.19). 

For the case with two horizontal planes of lateral bracing, we use the 
data of Example 1 together with d = 38, b = 68, A, = 15.62, A, = 23.60, 
k=0.472. With these data, Eq. 67 gives: 

K = 132.5 + 3.6 + 1588.4 = 1724.5. 
Hence, by Eqs. 57-61, 
K, = 1,993,500 (N, = 25.96). 

Compared with K, = 185,800, this represents a 973 per cent increase 
in K, (and a 227 per cent increase in N,) due to the provision of the double 


system of lateral bracing. 
With center ties, causing side-span participation (Fig. 1b), in addition 


to the double system of lateral bracing, Eq. 31 applied over the half-span 


he 
aa 


168 


D. B. STEINMAN 


and one side span, gives: 


264.0 4 2(125.0)  _ 
K — 1724.5 ' K — 1922.0 — 


This yields, for n = 2 dominant, 


K = 1807.5 (and K = 2353.0 for mn; = 1 dominant). 
Hence, by Eqs. 57-61, 
K, = 2,089,500 (N, = 26.58). 


Compared with the corresponding value of K, = 198,500, this repre- 
sents a 953 per cent increase in K, (and a 225 per cent increase in N;,) due 
to the provision of a double system of lateral bracing. 


Example 13 
Symmetric Mode (n = 1) in Torsional Oscillations 
For the composite mode m = 1, m; = 1, Fig. 2a, Eq. 31 yields K = 100.3 
(N = 5.60). Hence, by Eqs. 57-61, 
K, = 115,900 (NM; = 6.26). 
For the composite mode n = 1, 3, with m; = 1, Fig. 2b, Example 2, by 
Eq. 31, gives K = 99.5 (N = 5.57). Hence, by Eqs. 57-61, 
K, = 115,000 (N; = 6.23). 
This value is slightly lower than K, = 115,900, and represents the 


desired minimum. 
For the case with two horizontal planes of lateral bracing, with the 


data of Example 12, Eq. 67 gives: 
K, = 33.1 + 0.2 + 397.1 = 430.4 
for n = 1, also K, = 3890.4 for n = 3, and K, = 1922.0 for n,; = 
With these values of K,, for the composite mode n = 1, 3, with m 
Fig. 26, Eq. 31 becomes: 
264.0 264.0/9 2(125.0) _ 
K— 4304 * K — 38004 


K — 1922.0 © 

This yields, for m = 1 dominant, K = 649.4. (The relative ampli- 
tudes are given by Eq. 32.) Hence, by Eqs. 57-61, 
Compared with K, = 115,000, this represents a 553 per cent increase in 
K, (and a 156 per cent increase in N,) due to the provision of a double 
system of lateral bracing. 

For the composite mode m = 1, m, = 1, Fig. 2a, Eq. 31 takes the 
simpler form: 
264.0 2(125.0) _ 
K — 430.4 © K — 1922.0 — 
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This yields, for » = 1 dominant, K = 651.0. Hence, by Eqs. 57-61, 


The effect of omitting the component m = 3 thus appears to be prac- 
tically negligible. The lower value (K; = 750,700) represents the de- 
sired minimum. 
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EFFECT OF TORSIONAL STIFFNESS OF THE TOWERS 


Suspension bridge towers are normally designed as flexible canti- 
levers, fixed at the base. Their contribution to K, thefstiffness of the 
bridge against vertical oscillations, is small and usually negligible. 

The torsional stiffness of the towers, however, offers a larger contri- 
bution to K,, the resistance of the bridge to torsional oscillations. This 
contribution may be material, even without special design of the towers 
for greater torsional stiffness. In comparison with other methods, such 
special design of the towers would usually be an uneconomical means of 
increasing the torsional rigidity of the bridge. 

A couple formed by a horizontal cable force P at each saddle will 
twist the tower about its vertical axis to produce a longitudinal dis- 
placement e at each saddle. The torsional stiffness of the tower is 
denoted by 
P 


Kr (68) 


Let 
= L, for main span. 
L,, = L, from anchorage to tower. 


(Lio + Lar 41.) 
AH, = AH from anchorage to tower. 


AH, = AH, + P = AH in $ main span. 
Considering the main span, 
_ 8 AH Lio 
~A E.A. (69) 


Considering the side span, 


e=-— B (70) 


in which 


av 
= f idx = (cos 1). 


|| 
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Equating the two expressions for e, the foregoing equations yield : 


al, = — + (BAe + KrLw)B) 


AH, + P C(EA)B- 


in which 
E.A. 


U=L,+ 


The total potential energy due to cable stretch (A7/) in the whole 
bridge and including the torsion in the fwo towers will then be: 


(74) 


The first three terms in the bracket represent W, due to cable 
stretch, and the fourth term represents W, due to torsion in the two 
towers. 

For Kr = 0, Eq. 74 reduces to the AH term in Eq. 24. 

The potential energy due to truss distortion has already been 
determined, as represented in Eq. 67, and gives the first term of the 
numerator of the following equation. The kinetic energy of the towers 
is negligible. From Eq. 74 and W, = W, = 3K ¥ a’l, we obtain the 
equation for K: 


in which K,, is given by Eq. 23 or 67 for each component mode. 

Equation 75 gives the numerical value of A in terms of the relative 
amplitudes (a,/a). For a composite mode consisting of one main- 
span mode (amplitude a) and one side-span mode (amplitude a,), the 


. .. OK . 
partial derivative _* 0 yields K and a,/a for the mode n dominant ; 
a 


‘ 
and the partial derivative = 0 yields K and a/a, for the mode n, 
a; 


dominant. These values will be virtually unchanged if an additional 
mode of higher ” is superimposed. For a composite mode consisting of 
three or more simple modes superimposed, the resulting simultaneous 
equations may be solved graphically. 


a+. 
2 
: 
; 
: 
(73) 
: 
= 
(75) 
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In applying Eq. 75, K, and A. should be taken as the values for the 
complete cross-section of the bridge (two cables and two vertical trusses 
or girders), in order that Eq. 58 may then be applied to the resulting 
value of K to give K,. 

From K given by Eq. 75, the values of K, and N, are then given by 
Eqs. 58 and 59. 

Equation 75 assumes the two vertical stiffening trusses and cables 
oscillating in opposite phase. Consequently Eq. 75 does not apply to 
ordinary vertical oscillations, and does not represent vertical rigidity of 
the structure. 


Example 14 
Antisymmetric Mode with Torsional Stiffness of Towers 
(Fig. 1b) 
For the Mackinac Bridge: 
9140 ft. E.A. = 20.82 (10)9 Ib. 
= 2029 ft. Kr = 873,000 Ib./ft. 
= 2541 ft. U = 10,005 ft. 
48 478,000 Ib. f/l = 0.0921 fis 0.0436 
> Eq. 67, A, = 1724.5 for nm = 2 and K, = 1922 for m 1 
(Compare Example 12.) 
For the antisymmetric mode, nm = 2, m; = 1, Eq. 75 then yields 
K= 7.764a? + 1.584a,a + 7.9360," 
0.0038a? + 0.0036a;? 


OK 
= vields: 
da 


For n = 2, dominant, 
K = 1895, and a,/a = — 0.711. 
Hence 
K, = 2,190,600 (N, = 27.21). 


For m; = 1, dominant, = 0 yields: 


OK 
0a; 


K = 2352.8, and a/a, = 0.674. 
Hence 
K, = 2,719,800 (N, = 30.32). 


Comparison with the corresponding values in Example 12 shows that 
the effect of including the torsional stiffness of the towers is an increase 
of 4.8 per cent in K, (and 2.4 per cent in N,) for the mode n = 2, 
dominant. 
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Example 15 
Symmetric Mode with Torsional Stiffness of Towers 


(Fig. 2a) 
Data as in Example 14. 
By Eq. 67, K, = 430.4 for n = 1, 
and K, = 1922.0 for n, = 1. 


(Compare Example 13.) 
For the symmetric mode m = 1, m, = 1, Eq. 75 then yields 


K 2:846a* + 1.584010 + 7.9360;" 
~ + 0.0036a:? 


For n = 1, dominant, on = 0 yields: 


K = 718.2 and a;/a = — 0.148 


Hence 


K, = 830,200 (N, = 16.75). 


For n; = 1, dominant, oe 0 yields: 
0a, 


K = 2235.2 and a/a, 


0.140 


Hence 


K, = 2,583,900 (N; = 29.55). 


Comparison with the corresponding values in Example 13 shows that 
the effect of including torsional stiffness of the towers is an increase of 
10.3 per cent in K, (and 5.0 per cent in N,) for the mode m = 1, dominant. 

If we superimpose the harmonic » = 3 on the oscillation m = 1, , = 1, 
the numerical results obtained above are very slightly modified. The 
value of A, for m = 1, dominant, is reduced from 830,200 to 826,800; 
and the value of K, for ny = 1, dominant, is reduced from 2,583,900 to 


2,579,000. The reduction in K, and N, is insignificant, and such refine- 


ment is not justified. 


EFFECT OF CENTER OF SUSPENSION 

The effective center of suspension of a suspension bridge cross 
section in torsional oscillation is the midpoint of the line joining the two 
points of connection of the suspenders to the outside girders or trusses. 
The torsional stability of the section is increased by raising this center 
of suspension relative to the center of gravity of the suspended section. 
This is confirmed by simple model tests. The spring constant K, is 
augmented by a torque of gravity restitution ; a gravity pendulum effect 
is added to the torsional spring constant. 

If h is the vertical height of the center of suspension above the center 
of gravity of the suspended section, the vertical rise of the center of 


> 
: 
a 
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gravity of the suspended section, as it swings in an arc like a pendulum, 
is, by Eq. 57, 
@ 2h 
= —_ = - 2 7 
v=h (76) 


and the potential energy is augmented by 


AW, = w,-v (77) 
in which w, is the suspended dead load (exclusive of the cables) per unit 
length. Accordingly, by Eqs. 16-20, and 57-61, the value of K, given 
by Eq. 58 is augmented to 


+ wa. (78) 


The total energy W required for any angular amplitude is increased in 
the same proportion, and consequently the aerodynamic instability, 
measured by the rate of amplification, is reduced in the same ratio. 

In computing the polar radius of gyration r by Eq. 61, the center of 
oscillation of the suspended section was specified. If the polar moment 
of inertia mr* has been computed about the center of gravity of the 
suspended section, a correction must be made when the center of oscilla- 
tion is at a height 4 above the center of gravity. The polar moment of 
inertia is augmented by 


A(mr*) = m,h’. (79) 


By combining Eqs. 78 and 79 we may write the equivalent value 
of K;: 


b 2 

(5) 

K, = mrtw? = 
*) 


(80) 


In Eq. 80, the numerator represents the augmented potential energy 
and the denominator represents the augmented kinetic energy. The 
kinetic energy is augmented by the horizontal velocity of the center of 
gravity along the arc of the pendulum swing. 

Instead of increasing A by raising the points of suspender connection, 
as by deeper girders or by brackets, the same improvement may be 
accomplished by lowering the center of gravity of the suspended section. 
Accordingly, torsional rigidity and stability may be improved by chang- 
ing from deck or half-through construction to through construction. 
The author has demonstrated this conclusion by model tests. An un- 
stable deck truss section became aerodynamically stable when simply 
inverted. 


; 
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Example 16 

Effect of Center of Suspension 
For the Mackinac Bridge: 

b = 68 ft. 

r = 30.4 ft. 

w, = 6530, h = 6.96 ft. 


Hence Eq. 80 becomes: 


w = 9400, 


(b/2)*K + 45,449 


For n = 2 (See Example 12), A, is thus increased from 185,800 to 


223,100, an increase of 20 per cent. 
For n = 1 (See Example 13), K, is thus increased from 115,900 to 
155,700, an increase of 34 per cent. 
For the much higher values of K,, with two horizontal planes of lateral 


bracing, the effect of the raised center of suspension (h = 6.96) is found 


to be negligible. 
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VARIATIONAL PRINCIPLE FOR TRANSPORT PHENOMENA 


BY 
G. E. TAUBER!' 


SUMMARY 

The role of the variational principle in transport theory and its application to the 
calculation of the transport coefficients are described. A generalized variational 
principle including the effect of an external magnetic field and lattice vibrations is 
discussed. The variational principle for inter-band scattering is derived, and is 
generalized by the introduction of an external magnetic field and lattice vibrations, 
leading to the most general variational principle including inter-band scattering, in 
addition to an external magnetic field and lattice vibrations, or any combination of 
these. The transport coefficients are expressed in terms of energy integrals, which 
for the various cases take the form of ratios of infinite determinants, the number and 
form of their matrix elements depending on the interactions. A formal solution is 
given by the Ritz method for the different cases of interest. 


I. INTRODUCTION 


Many physical problems can be formulated in terms of a variational 
principle, the solution of which is often simpler and more accurate than 
the original formulation. This is due to the fact that the variational 
function is not as sensitive to fluctuations and errors as the differential 
(or integral) equation obtained upon variation. For practical calcu- 
lations use is made of the Ritz variational method in which the unknown 
function is expanded as an infinite series with coefficients determined by 
the variational method (1). 

In the case of transport phenomena the fundamental equation de- 
termining the distribution function, and thus the various transport 
phenomena, is an integro-differential equation known as the Boltzmann 
equation (2). We can write this equation as an operator equation 
of the form 

L(F) =X (1) 


where L represents the integral operator due to the change of the distri- 
bution function by collisions, and X the force due to the external field 
and temperature gradients. It has been shown by Kohler (3) that this 
equation is equivalent to a variational principle, for which the varia- 
tional function is simply the inner product 


(F: F) = (2) 


' Department of Physics, Western Reserve University, Cleveland, Ohio; and The Franklin 
Institute Laboratories, Philadelphia, Pa. 
? The boldface numbers in parentheses refer to the references appended to this paper. 
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subject to the subsidiary condition 


(F,X) = ff xd = (F; F) on account of (1). (3) 


It has been pointed out by Ziman (4) that the variational function 
(F; F) is not simply a mathematical device, but of physical interest as 
it is proportional to the rate of entropy production by collision processes, 
and thus affords a connection between microscopic statistical theory 
and irreversible thermodynamics. The positive definiteness of that 
function based on microscopic reversibility of the operator L assures 
us also that the Onsager reciprocity relations (5) among the various 
transport quantities are always satisfied (6). 
If we now expand the function F in a series of the form 


F= > ae (4) 


where « is the energy, variation with respect to the coefficients a, will 
result in a system of linear equations for these quantities 


a, = > d,a, (5) 


where the coefficients a, and d,, can be evaluated from a knowledge of 
the form of L and X. The various transport coefficients, such as elec- 
trical and heat conductivity, can then be obtained in terms of energy 
integrals K which turn out to be just the subsidiary conditions (3). It 
has been shown by Enskog (7) that these energy integrals can be ex- 
pressed in terms of infinite determinants, the elements of which can be 
evaluated for a given energy-momentum dependence. 

If deviations from equilibrium of the phonon distribution function 
are also taken into account this results in a system of coupled integral 
equations, which in operator form are given by 


X = Li(F) + L.(N) (6) 
Y = L3(F) + 


Here N represents the phonon distribution function, in addition to the 
previously discussed carrier distribution function F. Furthermore, in 
addition to the operator L, = L (of 1) there also appear the operators 
L., Ls and Ly, describing the interaction of the electrons with the 
phonons and phonons with phonons. A variational approach to this 
problem has been suggested by Ziman (4), and recently Dorn (8) ob- 
tained a variational principle and applied it to include the effect of 
longitudinal polarized vibrations for spherical energy surfaces. The 
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variational function is given by the sum of terms of the form (2) 
(F; F)i + (F;N)2+ (N; F)s + (7) 


where the subscripts 7 = 1 to 4 refer to the appropriate operator L,. 
The symmetry relations amongst the integral operators assure us of 
the positive definiteness of the variational function (7) and the existence 
of the Kelvin relations amongst the transport coefficients. If the lattice 
distribution function N is also expanded in terms of an infinite series 
analogous to the corresponding one for F (4) an application of the Ritz 
variational method will then again result in a system of linear equations 
for the unknown coefficients. The energy integrals are again given in 
terms of infinite determinants, the elements of which now depend on a 
knowledge of the operators L; (¢ = 1 to 4) and the functions X and Y. 

In the presence of an external magnetic field the above theory has 
to be modified. The Boltzmann equation for the carrier distribution 
function acquires an extra term, which represents the effect of the 
“magnetic scattering.” Taking the symmetry property of the mag- 
netic scattering operator under reversal of the field into account, it is 
still possible to obtain the Kelvin relations for that case (9). Although 
the magnetic field does not contribute to the entropy production it is 
possible to formulate a variational principle (10) which leads to the 
correct transport equation. Its solution is again achieved by an applica- 
tion of the Ritz method, except that two separate expansions for the 
conjugate sets F+ and F- must be used, depending on the direction of 
the magnetic field H and —H. 

We have shown (11) how a generalized variational principle can be 
found, which is valid for an arbitrary direction of the electric field and 
polarization of the lattice vibrations in the presence of an external 
magnetic field, and does not depend on any special form of the energy 
surface. The various transport coefficients, both for thermoelectric and 
thermomagnetic phenomena, are then obtained by the Ritz method in 
terms of infinite determinants without requiring an explicit solution of 
the transport equation. 

In some cases of physical interest, notably in semiconductors, one 
has more than one energy surface, and it becomes necessary to consider 
the effect of transitions between bands. The appropriate operator 
equation is then given by 

Lan(F*) = Xm (8) 


n 


where the subscripts m, m refer to the various bands. The variational 
function for this case turns out to be a sum of terms of the type given 


earlier (2) 
= (Fs) (9) 


3 
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subject to the subsidiary conditions 


(F, X) = ¥ (F*, X,.) = (F; F) on account of (8). (10) 


This variational principle can again be solved by an application of the 
Ritz method where the coefficients a, (of 4) now acquire an extra sub- 
script indicating the band m and become a,,. Carrying out the 
variation with respect to these parameters one again obtains a set of 
linear simultaneous equations for these unknowns of the form (5), 
the difference being the appearance of extra indices in the quantities 
a, and d,, indicating the bands. The current density J and heat 
current density W may be considered as a sum over bands over the 
appropriate expressions for the individual bands 


J=ZJj. and W => W.. (11) 


The transport coefficients are again expressed in terms of energy in- 
tegrals K which can be obtained from infinite determinants without the 
necessity of explicitly solving the system of equations for a,,,, the only 
difference being that now there appear more matrix elements including 
the contributions from all different bands m, n. 

A further generalization of the theory can be achieved by considering 
the effect of an external magnetic field or/and the presence of lattice 
vibrations for this case. If the non-equilibrium phonon distribution 
function is taken into account a system of coupled equations of the 
form (6) is obtained 


Xn = +E 
(12) 


and a variational principle can again be formulated, the variational 
function being of the type (7) where each term is now a sum over bands 
of the type (9). The effect of an external magnetic field again intro- 
duces an extra scattering operator L’ into the formalism and produces 
a separation of the solution into two conjugate sets depending on the 
direction of the magnetic field H. The energy integrals, and hence the 
transport coefficients, can again be obtained by the Ritz method and 
expressed in terms of infinite determinants, containing now many more 
matrix elements describing the effect of the interaction of the electrons 
with the lattice, phonon-phonon interaction, and the various bands. 
The purpose of this paper is to give as complete a discussion as is 
possible of the variational principle in transport theory and its applica- 
tion to the calculations of the various transport coefficients. Starting 
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with the derivation of the variational principle developed by Kohler® 
we include the effect of an external magnetic field or lattice vibrations, 
and finally combine both these effects in a discussion of the generalized 
variational principle. The variational principle for inter-band scatter- 
ing is set up next, and is afterwards generalized by the introduction of 
an external magnetic field and lattice vibrations, so that we finally 
obtain a variational principle which includes inter-band scattering, in 
addition to an external magnetic field and lattice vibrations. In order 
to present a systematic presentation, Section II is devoted to the 
Boltzmann equation and contains a discussion of the various cases, v7z. 
effect of magnetic field, lattice vibrations and inter-band scattering, and 
a combination of any or all of these. In Section III, the appropriate 
variational principles are derived and their solution indicated by an 
application of the Ritz method. Finally, in Section IV, starting from a 
general expression of the transport coefficients in terms of energy in- 
tegrals, it is shown how these can be expressed in terms of infinite 
determinants, which differ from each other in the number and nature 
of their matrix elements. In order to evaluate these matrix elements, 
the form of the integral operators must be known, together with the 
dependence of energy on momentum. Work along these lines is in 
progress, and has been reported elsewhere (12), while here we have only 
been concerned with presenting a systematic discussion of the varia- 
tional principle applied to transport phenomena and in particular, 
expressing the various energy integrals necessary for the calculation of 
the transport coefficients, in terms of infinite determinants for the vari- 
ous cases of interest. 


Il. BOLTZMAN TRANSPORT EQUATION 
1. Nature of Transport Equation 


For a good analysis of transport phenomena it is necessary to de- 
termine the distribution of the particles both in regard to position and 
velocities. The fundamental equation determining the distribution 
function f(v, r) is an integro-differential equation known as the Boltz- 
mann equation, which describes the probability of a particle being in 
a given state, and is therefore based on the one-electron model. It is, 
however, quite simple to generalize the distribution function, and hence 
its change with time to a n-dimensional phase space in which it is now 
a function of the co-ordinates and momenta of the m particles. Re- 
cently, several authors have adopted the approach of the density matrix 
in quantum statistics to derive the transport equation from general 
statistical consideration. However, for the present presentation, we 
will assume the validity of the standard theory and postpone a more 
rigorous derivation of the transport equation to a future discussion. 


3 Loc. cit. 
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If the number of particles which lie in an element of dr of configura- 
tion space and in dk of k-space at a time ¢ is given by 


(4x3)-1f (kk, r, t)dkdr 


(12) 


then at time ¢ + dt according to Liouville’s theorem they will lie in an 
equal volume of six dimensional phase space about k + kdt, r + vdt. 
The difference, therefore, between f(k + kdt, r + vdt, t + dt) and 
f(k, r, 2) must be due to the particles scattered into the new positions 
in the time dt. The total change of the distribution function f due to 
collisions with the lattice (and each other) and external fields must be 


zero, so that 
(13) 


(Of corr. + (Of /At) tictas = 0 


where (0f/0t) represents the change due to collisions and (0f/0t) tietas 
the change due to external fields, which is given by 


df/dt = af/at +k- Vif +v- Vf. (14) 


Here ak is the change of the impulse and equal to the external force 


ak = cE + e(v/c) X H (15) 


where E denotes the electric field, H the magnetic field, e the charge and 
c the velocity of light. The velocity of the electron v is expressed in 
terms of its wave vector k through 


(16) 


av = Vue(kz, R,, 


(In the steady state one also has df/dt = 0.) 


a. Collision Operator 


Since binary collisions of electrons are neglected, the number of 
particles per unit volume with velocities in the range dv which have 
their velocities changed to lie in dv’ in time dt must be proportional to 
f(v, r,t) and to the interaction W(v, v’, r), which is independent of f 
and only depends on the interaction mechanism. Similarly, the number 
of particles which are scattered into the range dv are proportional to 
f(v’, r,t) and the interaction W(v’,v, r). The net difference between 
these two quantities, integrated over dv’, determines (0f/0t) .o11., so-that 


(Of /At) ou. = f(v’, r, v, r) — f(v, r, )W(v, v’, (17) 


The equilibrium distribution function fy must be such that (0f0/0t) eon. =9, 
and thus 


W(v, v’, r)/W(v’, v, r) = fo(v’, r)/fo(v, r). 
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Now according to classical mechanics, f» is a function of the energy only, 
and in fact, is the Maxwell-Boltzmann function. Also, the collisions 
are supposed to be elastic, so that energy is conserved and f(v, r) 
= fo(v’,r). Hence W(v, v’, r) = W(v’, v, r) and 


(Of/dt) con. = f [f(v’, r,t) — f(v, r, v’, r)dv’. (18) 


In the case of electrons, Pauli’s exclusion principle has to be taken 
into account. We describe the states of the electrons by the wave vector 
k (rather than the velocity v) and suppose the W (k, k’) is the probability 
that an electron makes a transition from the state k to the state k’. 
As a transition can only occur if the corresponding state is unoccupied, 
the number of electrons, per unit volume, gained by dk from dk’ in 


time dt is 

W(k’, k)f(k’)[1 — f(k) \dkdk’dt 
and the number of electrons lost from dk to dk’ is 
W (k, k’) f(k)[1 — f(k’) }dkdk’dt 


so that the net change in f is given by 


(Af/At) f We’, k’)f(k)(1—f(k’) dk’. (19) 


This expression must vanish in the equilibrium state, where there exists 
detailed balancing and therefore 


V(k, k’) = W(k,k’)fo(1 — fo’) = W(k’, k) fo'(1 — fo) = V(k,k’) (20) 


in which we have simplified the notation by writing fo’ to denote fo(k’), 
and which leads to the Fermi function 


fo = [expB(e-—S) B= I/kT. (21) 


In problems of conductivity we are interested, primarily, in very small 
displacement from the equilibrium state. Thus we are justified in 


writing 
f=fothr = fo — Qdfr/de = fo + Qfo(l — fo)B (22) 


which follows from the form of the Fermi function fo (22). Hence, if 
terms involving Q? are neglected, which is equivalent to neglecting 
squares and products of the electric fields and temperature gradients, 
(19) becomes 


(Of/9t) con. = ef V (k, k’)[Q(k) — Q(k’) Jdk’. (23) 
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b. Operator Form of Transport Equation 

The transport equation is then obtained by equating (23) describing 
the effect of the collisions with the lattice to (14) due to the external 
field. From (21) one also obtains the relations 
Vito 


Vilo 


(24) 


= hvodf de 
(0fo/OT)V.T + (0fo/d5) VE 


af./aT 


(Afo/de), = — Afo/de. (25) 


With the help of these relations, one can write (14) in the steady state 
up to first orders in f 


df/dt = — vf eE, + Td/ax,(t/T) + ¢/T(@T/ax.)] 


i=1 


+ (Afo/de) (vx (26) 


The form of this expression suggests the separation of the equation 
into two parts, one being proportional to « and one independent of it. 
Thus, if we write for 


Q = ¥ {[eE, + T0/dx.(¢/T) + (A/T) foi}, (27) 


equations (26) (with /7 = 0) and (23) can be replaced by 


Li fis) = = (27a) 


L( foi) = vOfo/de = (27b) 


where L is the integral operator defined by (23) 


L(F) = 6 f Vs, — Fe) (28) 


The symmetry property (20) leads to the symmetry relation 


(F;G) = (G; F) 


where we have defined the inner product 


(F;G) = af F(k)L(G)dk. 


: 
bara 
: 
30) 
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The truth of this result can be seen by writing 
(F;G) =8 f f V (ic, k’) F(x) [G (kk) — G(e’) 
= 8 f f Ge) 
on account of (20). Hence, by addition 
= 8 f f Vie, Fe) IEG — 
which is symmetrical in F and G and hence establishes the relation (29). 


c. Effect of the Magnetic Field 


In the presence of an external magnetic field Eqs. (27) acquire an 
extra term arising from the last term in (26) and which can be taken 
into account by adding to the operator L the operator L’ 


L'(Q) = (e/hc) (Afo/de)(v X H)ViQ. (31) 


The operator L’ is not symmetric, but satisfies the relation 


FL'(G)dk = — f GL’ (F)dk 


and, on reversal of the magnetic field, //, 
FL,"(G)dk = GL"(P)dk (33) 


where L,’ and L_’ denote the operator L’ with +H and —H respec- 
tively in the expression (31). These results follow from the fact that 
the integral 


f v-[(v X 


vanishes. It is also to be noted that both v and H are solenoidal. 

From (33) it is seen that the solutions of the transport equation 
depend on the direction of the magnetic field. On reversing the mag- 
netic field the electrons will curve in the opposite direction and the 
condition of microscopic reversibility (20) then becomes 


V+(k, k’) = V-(k’,k). (34) 


It is thus possible to divide the set of solutions f,; and f.; into conjugate 
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sets f;,+ and f;;- (j = 1,2) for which we have analogous to (27) the 
following system of equations 


= L(fii*) + (fist) 
= L(fai*) + 


(35a) 
(350) 


together with the corresponding ones for a magnetic field in the opposite 
direction, which are obtained from (35) by changing all + into — 
quantities. The system of equations (27) or (35) can be obtained by 
means of a variational principle, which will be discussed in Section ITI. 


2. The Lattice Distribution Function 


The transition probability, W(k’, k), depends mainly on the inter- 
action of the electrons with the lattice. If m,; denotes the distribution 
function of the lattice vibration, where q describes the state of the 
lattice vibrations, 7 the polarization (j = 1, 2 for transverse and 
j =3 for longitudinal) the transition probability for a transition 
k — k’ = k + q with absorption of a lattice quantum is of the form 


(36) 


W (k, k’) = V(k, k’)n,6(ex — ex» + Av) 


and for a transition k — k’ = k — q with emission of a lattice quantum 


(37) 


W (k, k’) = V(k, k’)(n, + 1)6(ex — ex: — 


where the Dirac 6-functions ensure that energy is conserved in the transi- 
tion. In other words, (36) vanishes unless ex — ex + hy, = 0, and 
(37) vanishes unless ¢, — ex» — hy, = 0. Here », is the frequency of 
the lattice wave in the state q. 

If the effect of lattice vibrations is neglected, m, may be replaced 
by its equilibrium value 


(38) 


y = hv/kT. 


no(q) = (ev — 1)", 


However, in the general case considered here, both f and n,; depart 
from their equilibrium value, and in addition to (12), one has a trans- 
port equation for n, 


(an,/at)eon. = (AT (39) 
i=l 


where U’; are the components of the lattice wave velocity. Hence, the 
two functions f and m, must be determined simultaneously as the solu- 
tions of the two Boltzmann equations (12) and (39). 


: 
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a. Collision Operators 


For the two collision operators (Af/0t).ou. and (An/dt).ou., respec- 
tively, one obtains (taking into account the exclusion principle and the 
relations (36) and (37) and the reverse transitions k’ — k) : 


(Of /At) corr. = fitvae, k) fx (1 (mg +1) — V(k, fr (1 


X thy) + LV (k’, k) fe (1 
— Vik, k’) (me +1) (40) 


(On cot. = k) fe: (1 (mg +1) — V(k, kk’) 
X +hy,)dk. (41) 


In the equilibrium state, these expressions must vanish. Hence, putting 
f =fo and n, = mo, and using the energy equation ¢«’ = «+ hy we 
must have 


V(k’, k) + hv)[1 — fo(e) }(mo + 1) 
= V(k,k’) fo() [1 — fole + h»)] (42) 


V(k’, kt) fo(e — hv)[1 — fo(e) 
= V(k, k’)fo(e)[1 — fole — hv) }(mo + 1). (43) 


Both these restrictions are satisfied if V(k, k’) satisfies the symmetry 
relation (20) which expresses the principle of microscopic reversibility 
for the present problem and forms the basis of the Kelvin relations to be 
discussed later. 

If we now expand n,; in the same manner as has been done for f 


Ny; = No — N;Bdn/dy = no + Bno(no + 1)N; (44) 


the exact expressions (33) and (34) may be simplified. Keeping only 
terms linear in Q and N, we obtain 


(Of eon. -af V (k, k’)[m0(q) fo(k) (1 —fo(k’)) — +hv4) 
[Q(k’) (de) — N (a) fo(e’) (1 —Fo(e)) 
X 5(€x — —hv,)[O(k’) —Q(k) +N (q) (45) 
(46) 


where the left-hand side of these expressions is given by (84) and (88), 
respectively. They are now a system of integro-differential equations 
for the two quantities Q and N,j. 
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b. Operator Form 
The system of equations, (45) and (46) can be written in operator 
form 


= viOfo/Oe = + Leo(gis) (47a) 
= 0 = Lal(fis) + La(gis) (47d) 
Xo, = v0fo/de = Li( fei) + (47c) 
yoo = — = Lil foi) + La(ges) (47d) 


where the function N has been separated into two parts similar to Q 
(cf. 27) 
N = {[eE, + Td/dx'(¢/T) + (48) 


The integral operators L; (L = 1 to 4) are given by 


L\(F) = ef V (k, k’)[S(k, k’) + S(k’, k) J[F(k) — F(k’)]dk’ (49a) 


= B f V k’)[.S(k, k’) F(q) — S(k’, F(q) (495) 
L,(F) = B J V(k, k’).S(ke, — F(k’) Jak (49c) 
L,(F) = 8 f V (kc, k’).S(k, k’) F(q)dk (49d) 


(It is important to note that q = k’ — k is to be kept constant in (49c) 
and (49d) which are therefore functions of q.) S(k, k’) is given by 


S(k, k’) = mo(q) fo(k)[1 — fo(k’) — ex» + (50) 


corresponding to an absorption of a lattice quantum k > k’ = k + q, 
(36), while S(k’, k) corresponds to the emission of a lattice quantum 
k —k’ = k — q and has the transition probability given by (37). 

For compactness of notation, let us now introduce the generalized 
inner product analogue to (30) 


where the subscript ¢ indicates the operator L; defined by (49). The 
symmetry property (20) leads to the symmetry relations 


(F;G)i = (G; (52a) 
(F;G). = (G; F)s (520) 
(F;G), = (G; (52c) 


: 


Sept., 1959.] VARIATIONAL PRINCIPLE IN TRANSPORT THEORY 187 


The property (52a) follows from (29), while (52c) follows at once from 
the form of L,4(F) (49d). To establish (525) we write 


(F:G):= f F(k)L3(G)dk 


=p f f V (k, k’) F (ie) (ik, k’)G(q) — k’)G (q) 
since V(k, k’) = V(k’,k). Hence 
= 6 = G:F): 
which completes the proof. 


c. Effect of a Magnetic Field 


As in the case without lattice vibrations taken into account, the 
presence of an external magnetic field introduces an extra term in the 
Boltzmann equations (47) which again can be represented by an 
operator L’ defined by (31) and possessing the properties (32) and (33). 
The modified equations are then given by 


= Li(fist) + Lo(gii) + Le'(fii*) (53a) 

= L3(fist) + Lalgii) (53d) 

Li(foit) + Lo(gei) + Ly" (fast) (53c) 

= L3(foi*) + L4(g2:) (53d) 

with corresponding ones for a magnetic field in the opposite direction, 

which are again obtained from (53) by changing all + to — quantities. 

As L’ does not operate on N, it is not necessary to separate g;; and g2; 

into two conjugate sets, as must be done for f;; and f2;. The solution 

of this more complicated system of equations can also be obtained by a 
variational principle to be discussed in Section IIT. 

3. Interband Scattering 


In some cases of physical interest, notably in semiconductors, one 
has more than one energy surface, and it becomes necessary to consider 
the effect of transitions between bands. Before generalizing (27) to a 


} 
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multiple band model, we note that the operator equations (27) can be 
written as 


X ad Lux J (54) 


showing the connection between the forces X, and currents J,;. From 
the properties of L(Q) we can also infer that the coefficients Lyx form a 
positive definite symmetric matrix, with the symmetry relations 
Lux: = Linx. Indeed, this last result forms the basis of the Kelvin 
relations amongst the thermoelectric coefficients to be discussed later. 
The summation in (54) is to be carried out over all the states k. In 
order to take into account the other bands we simply extend the sum- 
mation in (54) over the appropriate states. Thus the system of equa- 
tions (27) has to be replaced by 


foi*) = ViemOfo/Oem = X2i™ (55d) 


where the operator L,,, is defined by 


Lnn(F) = ef V(kn, — F(kK,) dk, (56) 


and m and m denote the possible energy bands. Thus the effect of a 
scattering from band m to n is taken into account by the operator La, 
acting on the various states /" = F(k,) and producing a coupling 
between them. 

The symmetry property (20) in our present notation becomes 


V(k.,, = V(k,, Kn) (57) 


and again leads to the symmetry relation 


(F;G) = (G; F), (29’) 


where, however, now the inner product (F;G) is a sum of terms of the 
type (F";G") such that 


(F;G) = ¥ (F*;G*) (58) 


where we defined the product (F"; G") by 


(Fe. G*) = f 


= 
n 
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The proof of the generalized symmetry relation is obtained by writing 
again 


(F;G)=BD ff V(kn, k,) F(k,.)[G(k,.) — G(k,) (59D) 


a f k,,) F(k,.)[G(k,) — G(k,) Jdk,dk, (59c) 


on account of (57). Hence by addition 


ff V(Kn, K,)[F (Kn) — F(k,)] 
[G(k,.) — G(k,) Jdk,dk, (59d) 


which is symmetrical in F and G, as for each term (F"™;G") in the sum 
we have a corresponding one (F";G™), and thus establishes the truth 
of (29’). 


a. Effect of the Magnetic Field 

In the presence of an external magnetic field, the equations (55) 
acquire an extra term, or rather a series of terms, of the type (31). 
Thus the new modified operator equations are now given by 


> + Le (60a) 


+ Ln’ (fei**) (60) 


where the + and — signs again denote the direction of the magnetic 
field, and where the sets of functions f;,;" and f2," are now divided into 
conjugate sets f1;+ and f2;"* as before, and the operator L,,’ defined by 


= €/NC Of o/O€m(Vm X H) 


It should be noted that the operator L,,’ describes the effect of the mag- 
netic field on the band in question, and therefore only has a free index m, 
characteristic of the particular band. 


b. Presence of Lattice Distribution Function 


We have seen that the presence of a non-equilibrium lattice distri- 
bution function modifies the transport equation in two ways. First, 
the lattice distribution function satisfies a Boltzmann equation of the 
type (39) and secondly, the original Boltzmann equation for the carried 
distribution function acquires additional terms. The transition from 
band m to n is characterized by the emission (or absorption) of a lattice 
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quantum according to k, = k,, + q, and hence the interband scattering 
can be taken into account formally by replacing k and k’ in the expres- 
sions for the integral operators L; (49) by k,, and k,, respectively. We 
note also that in the case of L; and L, summation must be carried out 
over m and n, as the lattice function depends on q, and is therefore 
independent of m and n. 

The system of equations (47) is then replaced analogous to (55) by 
the following system of operator equations 


0 = + (£11) (616) 


where the integral operators L,,,' (¢ = 1 to 4) are now given by 


Lan'(F*) = B f V (ken, Kn) + Km) 

x [F(e) — F(k,) dk, (62a) 
Lunt(G) = ef V (Km, k,)[S(kem, ,)G(q) — Kn)G(q) dk, (628) 
La.(F*) = f V (Km, (im, — (62c) 
= af V (Km, k,)G(q)dk (62d) 


with 
= k,,) = no(q)fo(Km){1 — fo(Kn) ]6(€m €, + hv,) (63) 

for an absorption of a lattice quantum k,, > k, = k,, + q, and S(k,, k,,.) 
corresponding to the emission of a lattice quantum according tok,, > k,, 
=k, —q. We note that the operators L,,,' and L,,.4 have to be 
summed over both indices m and n, as the lattice distribution functions 
gi, and go, are independent of the interband interaction. This can also 
be seen by inspection of their operator forms, (59c) and (59d), which 
are integrated over both dk,, and dk,,. 

The appropriate generalized inner products analogous to (51) are 
now defined as 


(F;G); 


(F;G)/ = ¥ (F*; G) as” s=2,4 (645) 


~ 
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with 


2, 4. (65) 


(F";G) mai” = 


The need for differentiating between the two cases a and 6 arises from 
the fact that in the second case the functions G represent the lattice 
distribution function and are independent of the band, as can be seen 
by inspection of the corresponding operators L* and L‘* (626) and (62c). 
The superscript 7 gain indicates the operator L‘, while the subscript mn 
shows the bands. 

The symmetry property (57) again leads to symmetry relations of 
the type (52). 


II 


(52a’) 
(G; F)s (52’) 
(52c’) 


(F;G), 
(F;G): 
(F;G)s = (G; F)s 


the first of which follows from (29’) while the last can be reduced from 
the form of LZ‘. To establish (526’) we write again 


(F;G)2= F(Kn) Lm? (G)dKn 


X [S(Km, k.)G(q) — S(K,, k,.)G(q) Jdk,.dk,, 


ff 


X F(kn)G(q) — S(Kn, k,) F(k,.)G(q) Jdk,.dk,, 
since V(K,,, k,) = V(k,,k,.). Hence 


(F:G):= f f V (Km, (kn) — F(k,) 


which completes the proof. In the case of ¢ = 2 and 4, the functions G 
do not have a superscript, as L,,, operates on the lattice distribution 
function in these instances which follows from the definition of the 
appropriate operators (62) and (62d). 

Finally, if in addition to non-equilibrium lattice vibrations we also 
have an external magnetic field as described in (a) of this section, the 
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Boltzmann equations (61) will acquire extra terms of the type (57) 
and result in conjugate sets of solutions f;,"* and f2,"*, do not affect the 
functions g;; and gs;, as mentioned before. The modified system of 
equations is then given by 


(fii**) + + (fii**) (66a) 


+ Linn? (G28) + Lm’ (66c) 


+ Linn’ (£25) (66d) 


mn 


with corresponding ones for a magnetic field in the opposite direction, 
which are again obtained from (66) by changing all + to — quantities. 
The solution of the various systems of equations with increasing com- 
plexity, (55, 60, 61 and 66) can also be obtained by a variational prin- 
ciple which will be described in Section III. 


III. VARIATIONAL PRINCIPLE 


1. Variational Principle for Distribution Function 


It has been shown by Kohler that the Boltzmann equation for the 
particle distribution function (cf. 27) is equivalent to a variational 
principle, which preserves the inherent symmetry properties. An ap- 
plication of the Ritz method enables one then to obtain a series solution 
for the distribution function and the transport coefficients in terms of 
infinite determinants, which can be evaluated to any required degree of 
accuracy. 

Consider an equation of the form 

L(F) =X (67) 
where L is an integral operator of the type (28) with positive definite 
kernel V(k, k’), and X a known function, such as given by (27). It 
can be shown that (67) can be replaced by a variational principle. 
For compactness of notation, let us introduce in addition to the inner 
product (F; G) defined previously (30) 


(F;G) = f (68) 


the product 


(F, X) = f FXdk. (69) 


atte 
= 
Xai 
= 
was a 
: 
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Then, if F is a solution of (67) and G any function which satisfies the 
relation 


(G;G) = (G, X) (70) 


then of all functions G satisfying (70), F is the one which makes (G; G) 
a maximum. The proof of this statement depends on the relation 
(F;G) = (G; F) and (G;G) > 0 which are consequences of the sym- 
metry property (20) and the fact that V(k, k’) is essentially positive 
definite. These relations have been derived in the last section, (29 and 
29’). To prove the variational principle we now proceed as follows. 
If we multiply (67) by G and integrate over k-space, we find 


(G; F) = (G,X) = (G;@) (71) 
on account of (70). From (29’) we have (F —G; F —G) >0. But 
0 < (F —G; F -—G) = (F; F) -— (F;G) — (G; F) + (G; 6) 

(F —G; F —G) = (F; F) — 2(G; F) + (G;G) — on account of (29) 


(F; F) — (G;G) on account of (29’) 
so that finally 
(F; F) > (G;G) (72) 


which proves our assertion. 
Thus our variational function is the inner product (F; F) subject 
to the subsidiary condition 


(F; F) = (F, X). (73) 
Consider now the new variational function 
H = (F; F) + X(F, X) (74) 


where \ is a Lagrangian multiplier. Variation of H7 with respect to 6F 
results in 


= f +X = 0 


which is equivalent to (67) provided we set \ = — 2. 
Applying these results to the system of equations (27) we introduce 
the variational functions H; (/ = 1, 2) (74) 


H, = fis) + X15) (74a) 
H, (fos fei) + fai Xj) (74D) 


which on variation with respect to f:; and fs; respectively, lead to the 
system of equations (27) provided we set 4; = A» = — 2 


: 
4 
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a. Solution of the Variational Principle 


The variational principle (74) is then solved by the Ritz method. 
In order to do this we expand f,,; and f.; in a power series in € (or 
x = e€ — ¢ which might be more convenient) 


fu = Da, xk, (75a) 


where the coefficients a, and 6, are to be determined by the Ritz method. 
From (68) we have 


fu) = (76a) 


fos) = (76d) 


where 


= (X°Rk;; L*k;) = 
while the subsidiary conditions (73) become 


> d,,'a,a, > a,‘a, 
r.@ r 


> = > 
where we have defined 
X15), (xrk,, X2;). 
The variational functions 7/7 (74a and 746) then become 


= a,a.d,,*? + a,'’d, (80a) 


Hi = + By (80b) 


Variation with respect to a, and }, then results in 


2 ¥ + = 0 


2 b.d,,"? + d2B,"? 


2 
J. 
: 
77 
(785) 
(79) 
r 
r 
2 
(81a) 
at 
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The Lagrangian multipliers \; and \,: are found to be A; = A: = — 2 by 
multiplying (81) by a, and b,, respectively and comparing with the 
subsidiary conditions (78). This, then, results in the following system 
of equations 


8," (826) 


for the unknowns a, and 6,._ The distribution function /;, and f.; can 
then be obtained in terms of infinite determinants. Let us denote by 
D the determinant whose elements are d,,, thus D = |d,,| and D,‘ and 
Ds‘ the determinants obtained from D by adding a column of a, and 8,, 
respectively, and a row consisting of x’k; the last diagonal element 
being zero. Then, it follows from (82) and (75) and the property of 
determinants that f;; and f2; are given by 


fis 
fri 


a, xk; = D,/D (84a) 
a, x°k; = D,/D. (84d) 


b. Entropy Production 


There exists a simple relationship between the variational function 
(F; F) and the rate of entropy production. The transport equation 
for the carrier distribution function can be written in the general form 
(cf. 13) 

Of /At = Of + Vf + e/hE- Vif. (85) 


Multiplying this equation by Inf and adding to it the corresponding 
equation, where f has been replaced by (1 — f), (85) can be written 
after integrating over dk 


+ div.s = —k f [In f — In (1 — (86) 


where S and s are the entropy density and entropy current density, 
respectively, and defined by 


S 


f (fing + ~ (87a) 


kf (ring + (1 — In (1 — f) (876) 
In obtaining (86) we have made use of the vector indentities 


div (Av) = v-VA +A divv. 


hor if 

; 
(88) | 

ay 
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If we now use the expansion for f given by (22), then the right-hand 
side of (86) can be written up to first orders in Q as 


~1/T f [O(k) — «/kTIL(Q)dk. (89) 


As the energy remains constant during the type of collisions considered 
here, the term proportional to the energy can be omitted, so that we 
have finally 


aS/at + divs = —1/T f Q(k)L(Q)dk = — 1/T(Q;Q) (90) 


which, apart from a factor (—1/7), is just the variational function 
(Q;Q) forthe problem. Moreover, as we can prove from the H-theorem 
that is always smaller than zero, it follows that (Q;Q) is positive 
definite, as shown earlier. For the steady state, of course, 0S/dt 
+ divs = 0 and hence the over-all change of f is also zero. 


c. Effect of Magnetic Field 


In the presence of an external magnetic field, it has been shown that 
the Boltzmann transport equation acquires an extra term which we 
represented by the operator L’ defined by (87), so that we can write the 
modified equation in operator form as 


X+ = L,(F+) + L,'(F*). (91) 


The entropy production due to magnetic scattering according to (90) 
is then given by 


lure. = 1/T FL'(F)dk = 0 (92) 


which is seen to vanish on account of (32), so that there is no explicit 
contribution from magnetic scattering to the production of entropy. 
The only effect of the magnetic field is to bend the electrons without 
changing the order of the distribution. 

On account of the symmetry relation (33) on reversal of the mag- 
netic field H, it is possible to obtain a variational principle for this case 
also, which will give rise to the correct transport equation (35), although 
the variational function does not contribute to the entropy production 
as shown above. Analogous to (74) let us define a variational func- 
tion /7 by 

H = (F-; F+) + (F-; Ft)’ + v(F-, X) (93) 


ge 
oe 
4 
ne 
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where (F-; Ft) is the inner product defined by (68) and the + and — 
signs indicate the appropriate solution F* depending on the direction 
of the field, and (F-;G*)’ denotes the inner product 


= [ = (G+; F-)’ (94) 


where the equality follows from (34). Hence, on variation this results in 


5H = (6F-; Ft) + (F-; + (6F-; Ft)’ 
+ (F-;6F+)’ + r(6F-, X) = 0. (95) 


Now, since L is independent of //, the first two terms in (95) are equal, 
while the next two are equal on account of (94), so that (95) will reduce 
to (94) provided we choose again \ = — 2. 

Applying this result to the system of equations (35) we now have 
for our new variational functions /7, (k = 1, 2) 


Ay = (fis fist) + 3 fast)’ — (96a) 
= 3 fei*) + (faim 3 X2;*). (96d) 


The solution is again obtained by the Ritz method, where now the 
functions f,,;+ and f2,;+ are expanded according to 


fit = Catek, fot = (97) 


and the new coefficients d,, are defined by 


Everything goes through as before, and we finally obtain two systems of 
equations of the type (82a and 826) depending whether the + or — 
sign has been used for the unknowns a and 6 (and coefficients d, a and 8). 


2. Generalized Variational Principle 


We have seen that the Boltzmann equation for the particle distribu- 
tion function is equivalent to a variational principle. It can be shown 
that in the presence of a lattice distribution function the resulting 
equations (47) can be complicated in this case due to the presence of 
additiona] integral operators and the lattice distribution function. 

Let Q and N be two functions satisfying the integral equations 


X = 1i(Q) + L2(N) 
Y = L,(Q) + L.(N) 


(99) 


ork. 
/ 
8 
a 
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where L, (i = 1 to 4) are the integral operators defined by (49). Fur- 
thermore, let R and M be two functions satisfying the relations 


(R, X) = (R; R)i + (R; M): 
(M, Y) (M;R); + (M; M), 


where the inner product (F; G),; has been defined by (51) and (F, X) by 
(69). As Q and N satisfy the equations (99) we have on multiplying 
by R and M, respectively, and integrating over k 


(R, X) = (R; + (R; 
(M, Y) (M;Q); + (M; N). 


from which it follows that the right-hand sides of the system of equations 
(100) and (101) are equal. 

Consider now the linear combination of expressions of the type 
(F;G); which can be formed for these functions 


(Fi F2; GiG2) = (Fi; + (Fi; Ge)2 + (Gis + (Gis Go) 


where F; and G; (j = 1,2) are any functions. From the symmetry 
relations (52) it follows that it is symmetric on interchange of the index 
1 with 2, that is: 


(100) 


(101) 


(102) 


GxG1) = (103) 


Moreover, on setting F; = F; = F and G,; = G; = G, the resulting ex- 
pression is positive definite, that is, 


(F;G) = (F; F)i + (F;G)2 + (G; F)s + (G;G), > 9. 


(104) 


The truth of this statement can be established by the use of the sym- 
metry relations and rearranging some of the terms. For example, due 
to the symmetry in k and k’ the first term in (104) can be written 
following the argument leading to (29’) as 


(F; F), = sf f V (k, k’)S(k, k’)[ F(k) — Paka’, 


As G is an odd function of its argument, the second term in (104) be- 
comes, on making use of the form of the operator L, (497), 


(F:G)s = eff V(k, k’)S(k, k’)[F(k) — Fk’) 1G (q)dkdk’ 


and so finally, the whole expression (104) is given by 


(F:G) = aff [F(k) + G(q) PP. (k, k’)dkdk’. (105) 


a 
2 
= 
| 
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This is clearly a positive definite function, as the probability P,(k, k’) 
= V(k, k’)S(k, k’) is of necessity a positive definite function and the 
integrand in (105) consists of a square multiplied by it. This, then 
completes the proof of our assertion. It might be pointed out that 
(105) is simply related to the total entropy production for this process. 
If we now replace the function F by Q — Rand G by N — M in (104) 
it can be shown that the following inequality holds 


(Q;Q)1 + (Q; + (N;3Q)3 + (NN), 
> (R;R): + (R; M)2 + (M;R)s + (M; M), (106) 


where we have made use of the relations (101) and (102) and the sym- 
metry properties (52). This enables us now to state the following 
theorem: 


If Q and N are two functions satisfying the integral equa- 
tions (99), and if R and M are two arbitrary functions satisfy- 
ing the relations (100), then of all these functions Q and N are 
the ones making (104) a maximum subject to the subsidiary 
conditions (101). Thus our variational function is the gen- 
eralized inner product (Q; NV) subject to the subsidiary con- 
ditions 

(Q, X) = (0; Qi + (Q; N)2 
(N, Y) = (N;Q)s + 


(107) 


This theorem can also be applied to the system of equations (47) 
simply by replacing Q by f,, (J = 1,2). The expression to be maxi- 
mized is then given by 


(firs £13) = (fiss + + fis)s + 
=1,2,3 (108) 


and the subsidiary conditions (107) replaced by 


(109a) 
¥15) = (giss + (21:3 


(fos, X2;) (fois + (fost 


(109d) 
(£2, = (goii fos)s + (Zoi; Loja 


From the form of these equations, it can be seen that we simply have 
two systems of equations (corresponding to the value of 7 = 1 or 2) 
for the functions f;;and g;;._ Consider now the new variational function 


H = (Q;N) +A(Q, X) + w(N, FY) (110) 


ah 
es 

and 
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where \ and yw are two Lagrangian multipliers. Independent variation 


with respect to Q and N yields the equations 


0) 
0) 


2L:(Q) + 2L2(N) + AX 
2L;(Q) + 2L,(N) + pY 


which are equivalent to (99) provided we set \ = wp = — 2. 
Applying these results to the system of equations (47) we introduce 
the variational functions 


Hy = B15) X15) ¥15) (111a) 
H, (fess Ao( foi, X2;) -+- (111d) 


which, on variation with respect to f;; and g,,;, and fe; and gi, respec- 
tively, lead to the system of equations (47) provided we set A; = Az 
Mi = Me = 


a. Solution of the Generalized Variational Principle 


The generalized variational principle (108) is again solved by an 
application of the Ritz method. For this purpose we expand the func- 
tions g,; and ¢», in a power series in y = hy/kT, in addition to the ex- 
pansion of the f’s carried out previously (cf. 75). 


= dy (112d) 


a 


It should be noted that there are three independent directions of the 
lattice waves in general, and it is therefore not permissible to omit the 
index j from the quantities c,;, and d;,, which are different for the various 
directions of the lattice polarization (j7 = 1 for longitudinal and 7 = 2, 3 
for transverse waves). wu; are the components of velocity of the lattice 
waves. The coefficients a, 6, ¢ and d are then determined by the Ritz 
variational principle subject to the subsidiary conditions (109). 

With the help of the expansions (75) and (112) the functions to be 


maximized can then be written as 


fess = (OD + 206 + did (113d) 


ake 
= 
\ 
4 
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where the coefficients and are given by 


(x°k;; = (114a) 


- 


The subsidiary conditions (109) can also be transformed into linear 
equations in the unknowns a, b, c, and d and yield 


r 


(115a) 


r r,@ 


r,8 


(1150) 


r,@ 


where, in addition to the a, and 8, defined previously (79), we also 
define 


The variational functions H/7; (J = 1, 2) (111) then become 


+r: + wi (1170) 


r r 


+ re + pe ;,. (117d) 


Variation with respect to the quantities a,, 6,, c;, and d;, then results 
again in a system of linear equations for these unknowns, where the 
Lagrangian multipliers \; and uw, (/ = 1, 2) are again found by multi- 
plying with a, etc., respectively, and comparing with the subsidiary 
conditions (109). The final system of equations is thus given by 


(fre. + — = Q (118) 


& 


: 
i 
rk 
201 | 
d : 
rs 

3 
2 

r = 

r = 
| 

=§ — 

r.8 

r r 
i 
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> + h,,**d j,) 6,"? . (118d) 


The distribution functions f; and g; can again be expressed in terms 
of infinite determinants. Let D denote‘ the determinant obtained from 
the matrix 


D,, = (119) 


and D,,, and D,s; be the determinants obtained from D by adding a 
column of ay or 86 and a row consisting of x’k; or y'u,; in the appro- 
priate places. The other elements and the last diagonal element are 
zero. Then it follows from (118), (75) and (112) and the property of 
adding determinants that the distribution functions f,; and g,; (2 = 1, 2) 
are given by 

= a, = (120a) 

fos b, Diss D (1206) 

= XU; = (120c) 


dj, Dyas, D. (120d) 


In the absence of a lattice distribution function (and the corresponding 
interactions characterized by the operators L; (¢ = 2, 3 and 4)) this, of 
course, reduces to the previous result (84). 


b. Entropy Production 


The connection between the variational function (Q; V) and entropy 
production can again easily be demonstrated. The transport equation 
for the phonon distribution function can be written in general form as 


(cf. 39) 
On/dt = ur Vr + On/dt. (121) 


Multiplying this equation by In m and subtracting from it the corre- 
sponding one for which has been replaced by » + 1 (in accordance 
with the Bose-Einstein statistics) we obtain after integrating over dq 


dS,/dt + divs, = —k f [In m — In (m + 1) ](On/0t).ou.dq (122) 


where S, and s, are the entropy density and entropy current density 
of the phonons, respectively, and defined by 


* To simplify the notation, the superscripts 1 and j are being omitted in (119) and similar 


expressions. 


202 
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S$, ft Inn — (nm + 1) In (m + 1) (123a) 
kf [nInm — (m + 1) In (m + 1) judg. (123d) 


For the corresponding entropy density and entropy current density of 
the carriers we again use the expressions (87), so that the total entropy 
density and entropy current density are given by their sums 


S=S.+ 5S, and Ss = 8. + 8;. (124) 


If we now use the expansions for f and n, (22) and (44) and limit our- 
selves again to terms linear in Q and N, both conservation laws can be 
written as 


aS/at + divs = — f [Q(k) — (e — + L.(N) 


rf [N(q) — hv + La(N) dq 


where we have made use of the operator equations for these functions (47) 


X = L,(Q) + L:(N) (125a) 
Y = L;(Q) + Li(N) 


to simplify the right-hand side. As the energy remains constant during 
the type of collision considered here, the terms proportional to « — ¢ 
and hy give no contribution to the entropy production, so that we have 
finally 

dS/dt + divs = — 1/T(Q; F) (126) 


which is (apart from the factor 1/7) the variational function (Q; F) 
(105) and thus demonstrates independently that (Q; F) is a positive 
definite function as shown earlier. 


c. Effect of the Magnetic Field 


As has been shown previously, the presence of the magnetic field 
introduces an extra term into the transport equations characterized by 
the operator L’ defined by (87) which obeys the symmetry relations 
(94). The appropriate transport equations in operator form (cf. 53) 
can be written as 


X+ = L,(Q+) + L.(N) + L'(Q+) (1272) 
= L,(Q+) + (1276) 


| 
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Analogous to (110) let us now define a variational function H by 


H = (Q-,Q+;N) + X) + w(N, Y) (128) 
where the inner product (Q-Q*; NV) is given by 


(Q-Q+; N) = (0-5 Q*)1 + + (N5Q*)s 

+ (N; N)a + + 3(Q*5Q°)’. (129) 
The first four terms do not involve the magnetic field, and hence it is 
immaterial whether we use the + or — sign in the functions Q. Fur- 
thermore, on setting Q- = Q*, the last two terms in (129) disappear on 
account of (32) and (128), and thus reduce to (105) in the absence of 
the magnetic field. 

Let Q*+ and N be two functions satisfying the integral equations 
(127) and Q- and N the ones obtained on reversing the direction of the 
field H. Then on forming the scalar products (Q, X) and (N, Y), one 
obtains the following subsidiary conditions 


(Q-, X) = (0°53 Q*)1 + N)2 + (1302) 
(N, Y) = (N;Qt)s + (N; (1308) 


Independent variation of 7 with respect to the functions 6Q~ and 6N 
then yields the equations 


(6Q0-;Q*), + (6Q-; +1 2(6Q-; Q+)’ 
+ 1/2(Q+: 6Q-)’ + x(80-, X) = 0 (1314) 


(Q-;5N)2 + (6N, Qt)s + (6N; + (N;36N), + Y) = 0 (1310) 


which are satisfied by (127) on eliminating the Lagrangian multipliers 
and uw through the subsidiary conditions (120), which for this case are 
seen to be equal to’ = uw» = — 1. (The difference from the previous 
result is due to the fact that we now consider Q- and Q* to be distinct 
functions.) 

This result can now be applied to the system of equations (53) and 
their conjugate ones by replacing Q* by f,:+ or fsit, and N by gi; or g»). 
The appropriate variational functions are given by 


iH, (fis, fist 1;*) (2145 ¥15) (132a) 
H, (fo » 3 225) (fos X2;*) rie (gai, (1326) 
The solution of the variational principle is again achieved by the Ritz 


method, where the functions f,;+ and f:;+ are expanded according to 
(97), and coefficients d, g and A are defined by (114), except that the 


+ 
ay 
a 
| 
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coefficients d,,* consist now of two terms according to (98) due to an 
additional term representing the magnetic field. The final system of 
equations is of the type (118), where, however, we have a different 
system of equations depending whether the + or — sign has been used 
for the quantities a+ and b* (and coefficients d, a+ and 8+), depending 
on the direction of the magnetic field. 


3. Variational Principle for Interband Scattering 


In the case of interband scattering, we have seen that the inner 
product (F; G) is defined as a sum of terms of the form (F"; G"), where 
m and n indicate the bands. If, in addition to that inner product 
(cf. 58 and 59) 


(F;G) = (F*;G) = > F"Linn(G")dK» (133) 


mn mn 


we also define the product 


m 


it is easily seen that the extremal function of the problem is given by 
(133) subject to the subsidiary condition 


As before, we can again define a new variational function by 


H= +A 5 (, X*) (136) 


min 


which, on variation with respect to 6F leads to the system of operator 
equations 


n 


where the Lagrangian multiplier is found to be equal to X = — 2 by 
comparison with the subsidiary conditions (135). Applying this to 
our system of equations (55) we again obtain variational functions //, 
and I, of the form (74), the only difference being that now each of the 
terms (f; f) and (f, x) appears as a sum over bands. 


a. Solution of the Variational Principle 


The variational principle can again be solved by the Ritz method. 
For this purpose we expand the functions f,," and f»," as an infinite 


m — 
m 
: 
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power series in €, (Or Xm) 


fim = en)" (138) 


= Ri Omn(€m)” (138d) 


where now the unknown coefficients a,,, and b,,, and the energy €, also 
depend on the band under consideration. 
With the help of the expansion for f the inner products (f1,; f1;) and 


(fei; become 
(fis firs) = XL (1392) 


mrn,s 


(fess fos) > (1396) 


mir 
where we have introduced the notation 

while the subsidiary conditions become 


mirn.s 


mrnie 


where we have defined 


Amr"?! (Riém’, X15") (142a) 
(Ri€m”, X2;"). (1426) 


The variational functions 7, (1 = 1, 2) are thus given by 


mr ne 


mr ne 
To find the values of the parameters we now vary these equations with 
respect to ad, and 6,,. From the conditions for an extreznum we obtain 
the system of equations 


2 + = 0 (144a) 


ne 


2 Gnd + = 0. (144d) 


UJ. F. I. 
r 
: 
mr 
mr 
| 
: 
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The unknown Lagrangian multipliers \; and \, are found on multiplying 
these equations by a,,, and b,,,, respectively, summing over m and n 


2 > > + = 0 


mrins 


2 > + Ae > = 


mrine mr 


and comparing with the subsidiary conditions (141). 

Thus from (144) we get two systems of linear equations for the un- 
knowns a,, and b,, which satisfy the stationary condition and the sub- 
sidiary conditions (141) 


We note that this result is analogous to our previous result (82) for the 
unknowns a, and b,, except that now we also have to carry out summa- 
tions over the various bands. The two distribution functions f;,; and fo, 
can again be obtained in terms of infinite determinants. Let D denote 
the determinant obtained from the matrix d»,... Where the elements are 
grouped according to increasing m (and m) fora given r (ands). Thus, 
for example, the determinant corresponding to m =n =r=s = 0,1 
is given by 

dov.00 doo.10 doo.o1 doo.11| 

| do1.00 doi,10 do1,01 | 
di1,10 


D= (146) 


Furthermore, let D,. and D,,3 be the determinants obtained from D by 

bordering it with a column of a,, and §8,,, respectively, and a row of 

k.(e)’ in the appropriate places, the others and the last diagonal ele- 

ment being zero. Thus, for example, Dos corresponding to (146) is 
given by 

| doo.10 doo.o1 doo.11 Boo| 


Dos 


00 Gia Bir 
kie,° 0 0 0 | 


Then, it follows that the functions f;;" and f2," are given by 


= — (148a) 
fer = — (1485) 


mr 
ns 
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b. Effect of a Magnetic Field 


In the presence of an external magnetic field the operator equations 
(55) have to be replaced by (60) which have an additional term L,,’(F™) 
describing the effect of the magnetic field. It has been noticed that 
this operator only has a free index describing the effect of the magnetic 
field on the band in question. Hence, let us define the product (F; F)’ 


(F; F)’ = (F*; F*)’ = (149) 


which, on account of the symmetry properties of L,,’, satisfies the relation 


(F-: Gt)’ = (F*-;G=+)’ = (G+; Fe-)’ = (Gt; F-)’. (150) 


In accordance with the procedure of Section 1, we take for our varia- 
tional function H7 the combination 


H = (F-; F+) + (F-; F+)' + X) (151) 


where (F-; F*) is the inner product defined by (133) and (F-; F+)’ the 
one by (149). The + and — signs again indicate the appropriate 
solutions depending on the direction of the field. Independent varia- 
tion with respect to F™ then results in the operator equation 


Xe+ = Lan(F**) + (152) 


for the distribution function, where the Lagrangian multiplier is found 
again to be \ = — 2 in order to satisfy the subsidiary conditions 


mr 


Applying this result to the system of equations (60), we now have 
for our new variational function /7, (/ = 1, 2) expressions of the form 
(96), the only difference being that now each of the inner products has 
to be summed over the various bands in accordance with the definitions 
(133, 134 and 149). The solution is again obtained by the Ritz method, 
where now the functions f;,"* and f:," are expanded according to 


= ky Amr*(€m)" (154a) 


= Ry Dmr*(€m)" (154d) 


r 


and the new coefficients d,,,.. are defined by 


= (Rien; Ree?) + (Riem’; (155) 


ye 
= m 
n 
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The rest of the calculation proceeds as before, and we finally obtain 
two systems of equations of the type (145) for the unknown coefficients 
a,,+ and 6,,*, respectively, depending whether the + or — sign has 
been used for the unknowns a and 3b and the coefficients d, a and 8. 
The distribution functions f can again be expressed in terms of infinite 
determinants, where again one has to differentiate between the + and — 
sign characteristic of the two solutions. 


c. Presence of Lattice Vibrations 

The carrier distribution function F™ depends on the band in ques- 
tion, while the lattice distribution function G does not. In view of this, 
it is necessary to differentiate between operators which operate on F™ 
and those which operate on G. and define the corresponding inner 
products accordingly, as can be seen from (64 and 65). Formally, 
however, we may expect that our previous results will still apply, pro- 
vided that we take account of the appropriate definitions of the inner 
products. 

Let us take then, for our variational function 


H = (Q;N) + X(Q, X) + w(N, Y) 
where the bracket (Q; N) denotes the linear combination 
Q;N) = (Q;Q)1 + + Q)s+(N; (157) 


whose components have been defined by (64). Furthermore, the 
products (Q, X) and (N, Y) are given, respectively, by 


(Q, X) = = f (158a) 
(N, Y) f NYdq, (1580) 
only the first of which contains the dependence on the band. \ and u 


are, of course, again Lagrangian multipliers. Independent variation 
of H with respect to Q™ and N yields the equations 


X Lan'(Q*) + = X* (159a) 


Lan?(Q*) + = Y (1595) 


where the Lagrangian multipliers \ and uw are found to be equal to 
= w = — 2 by comparison with the subsidiary conditions 
+ (Q; N).’ = (Q, X) (160a) 
(N;Q)3 + (N; = (N, Y). (1600) 
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Applying our result to the system of equations (61), we use the 
variational functions of the form (111), where the various terms are to 
be interpreted in terms of our definitions for inner products, such as 
(64) and (158). Variation with respect to f1," and g;; and f2;" and 
£2:, respectively, then leads to the system of equations (61), where the 
Lagrangian multipliers have been eliminated with the help of the sub- 
sidiary conditions (109). In order to solve the variational principle 
we expand the functions f;," and f2;" according to (138) and g,, and 
», according to (54). With the help of these expansions the functions 
to be maximized can then be written as . 


mr ne 


(fos 22;) > > +20 med ied je | (161d) 


mr ns 


where the coefficients d,,,,. are given by 


darned = €n'R;) mn! = (162a) 
(162c) 


The subsidiary conditions can also be transformed into linear equations 
in the unknowns and are given by 


mr nea 


ne 


mr 


mrne 


where the a,,, and 8,, are defined by (142) and y, and 6, by (116) in 
accordance with (158). 
The variational functions 7; (J = 1, 2) then become 


= (fics D amr me + jr (1642) 
(fois £25) + r2 + > jr. (1645) 


Variation with respect to the quantities @,,,, b,,, Cj, and d,;, then results 
again in a system of linear equations for these unknowns, which after 
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elimination of the Lagrangian multipliers finally results in the system 
of equations 


LL + je) (165) 


Bas + je) (165c) 


t+ js) — 0. (165d) 


m nes 


The distribution function f,;" and g,; (/ = 1, 2) can again be expressed 
in terms of infinite determinants. For this purpose, we define D to be 
the determinant obtained from the matrix 


Zmr,ne 


where the elements have been ordered according to the convention 
leading to (146). Furthermore, we define D,n«, and D,ngs to be the 
determinants obtained from D by adding a column of a7, Or Bd, and a 
row consisting of ¢,,’k; in the appropriate places, the others and the last 
diagonal element being zero. Also let D,., and D,ss be again the de- 
terminants obtained from D by bordering it with a column of a@»,y, or 
Bnd, and a row of y"u,, the last diagonal element being zero. Then it 
follows from the properties of determinants that the distribution func- 
tions f” and g are given by 


= = Dimgs/D (167b) 


gy = Ldj yu; = (167d) 


If, in addition to the lattice distribution function g, we also have an 
external magnetic field, the appropriate variational principle is given by 
one of the form (132), where the various terms have to be interpreted in 
accordance with previous results. For example, the products (f1,~, x,;*) 
and (f2;~, X2;+) are to be considered as sums over the various bands as 
in (134) with the + and — signs indicating the direction of the field. 
The inner product (f-f+;g) is of the form (129) where, in addition to 
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the terms defined by (122), we also have to include those characteristic 
of the magnetic field and given by (149). Independent variation with 
respect to f,;"* and g,,; then leads to the appropriate equations. The 
solution of this more complicated variational principle is again achieved 
by the Ritz method, where the functions f,"* and f2."* are expanded 
according to (154) and g; and g. according to (112). The coefficients 
d, g are defined as in (162) except that the coefficients d,,,,,, Consist now 
of two terms and are given by (155). The final system of equations is 
of the type (165), where, however, we now have two sets of equations 
depending on whether the + or — sign has been used for the quantities 


Am,* and b,,,* (and coefficients) depending on the direction of the mag- 


netic field. 
This, then, concludes our discussion of the variational principle as 


applied to transport theory, where we have tried to discuss the various 
cases arising whether (7) an external magnetic field, (77) non-equilibrium 
lattice distributions or (777) interband scattering are present, singly or 
in combination. In the next section we shall apply our result to the 
evaluation of the various transport coefficients. 


IV. TRANSPORT PHENOMENA 


1. General Formulation of Transport Coefficients 

The various transport coefficients can be obtained from the general 
expressions for the electric current density J and energy current density 
W. The former density is given by 


(i = 1, 2, 3) (168) 


while the latter consists of two contributions, one due to the electronic 
energy density, and one due to the lattice current density, and is given by 


W, = (1/82?) + (1/87?) f uhvndq (¢ = 1,2,3). (169) 


Here f and » denote again the carrier and lattice distribution functions 
for which we use their respective expansions (22) and (44). The in- 
tegrals containing f, or n, will give no contributions to the current and 
energy current densities, while the terms containing Q and N can be 
separated into two parts according to (27) and (48). Of course, in the 
simple case where we neglect the deviation of the lattice distribution 
function from equilibrium, the latter will yield no contribution. Using 
these expansions we obtain for J; and W,; (where 7 = 1 to 3 denote the 
components of J and W) 

J, = [Ki + + Ki; (e/T)0T/dx‘} (170) 


hve: 
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respectively, where we have defined the energy integrals 


= f = Guy 


82°K 


= (fei, X11) 


(fis, X24) + (213, Yai) 


= ff — f 


and where we have also made use of the definition (69). 


(fei, + (£25, Yo) 


Hence it is 


(172a) 


(1726) 


(172c) 


(172d) 


seen that the energy integrals (172) are nothing else than the subsidiary 
The K’s defined above are not independent, but 


conditions (109). 


K,; = 


satisfy the so-called Kelvin relations 


K,;® = K ;;, K,;™ = K;;™. 


The proof of (173) can be carried out with the help of the symmetry 
relations (52) and the operator equations (47). 


(173) 


From the fact that the electric current density J; can be written as 


resistivity p,; is given by 


(pi J; — — (1/e)dt/dx' 


p 


J, = 6,;E; 


we deduce from (168) that the conductivity tensor o,,; is given by 


05; = eK 


= (1/e) 


— «;,0T/dx') — (¢/T)J; 


(174) 


(175) 


If we denote the reciprocal tensor of K;;“ by K,;“, then the electrical 


(176) 


Equations (170) and (171) can be solved for the electric field and the 
heat current density 


213 
j 
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where .S;; is the tensor of the Seebeck coefficient, I1,; that of the Peltier 
coefficient, and «;; the heat conductivity tensor. Comparisons of (177) 


and (178) with (170) and (171) yields for these tensors 


= + + Kin (179) 


= — — Kans (180) 


= — Kan Kaj” (181) 


where 6,;; denotes the Kroenicker delta, 6;; = 1 if 7 = j, 6,; = 0 if 
t ~ j. Moreover, from the Kelvin relations (173), it also follows that 


Il, = — TS;; (182) 


in accordance with known results. 

The expressions for the various transport coefficients (175, 176, 
179-182) are quite general and do not depend on any particular model 
used. Their value depends only on the energy integrals which we shall 
evaluate for the various cases of interest by the Ritz variational 
principle. 


2. Effect of Magnetic Field 


We have seen that in the presence of an external magnetic field the 
transport equations acquire an extra term due to the interaction with 
the magnetic field, and that the distribution functions f; and f, have 
to be divided into two conjugate sets f;+ and f,- (J = 1, 2) correspond- 
ing to the direction of the field. The effect of this separation will be 
felt in the energy integrals which are now a function of the magnetic 
field H. The electric current and heat current densities are still given 
by their respective expressions (170) and (171), but the energy integrals 
K,; are now defined by 


82°K,; (H) (fis-, (183a) 
(H) = (1838) 
(H) = (fe; ‘ X2;*) + (g2j, (183d) 
where the products (f-, x*) are of the type (69) with the appropriate 


direction of the magnetic field H. From (183) we can obtain the con- 
jugate set corresponding to —H by changing all + to — signs and vice 
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versa. It can be seen that the energy integrals (183) are again the 
subsidiary conditions (130) which differ from the previous ones (109) 
by having an extra term of the form (f;-; f:*+)’ and two sets of the dis- 
tribution function f; (/ = 1, 2). 

With the help of the symmetry relations (108) and (94) for inner 
products involving the operator L’ (31) it can be shown that they satisfy 
the Kelvin relations (173), where now, however, the change of direction 
of the magnetic field must be taken into account. We thus obtain 


K,,;“(H) = (184) 


The various transport coefficients are still given by their appropriate 
expressions, (175, 176, 179-181), but where now the energy integrals 
defined by (183) must be used. From their symmetry relations under 
reversal of the magnetic field (184), one sees readily that the following 
relations hold among the various transport quantities 


o;;(H) K,;(H) K;,(—H) (185) 
TS,,(H) = — I,.(—H). 

If one separates the above tensors into their symmetric and anti- 
symmetric parts, one finds that the symmetric parts of these tensors are 
even functions of the magnetic field, while their antisymmetric parts are 
odd functions of the field, that is, they change sign on reversal of the field. 


3. Calculation of Energy Integrals 


We have seen in the previous sections that all thermoelectric coeffi- 
cients can be expressed in terms of the energy integrals K,;._ Thus, te 
obtain expressions for these integrals, it is only necessary to calculato 
the various thermoelectric and magnetic coefficients. By applying an 
elegant method due to Enskog to the various cases, these can be ob- 
tained in terms of infinite determinants without the explicit solution of 
the transport equations. 


a. Carrier Distribution Function 


Let us first consider the simple case for which the deviation of the 
lattice distribution function from equilibrium has been neglected, that 
is, ma = n,. The appropriate energy integrals are then given by 


= f = (fr;, (1862) 


= f = (fas. X10) (1860) 


; 
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= J = (186c) 


= J = (fig (186d) 


where the distribution functions f;; and f2; are given by (184) and x, 
and x2; are defined by (27). From these definitions and the ones for 
a, and 8, (79), it is seen that the various energy integrals are just sums 
of products of the coefficients a, and 6, and the integrals a, and §8,. 
Thus, if we denote by D,s the determinant obtained from D = !d,,| by 
adding a column of 8, and a row of a, (with similar additions for the 
other types), we obtain for the energy integrals (186) 


82°K,;) = > — Daa/D (187a) 


r 


= = — Das/D (187d) 


= a6," Dsa/D (187c) 


r 


82°K, = > b,B,*? Da3/D. (187d) 


The addition of an external magnetic field leaves that result unchanged 
in principle, except that the matrix elements d,, (77) have to be replaced 
by the corresponding ones including the effect of the field (98), and 
+ or — signs have to be used for the coefficients a,+ and 8,+ depending 
on the direction of the field. In other words, we will now have two 
sets of energy integrals K,;(H) and K,;(—H) of the type 


82°K,7(H) = — D.-s+/D; &2°K,;°(—H) = — D,ats-/D (188) 
and similar ones for the other cases in (187). 


b. Effect of Lattice Vibrations 


In the general case where the deviation of the lattice distribution 
function from equilibrium is also taken into account, the energy in- 
tegrals K,; are given by (172), the distribution functions f;; and gi; 
(1 = 1,2) by (120), and the coefficients a,, B,, y, and 6, by (79) and 
(116), respectively. It can be seen again that the energy integrals are 
just sums of products of the coefficients a@,, b,, c, and integrals a,, 8,, Y, 
and 6,, and therefore can be expressed again in terms of infinite de- 
terminants. With the help of the defining equations and the properties 


Nine 
: 
r 
Ds 
r 
| 
: 


Sept. 1959.] VARIATIONAL PRINCIPLE IN TRANSPORT THEORY 217 


of adding of determinants we obtain for the energy integrals A,; in 
this case 


;! — Darye,/D (189a) 


II 
ll 


+ 


— Dsiay/D (1896) 


r 


— Daysi/D (189¢) 


82°K ;;* (a,B,"? C 


D3s33/D (189d) 


where D,,; is the determinant obtained from D by a line consisting of 
ay and a column of 86 with similar definitions for the other determinants. 
(It should be noted that in this case D is the determinant obtained from 
the matrix D,, defined by (119).) 

As in the simple case discussed above, the addition of an external 
magnetic field provides only minor modifications of that result. For 
the matrix elements d,, we have to use the ones which include the effect 
of the field (98), and + and — signs have to be used for the quantities 
a,* and 6,*+, depending on the direction of the field. Of course, the 
integrals y,+ and 6,* are not affected by the field and remain unchanged. 
The final expressions for the energy integrals are still of the form (189), 
with the modification introduced by (88) resulting in two sets of 
integrals K;;(H) and K,;(—H) as defined by (183), from which the 
various transport quantities are to be calculated in accordance with 
(175, 176, 179-181). 


4. Transport Coefficients for Interband Scattering 


a. Without Taking Account of Lattice Vibrations 


In the case of interband scattering the electric and energy current 
densities may be considered to be given as sums of the appropriate 
expressions for each band. Thus 


(190) 


‘where 
+ (191) 


and 


wet. (192) 
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with 


Wim = (1/82?) f = — + T8/ax,(¢/T)] 
+ (193) 


The K’s are the energy integrals for the various bands and analogous to 
(186) defined by 


= (fo;", X10") (194) 
= x2") (194c) 


= (f2;™, X2;"). (194d) 


As before, the indices m and n refer to the bands and 7 and j to the 
direction. From (190-193) it can be seen that the current and energy 
current densities can be written as 


Ji = + eT(0/dx;)(¢/T)] + Kiz(e/T)8T/dx;} (195) 


and 


= — + T(0/dx,)(¢/T)] + Kis(1/T)0T/dx;}, (196) 


respectively, where we introduced the notation 


K;;* = > (u = 1, 2,3 or 4). (197) 


The expressions (195) and (196) are, therefore, of the same form as the 
corresponding ones in the standard theory (170) and (171), the only 
difference being the definition (197) for the energy integrals. From 
(195) and (196), all transport coefficients can now be obtained in the 
usual manner (see Sections 1 and 2). 

The energy integrals can be obtained again in terms of infinite 
determinants by applying the results of Section III to this case. It 
has been shown there that the distribution functions f;," and f2," are 
given by (148), and the various K’s are then obtained by forming the 
inner products of the f’s with x; or x: in accordance with (197). Thus, 
denote by Das’? the determinant obtained from D,,3‘? by replacing the 
row of ye," by am, in the appropriate places, the others and the last 
diagonal element being again zero. In other words, Dyas‘? is of the 
same form as Ds except that it has more rows and columns and some 
zeros in the last row. For example, D.a3'’ obtained from D,3‘/ (147) is 
given by 
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| doo, 00 doo, 10 doo, ol doo, ll Bus 


(198) 
dii,01 Bi 
Qo1 0 0 | 
We then obtain for the various K,;,;“" 


LX = (1994) 
= / (199d) 


Dirge D (199¢) 


r 


Summation over m then yields the A,," (uw = 1, 2, 3 or 4) (197). In- 
spection of the determinants D,,.3 such as (198) shows that the summa- 
tion over m can be carried out directly in the determinants. Hence, if 
we define D,s to be a determinant obtained from D by bordering it with 
a row of a», and a column of 8,, (defined by (142)), where the last 
diagonal element is again zero, the K,;“ are then given by 


= — Daa/D (200a) 
= — Das/D (2000) 
= — Ds,/D (200c) 
= — Dss/D (200d) 


in accordance with our previous results (187), the only difference being 
that now we have more rows and columns corresponding to the different 
values of (mr, ns). The various transport coefficients are then obtained 
from (200) and the appropriate formulae. 

The effect of an external magnetic field is taken into account by 
replacing the matrix elements d,,;,,. by the ones including the effect of 
the field (155), and supplying + or — signs to the integrals a,,,+ and 
Bm-* depending on the direction of the field. As before, we will again 
have two sets of energy integrals K,;“(H) and K,;“(—H) of the type 
(188), the only difference being that each set is of the form (200). 


b. Presence of Lattice Vibrations 


In the case for which deviations from equilibrium of the lattice 
distribution function are also taken into account, the energy integrals 
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K,,“ for the individual bands are given analogous to (172) by 


X15") 


= (2014) 


= (f2;", X15") (2016) 
= (f15", + (215, Vos) (201c) 
(f2;", + (g2j, V2). (201d) 


The distribution functions f;," and g,;" (J = 1, 2) are given by (167), 
where the appropriate integrals and matrix elements have been defined 
by (116) and (142). The various K’s are then obtained from (167) 
by forming the inner products of sums of f’s and g’s with x and y. 
From (201), the various defining equations and the properties of addi- 
tion of determinants, we then obtain for the energy integrals in this case 


(202a) 


D 


,'™ > + C jr¥r'?) 


r 


— Dnstay/D (2026) 


r 


(OmrBmr'? + d ;,5,*?) (202d) 


where the determinant D,,3:2, is the sum of the determinants Dina, and 
Day replacing, however, the row of €,"k; by Bar Or an, and the row of 
y'u, by y, or 6, in the appropriate places, the others and the last diagonal 
element being zero. Summation over m then yields the K,;" for this 
case. It can be seen that this summation can be carried out again 
directly in the determinants. Thus, if we denote by Days; the de- 
terminant obtained from D (as given in (166)) by bordering it with 
a row of a,,,y, and a column of 8,,y,, where the last diagonal element is 
again zero, the K,;“ are given by the previous expressions (189), the 
difference being that we now have more rows and columns corresponding 
to the different values of (mr, ns) and that the various matrix elements, 
Amr, nsmr.ne are given by (162) and the integrals Bur, 5+ by (116) 
and (142). 

The modification introduced by an external magnetic field results 
again in two sets of energy integrals K,;“(H) and K,,;“(—H) with the 
matrix elements dm, a. given by (155) and + or — signs introduced in 
the quantities a,,+ and 8,,,+ depending on the direction of the field. 
From these expressions for the energy integrals, the various transport 
quantities can again be calculated in accordance with their defining 
equations (175, 176, 179-181) as before. 
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We have thus been able to demonstrate how a variational principle 
can be formulated which is equivalent to the transport equations for 
carrier distribution functions, lattice vibrations, external magnetic 
field, interband scattering and combinations of all or any of these. 
General expressions for the various transport coefficients have been 
given in terms of energy integrals, which have been evaluated in terms 
of infinite determinants for any of the above cases. To evaluate the 
different determinants the coefficients d, etc. must be calculated for 
which purpose it is necessary to make definite assumptions abcut the 
form of the interaction and energy dependence. Work along these lines 
for the calculation of the transport coefficients for warped energy sur- 
faces is in progress and has been reported elsewhere (12). In this paper 
we have only been interested to give a systematic discussion of the vari- 
ational principle as applied to transport theory and its application to the 
calculation of the transport coefficients. 


REFERENCES 


(1) P. M. Morse Ano H. Fescupacn, “Methods of Theoretical Physics,"" New York, McGraw- 
Hill Book Co., Inc., 1953, Part II, Chapter 9. 
(2) A. H. Witson, ‘The Theory of Metals,” New York, Cambridge University Press, 1953. 
(3) M. XonLeR, Z. Physik, Vol. 124, p. 772 (1948) ; ibid., Vol. 125, p. 679 (1949). 
(4) J. M. Ziman, Can. J. Phys., Vol. 34, p. 1256 (1956). 
(5) L. ONsaGErR, Phys. Rev., Vol. 37, p. 405 (1931); ibid., Vol. 38, p. 2265 (1931). 
(6) E. H. SonpHermmer, Proc. Roy. Soc. A, Vol. 234, p. 391 (1956). 
(7) D. ENskoG, dissertation, University of Utrecht, 1917. 
(8) D. Dorn, Zeits. f. Naturf., Vol. 12a, p. 739 (1957). 
(9) M. Kon_er, Ann. Physik, Vol. 40, p. 165 (1941). 
(10) F. G. MOLINER AND S. Stmons, Proc. Camb. Phil. Soc., Vol. 53, p. 848 (1957). 
(11) G. E. Tauser, Can. J. Phys., Vol. 36, p. 1308 (1958). 
(12) G. E. TauBER, AND L. Sorrer, J. Phys. Chem. Solids, Vol. 8, 138 (1959). 


3 


THE FRANKLIN INSTITUTE 


LIBRARY 


The Committee on Library desires to add to the collections any technical works that 
members may wish to contribute. Contributions will be gratefully acknowledged and placed 
in the Library. Duplicates received will be transferred to other libraries as gifts of the donor. 


307 items have been added the past month. 
Photostat service. Photostat prints of any material in the collections can be supplied 


on request. 
The Library and reading room are open from 9 A.M. until 5 p.m. on Mondays, Tuesdays, 


Thursdays, Fridays, and Saturdays; and 2 P.M. until 10 p.m. on Wednesdays. 


RECENT ADDITIONS 
AERONAUTICS 


AGARD. Combustion and Propulsion. 1958. 
Fart, J., Ep. Air Intake Problems in Supersonic Propulsion. 1958. 


AGRICULTURE 


The Scientific Principles of Crop Protection. Ed. 4. 1959. 


MARTIN, HUBERT. 


ASTRONOMY 


MEHLIN, THEODORE GREFE. Astronomy. 1959. 


BIOLOGICAL CHEMISTRY 
A Report on Biochemical Problems and on 


Havurowit1z, FEvIx. 
Biochemical Research Since 1949. 


Progress in Biochemistry. 
1959. 


CHEMISTRY AND CHEMICAL TECHNOLOGY 


AstLeE, MELVIN J. AND SHELTON, J. Retp. Organic Chemistry. Ed. 2. 1959. 

BERGLER, WALTER H., JR., ET AL. Structural Adhesives. 1958. 

BJERRUM, JANNIK, ET AL. Stability Constants of Metal-ion Complexes, with Solubility 
Products of Inorganic Substances. 1957. 

LANGFORD, KENNETH E. Analysis of Electroplating and Related Solutions. Ed. 2. 1958. 

Lykow, A. Transporterscheinungen in Kapillarporésen Kérpern. 1958. 

Meats, R. N. anp Lewis, F. M. Silicones. 1959. 

Noyes, ARTHUR AMOS AND MULLIKEN, SAMUEL P. Laboratory Experiments on the Class 
Reactions and Identification of Organic Substances. Ed. 3. 1915. 

Otsen, J. C. A Text-Book of Quantitative Chemical Analysis by Gravimetric, Electrolytic, 
Volumetric and Gasometric Methods. Ed. 2. 1905. 

Pattison, E. Scort, Ep. Industrial Fatty Acids and their Applications. 1959. 

Wueat-ey, P. J. The Determination of Molecular Structure. 1959. 

YounG, Jay A. Practice in Thinking ; a Laboratory Course in Introductory Chemistry. 1958. 


CRYSTALLOGRAPHY 


Vector Space and its Application to Crystal Structure Investi- 


BUERGER, MARTIN JULIAN. 
gation. 1959. 


222 


: 
: 


Sept., 1959.] THE FRANKLIN INSTITUTE 223 


ELECTRICITY AND ELECTRICAL ENGINEERING 


Amos, S. W. AND BirkinsHaw, D.C. Television Engineering. 1958. 

Davipson, I. A. Receiving Aerial Systems for Broadcast and Television. 1957. 

Donapson, P. E. K., ET AL. Electronic Apparatus for Biological Research. 1958. Electric 
Motors and Generators. 1959. 

Frost-Smitu, E.H. The Theory and Design of Magnetic Amplifiers. 1958. 

Gitte, ANGELO C. Transistors. 1959. 

LANGsporF, ALEXANDER Suss. Principles of Direct-Current Mechanics. Ed. 6. 1959. 

LoverinGc, W. F. Radio Communication. 1958. 

PumpHREY, Frep H. Fundamentals of Electrical Engineering. Ed. 2. 1959. 

SESHU, SUNDARAM AND BALABANIAN, NORMAN. Linear Network Analysis. 1959. 


ENGINEERING 


Bors, Sipney F. AND GENNARO, JosEPH J. Advanced Structural Analysis. 1959. 
CuHronowicz, ALBIN. The Design of Shells. 1959. 

Donkin, C. T. B. Elementary Practical Hydraulics of Flow Pipes. 1959. 
Potter, Puitie J. Power Plant Theory and Design. Ed. 2. 1959. 


GENERAL 

DipNeER, BERN. The Atlantic Cable. 1959. 

JunGck, Ropert. Brighter than a Thousand Suns. A Personal History of the Atomic 
Scientists. 1956. 

O'Dea, WittiaM T. The Social History of Lighting. 1958. 

VON FANGE, EUGENE K. Professional Creativity. 1959. 


LOGIC 


BARTLETT, FREDERIC. Thinking; an Experimental and Social Study. 1958. 
BERKELEY, EpMuND CALLis. Symbolic Logic and Intelligent Machines. 1959. 
BRIDGMAN, Percy WiLuiaMs. The Way Things Are. 1959. 

CarnaP, RupoiFr. Introduction to Symbolic Logic and its Applications. 1958. 
SCHRGODINGER, ERwin. Mind and Matter. 1958. 


MANUFACTURE 


BosTROM, SIEGFRIED. Kautschuk-Handbuch. Vol. 3. 1958. 
RopGers, THomas M. Handbook of Practical Electroplating. 1959. 


MATHEMATICS 


ALLENDOERFER, Cart B. AND OAKLEY, CLETUS O. Fundamentals of Freshman Mathematics. 
1959. 

BATEMAN, Harry. Partial Differential Equations of Mathematical Physics. 1959. 

James, GLENN AND JAMES, Rospert C. Mathematics Dictionary. Multilingual Ed. 1959. 

JEENEL, JOACHION. Programming for Digital Computers. 1959. 

LANG, SERGE. Abelian Varieties. 1959. 

McCormick, Epwarp Mack. Digital Computer Primer. 1959. 

Rick, HAROLD S. AND KNIGHT, RayMOND M._ Technical Calculus and Analysis. 1959. 

Taytor, ANGus Etuis. Calculus, with Analytical Geometry. 1959. 

ZYGMUND, ANTONI. Trigonometric Series. 1959. 


MECHANICAL DRAWING 


Biack, Eart D. Graphical Communication; Drafting, Sketching, and Blue-Print Reading. 
1959. 
LuzADDER, WARREN JAcoB. Fundamentals of Engineering Drawing for Technical Students 


and Professional Draftsmen. Ed. 4. 1959. 


| 
: 


224 THe FRANKLIN INSTITUTE 


MECHANICAL ENGINEERING 
Sanpers, T. H. Springs and Suspension. 1930. 


MEDICINE 


BEHRENS, CHARLES FREDERICK. Atomic Medicine. Ed. 3. 1959. 


METALLURGY 


Advances in Steel Technology in 1956. 1958. 
BETTERIDGE, W. The Nimonic Alloys. 1959. 


METEOROLOGY 


DeFAnt, ALBERT AND DEFANT, FrrepRricH. Physikalische Dynamik der Atmosphare. 1958. 
Massey, Harrie STEWART WILSON AND Boyp, R. L. F. The Upper Atmosphere. 1959. 


MINERALOGY 


Kraus, Epwarp Henry; Hunt, WALTER FRED AND RAMSDELL, LEWIs STEPHEN. Mineralogy. 
Ed. 5. 1959. 


MISSILES AND ROCKETS 
Rocket Encyclopedia. 1959. 
Space Weapons; a Handbook of Military Astronautics. 1959. 
ULANoFF, STANLEY. Illustrated Guide to U. S. Missiles and Rockets. 1959. 


NUCLEAR ENGINEERING 


GoLpsTEIN, HERBERT. Fundamental Aspects of Reactor Shielding. 1959. 
U. S. Atomic EnerGy Commission. Progress in Peaceful Uses of Atomic Energy; July 
December, 1957. 1958. 


PETROLEUM TECHNOLOGY 


Crart, BENJAMIN CoLE AND Hawkins, M. F. Applied Petroleum Reservoir Engineering. 
1959. 


PHOTOGRAPHY 
WHEELER, Leste J. Principles of Cinematography. Ed. 2. 1958. 


PHYSICS 


CarLos G., JR. Liquid Scintillation Counting. 1958. 
BoGo.iuBsoy, N. N.; ToLMacueyv, V. V. AND Suirkov, D. V. A New Method in the Theory 
of Superconductivity. 1959 
Hayes, WALLACE D. AND PrRoBsTEIN, RONALD F. Hypersonic Flow Theory. 1959. 
Hooker, K. H. aNnp WEIMER, K., ED. Lexikon der Kern-und Reaktortechnik. Vol. 1 and 2. 
1959. 
Horcukiss, HoMER J. AND FarrBANKs, FLoyp C. Principles of Experimental Physics. 1913. 
MArsHAK, RoBert E. Meson Physics. 1952. 
Scott, WiLL1AM TaussiG. The Physics of Electricity and Magnetism. 1959. 
SHORTLEY, GEORGE H1RAM AND WILLIAMS, DupLey. Principles of College Physics. 1959. 
UNIVERSITE DE GRENOBLE. The Many Body Problems. 1959. 
voN ARTHUR Molecular Science and Molecular Engineering. 1959. 
WIGNER, EUGENE P. Group Theory and its Application to the Quantum Mechanics of 
Atomic Spectra. 1959. 
RESEARCH 


AMERICAN CHEMICAL Society. COMMITTEE ON PROFESSIONAL TRAINING. Directory of 
Graduate Research. 1957. 1958. 


: 
: 
VJ. F. I. 
~ 


Sept., 1959.] THE FRANKLIN INSTITUTE 225 


NATIONAL SCIENCE FounpaTION. Scientific Research and Development in Colleges and 
Universities ; Expenditures and Manpower, 1953-54. 1958. 
SANITARY ENGINEERING 

SADDINGTON, K. AND TEMPLETON, W. L. Disposals of Radioactive Waste. 1958. 

SCIENCE 
McCrensky, Epwarp. Scientific Manpower in Europe. 1958. 
Roucek, JosepH S., ED. The Challenge of Science Education. 1959. 
THIRRING, HANS. Energy for Man; Windmills to Nuclear Power. 1958. 

TEXTILES 
HowEL.t, H. G.; Mieszxkis, K. W. AND TaBor, D. Friction in Textiles. 1959. 


YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 135-yvear history, there is vital need for 
The Franklin Institute to effectively promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 
We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute's educational programs are impressive, for they begin with students in 
the early grades of our elementary schools and continue throughout an individual's 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased. 

The Franklin Institute is not richly endowed. It is a non-profit organization, 
depending for encouragement and support on an understanding public. Capable 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology through education. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 
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MATTERS PERTAINING TO ELIMINATION OF BACKGROUND IN 
LOW ENERGY NUCLEAR PHYSICS MEASUREMENTS, 
WITH SPECIAL REFERENCE TO THE BARTOL 
ONR GENERATOR*+ 


BY 
C. P. SWANN, V. K. RASMUSSEN AND H. O. ALBRECHT 


ABSTRACT 


A new low background target room has been constructed for the Bartol-ONR 
Van de Graaff generator. This room is constructed of low radioactivity content con- 
crete and serves to shield the target area from the radiation arising from the soil, the 
beam analyzing magnet, and the small Bartol generator. Over the energy range from 
400 kev to 3 Mev, the natural radioactive background is at least a factor of 5 lower 
than in the original target room. The effects of neutrons from the analyzing magnets 
and the small machine have been reduced to a negligible amount. The cost of the 
special concrete used in this structure was less than twice the cost of normal concrete. 


INTRODUCTION 
Certain problems involving the reduction of natural radioactive 
background and background from other accelerators have arisen in 
connection with the Bartol-ONR Van de Graaff generator. While these 
are hardly unique problems, we found that the published literature did 
not seem to contain much that was directly applicable to our situation.! 
The following is written, therefore, with the hope that a somewhat 
detailed description of our solution of these problems might be useful 
to others in the field of low-energy nuclear physics. 


GENERAL NATURE OF THE PROBLEM AS APPLIED TO THE BARTOL-ONR ACCELERATOR 


The Bartol-ONR 5.5 million volt Van de Graaff generator built at 
the Bartol Research Foundation of The Franklin Institute has now been 
in operation for seven years. During this period, there has been a 
continuing program to reduce the y-ray and neutron background result- 
ing from the accelerated particles striking parts of the vacuum plumbing 
other than the intended target. This program, initiated as part of the 
nuclear resonance fluorescence studies, has been successful to the extent 
that in a typical resonance fluorescence measurement the major part of 
the background was independent of whether or not the accelerator was 
running; that is, this background resulted from the natural radioactive 


* Assisted by the joint program of the Office of Naval Research and the United States 
Atomic Energy Commission. 

t The authors wish to express their appreciation of the support and encouragement given to 
this project by the retiring Director of the Bartol Research Foundation, Dr. W. F. G. Swann, 
and of the advice and assistance by our colleague, Dr. F. R. Metzger. 

' The authors are indebted to Dr. H. E. Gove for a helpful discussion of the solution of 
similar background problems at Chalk River. 
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contamination of the building and of the underlying soil. The various 
possibilities for absorbing out this radiation were considered, but the 
technical problems and cost seemed far too great. 

The building housing the Bartol-ONR accelerator also houses the 
1.7 million volt Bartol Van de Graaff, and the target areas for the two 
machines are separated only by a 2-ft. thick concrete block wall (see 
Fig. 1). The smaller machine is used as a neutron source, and when it 
is in operation experiments with the larger machine involving detection 
of either neutrons or y-rays are interfered with severely. Induced 
radioactivities may also be troublesome for several hours after the 
smaller machine is turned off. Serious limitations have thus been 
imposed on the operating time of both machines. Means by which the 
shielding between the two target areas could be improved were con- 
sidered, but even the best, which would have been very expensive, 
would not have reduced the neutron effect to an amount small com- 
pared with the natural room background. 

The use of beams other than protons in the Bartol-ONR machine 
was also banned for two reasons. First, the vacuum system becomes 
contaminated with the gas used in the ion source, and it may take 
many hours of operation with hydrogen before the system is cleaned 
out. Second, the radioactivity induced by these other beams can be 
appreciable and may take days to decay to a low level. 

A solution to all of the problems appeared to be the construction of 
a new target room for the large generator. This room had to be ade- 
quately shielded against the neutron effects of the small machine and the 
radioactivity of the ground and had to be constructed of low Z, low radio- 
activity content materials. The low Z was desirable to reduce the 
intensity of elastically scattered gamma rays which might add appre- 
ciably to the background of the resonance fluorescence studies. 

Only proton beams should be used in this room while the old room 
could be used for other particles as well as for all charged particle 
spectrometer studies. 

In order to determine the feasibility of meeting the requirements for 
this new room, several studies were conducted. First, the radioactivity 
within the soil around the main building was measured, and the amount 
of shielding to absorb out this radiation was determined. Second, the 
neutron intensities resulting from the small machine were measured at 
several points around the building, and the required additional shielding 
necessary to reduce this to a low level was also measured. Third, 
studies were undertaken to determine the type of construction of the 
new room which would give a low background for a reasonable cost. 


ROOM LOCATION AND ABSORPTION 


The radioactivity of the soil was measured by lowering an Nal 
detector into a 3-ft. deep hole in the ground. The presence of 2.62-Mev 
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and 1.46-Mev gamma rays from thorium and potassium, respectively, 
was quite apparent, and the background of this origin was somewhat 
greater than in the target toom. It was further determined that 2 in. 
of lead or 12 in. of normal concrete would reduce this background to a 
reasonable amount. 

Figure 1 shows the arrangement of the main building and the points 
A, B, and C indicate locations at which the neutron intensities resulting 
from the small generator were measured. It also shows the location 
of the new target room. The counting rate dropped from A to B to C, 
as would be expected because of the better shielding, and the rate at C 
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Fic. 1. Floor plan of the old and new target areas. 
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appeared to be low enough to enable one foot or so of concrete to reduce 
the effect of the remaining neutrons to an acceptable amount. In 
designing a new target room, and in order to maintain the same floor 
level as in the main building and a reasonable ceiling height of 9 ft., 
about 3 ft. of the room plus the ceiling thickness had to be above ground. 
It was realized that because of uncertainties in the measurements, 1-ft. 
concrete walls and ceiling might not be sufficient, but it was expected 
that enough dirt could be piled around and on top of the protruding 
structure to take care of this possibility. Background measurements 
upon completion proved this to be unnecessary. 


4 
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CONCRETE STUDIES 


In the appendix is given a complete report on the concrete studies. 
In brief, sources of low radioactivity content sand and gravel were 
located and concrete blocks were prepared from these ingredients. 
From these a cell of about 1 cubic meter volume was prepared and the 
radiation within this cell measured and compared with that in the old 
target room and that in a lead cell.?. This block cell, over the energy. 
range from 440 kev to 3 Mev, gave a reduction in background by a 
factor of 3 to 4 over the corresponding background in the old target 
room. This factor was not corrected for the natural radioacivity 
within the Nal crystal and photomultiplier and for leakage through the 
cracks of the cell. However, even a factor of 3 in reduction seemed 
desirable considering the other advantages of the new room. Since 
the cost of the special concrete as compared to the cost of normal con- 
crete was not excessive, plans were drawn up for the new room and 
estimates of cost obtained. These figures came within the allowable 
amount. Therefore, the project was approved. Figure 1 shows the 
general features of this room. 


PULSE HEIGHT DISTRIBUTION OF 
ROOM BACKGROUND 


@ Old target room 
New target room 


HOURS 


COUNTS FOR 


PULSE HEIGHT - VOLTS 
Fic. 2. Pulse-height distribution of the y-ray background in the old and new target 
rooms. The indicated energies are those for the photo peaks of K® (1.46 Mev) and thorium 
(2.62 Mev). 


2 It will be observed that in view of the law of inverse squares, and of the fact that the 
area of the walls of the cell are proportional to the square of the linear dimensions, observations 
of background taken in a cell will, to a first approximation, be independent of its size for an 
assigned wall thickness. 
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Upon completion, the natural background within the new room was 
measured using a 100 channel analyzer. The results are shown in 
Fig. 2 along with the spectrum obtained in the old target room. It 
may be seen that the background in the new room is lower than that 
in the old room by a factor of at least 5, and without considering the 
natural radioactivity within the crystal and photomultiplier. Pre- 
liminary tests also indicated that the effect of the neutrons produced by 
the small generator is very low. 


BEAM EXTENSION SYSTEM AND OPERATION 


While the new room was under construction, the system required 
to bring the charged particle beam into this room was being developed. 
Since the targets to be used in the new room would be about 35 ft. 
from the exit control slits of the analyzing magnet, provisions had to 
be made for positioning and focussing of the beam. Also, it was evident 
that means might have to be developed to take out the bounce of the 
beam between the control slits. This bouncing is the result of the 
relatively slow regulation caused by the large distance between the 
control corona needles and the central electrode of the machine. The 
beam is positioned by a system of horizontal and vertical deflector 
plates just following the control slits. A double system is used in the 
horizontal plane because of a slight misalignment of the axis of the old 
and new beam tubes. The beam focussing is taken care of by an elec- 
trostatic strong focussing lens. The original plan was to eliminate the 
beam bounce by use of a pair of slits in the new room and the vertical 
deflectors, that is, signals taken from these slits would apply correction 
voltages to the vertical deflectors. Actually it has been found that this 
circuitry is not required since the strong focussing lens seems to take 
care of the bounce. 

The beam as it appears in the new room is reasonably well focussed 
into a spot about § in. in diameter with 70 to 80 per cent transmission 
as compared to the current immediately following the generator control 
slits. With a 3-Mev 4ya beam of protons on a wolfram flap within the 
new room, neither the neutron nor gamma ray background is increased 
appreciably. 

CONCLUSION 

It appears that the new arrangement for the large generator will 
meet all of the requirements set up initially. In addition, there are 
other advantages which have become apparent. Actually there is 
considerably more space in this new room and less scattering material, 
that is, the beam analyzer magnet and the charged particle spectrometer 
are absent. Also, the long beam path provides for convenient methods 
of pulsing the beam. 
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APPENDIX 


The information here given may be of service to others who are faced with problems 
analogous to that here involved. To recapitulate, the elimination of background envisaged 
the construction of a new target room for the large generator, a room isolated by distance and 
shielding from the effects of the two machines located in the main building. This room had 
to be partially underground, and since the soil is quite radioactive, the walls and floor of the 
room had to be sufficiently thick to absorb out the soil radiation and yet be low in activity 
themselves. 


Low Background Concrete 
The following types of construction have been considered : 


1. Poured concrete 

2. Timber with backfill of sand 

3. Marble lining of (1) 

4. Silica brick (asphaltic concrete) 


The second method seem impractical and was not considered further. The third method 
seemed feasible if the first should not turn out too well; in any case, marble could be added later. 


Cost of Silica Brick 
Because silica brick might well be somewhat lower in radioactivity than a special concrete, 
the following cost figures were obtained from General Refractories: 


Size in inches $ per cu. ft. 
xX 3 4.06 
4.30 
iz x%6x%3 3.50 
IZ 3.78 


The density of silica brick is about 103 Ibs/cu. ft. Normal poured concrete would cost 
about $0.54 per cubic ft. and the density would be about 150 lbs. per cubic ft. Taking into 
account the densities, the cost of the silica bricks would be 8 to 10 times the cost of the poured 
concrete. This does not take into account the relative cost of installation which again would 
favor the concrete. Furthermore, it was expected that the cost of the special concrete we 
would require would not be more than 14 times the cost of normal concrete. Consequently, 
the use of silica bricks was not considered further. The remainder of the report concerns the 
studies on poured concrete. 


Materials of Low Radioactivity for Concrete 

Potassium and the thorium and uranium series provide all or nearly all of the natural 
radioactivity to be expected in concrete. Gray Portland cement normally contains }—} per 
cent potassium oxide and can only by good luck be chosen from a low potassium batch without 
extra cost. White cement, according to the one sample treated, runs somewhat lower in 
radioactivity but at about twice the cost of the gray cement. Thus the cement, which makes 
up about } of the concrete by weight, sets a rather firm lower limit on the activity of low cost 


concrete. 

Sand and gravel or crushed stone offer two possibilities for selection because silica and 
calcium carbonate do occur fairly pure in commercially available quantities. Consequently, 
sources of pure silica sand and gravel and of crushed limestone have been sought out within a 
reasonable distance. The local gravel and stone is, for the most part, river gravel or ‘‘granite.” 
The former contains pebbles of all sorts, many containing feldspar which in turn contains 
potassium. ‘Granite’ contains much feldspar and often traces of uranium and thorium 


minerals. 
Sand which is nearly pure silica is comparatively common since silica is resistant to 
leaching and attrition. Gravel, even if white and essentially silica, will contain more inclusions 


of radioactive grains. 
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Samples were obtained from local areas which supplied all of the local concrete mixing 
plants and from more distant sources towards the glass-producing area of New Jersey. 


Testing of Materials 
The various materials were tested by filling two tin pans with the samples and placing 
these pans in a lead enclosure in constant relationship to a 2 X 1j-in. Nal detector. The two 
pans held about 11 lb. of cement, 17 Ib. of sand, or 16 Ib. of gravel. The justification for using 
comparatively smal! spot samples of concrete ingredients are: 
1. The materials have been fairly well mixed, both by local geological processes, and 
screening and washing operations. 
2. It would be practically impossible to arrange for the use of exactly the material from 


which a more elaborate sample was taken. 
3. A few duplicating samples taken at different times during this work gave similar results. 
Of course it is always possible that the sample may include an exceptional lump of uranium 
or thorium material in a single stone. The use of duplicating samples however, eliminated 


such a chance occurrence. 

The pulses from the cathode follower of the photomultiplier were fed into a modified 
Atomic Instruments Model 218 linear amplifier. This amplifier has two discriminators. One 
discriminator was set for pulses greater than 440 kev and the other for pulses greater than 880 
kev. The high level output was fed into a single channel pulse height analyzer which was set 
for pulses greater than 1540 kev. With these pulses fed into three scalers, the number of 
pulses in the intervals of 440 to 880 kev and 880 to 1540 kev were determined. The gain was 
reduced by a factor of two to give the pulses in the intervals from 880 to 1760 kev and 1760 to 
3080 kev. 

All crucial samples were run sufficiently often in juxtaposition with blank runs to give 
assurance of the meaningfulness of the results for present purposes. The errors in sampling 
and procuring of these materials would not justify more elaborate precautions or the expendi- 


ture of more time. 


The results for a number of samples are given in Table I. Several samples are included 


TABLE I.—Radioactivity in Various Materials. (Counts with Sample in 
Place Minus Counts without Sample.) 


Counts per 10-Min. Interval 


0.44 to 0.88 0.88 to 1.54 0.88 to 1.76 1.76 to 3.08 
Sample and No. Mev Mev Mev Mev 
16 Allentown cement 133 63 62 0 
4 White cement 74 23 25 8 
1 Wenonah sand 108 98 97 —3 
3 Pine Valley sand 333 29 34 8 
12 Millville bar sand 13 2 24 —3 
13 Millville coarse sand 55 13 13 —3 
2A Tallytown 1-in. gravel 63 87 77 19 
3 Pine Valley gravel 30 45 21 —2 
6 Valley Forge limestone 102 10 31 —3 
7 Cedar Hollow limestone 113 11 - 
9 Wenonah }-in. gravel 309 61 73 8 
11 Millville }-in. gravel 63 —8 —9 1 
14 Millville 2-in. gravel 21 —8 -1 8 
18 Chickies quartzite 248 118 146 14 
21 Glen Mills 1-in. stone 277 152 176 8 
15 Ironstone 475 181 206 72 
17 Pegmatite 1067 1029 1044 46 
20 1.4 Ib. potassium salts 742 1130 1062 16 
Soil south of generator 246 109 102 10 
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at the end for general information. It appears from these last results that the precipitated iron 
from the leaching solutions acting on the coastal plain sediments may concentrate considerable 
radioactivity, and that the lightest possible colors of sand and gravel are therefore desirabie. 
The desirability of eliminating pegmatite (light colored veins) from crushed stone is also 
indicated. 

The best sand and gravel appear to be those from the Millville Sand and Gravel Company 
at Port Elizabeth, New Jersey, whereas the best cement is the white variety. For the concrete 
test blocks studied in the next section, the Millville bar sand and j-in. gravel and the gray 
cement were used, the latter for cost reasons as compared to the ultimate effect. 


Background in Concrete Enclosure 
Twenty concrete blocks 8 X 8 X 12 in. in size were made using a 1:2:3 mix by volume. 
Six of these test blocks were laid on a platform of 3 in. of lead. Eight more blocks on end were 
placed on this square to form an 8 X 8 X 12-in. enclosure. Six more blocks, arranged like the 
bottom, constituted the cover. 
The results of placing the scintillation counter in this enclosure as compared to a lead 
enclosure and the open room are given in Table Il. The open room correspends to the activity 


TABLE II.—Background Measurements.* 


Counts per 10-Min. Interval 


0.44-0.88 Mov 0.88--1.54 Mev 1.76-3.08 Mev 
Counter on floor of main building 9855 5346 994 
Counter in lead shield 509 262 69 
Counter in 8 X 8 X 12-in. concrete 3146 1199 267 


enclosure 


* No correction has been made for activity within the crystal. 


of the concrete of the main building plus the amount of radiation from the soil coming through 
the floor and walls. The results were well reproducible on successive 10-min. runs, the 
errors on the larger figure being about +2.5 per cent and those on the lower being about 10 
per cent. Two effects must be considered in determining the real effect of the concrete en- 
closure. First, cracks exist between the blocks and around the leads to the photomultiplier; 
in the case of the lead enclosure there was better overlapping of the lead blocks. Second, 
undoubtedly much of the counting rate within the lead enclosure is from the Nal crystal itself. 
This effect, whatever it may be, should be subtracted from all of the results, and the result of 
such correction is to lower the estimated background to an estimate which, percentage wise, is 
greater for small background than for large ones. 


Material Requirements and Cost 

Before making a final decision on the construction of this new room, some idea of the cost 
of the special concrete must be made. It had already been decided that the walls, floor and 
ceiling must be one foot thick and that the inside dimensions should be 20 ft. square by 9 ft. 
high. The amount of concrete required was then: 


Vol., cu. ft. 
Floor and roof 22 X 22 ft. 968 
Four walls 9 X 21 ft. 756 


1724 cu. ft. 


With a density of 149 lb./cu. ft. this gives a total weight of about 129 tons of concrete. 
The mix decided upon—one part cement, two parts sand, and three parts gravel—has a com- 
pressive strength rating of 2800 psi. and costs between $14 and $15 per cubic yard. This would 
give as the concrete cost for normal material of about $1000. 

The actual amounts of the individual components are about 17 tons of cement, 43 tons of 
sand, and 76 tons of gravel. The following costs were quoted by the Millville Sand and 
Gravel Company at the Port Elizabeth, New Jersey pit. 
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$/ton 
0.60 
1.80 


Pure bar sand 
?-in. gravel 


Freight to Philadelphia 3.20 


Our presumed supplier of this concrete gets his sand from Wenonah, New Jersey, and his 
gravel from Glen Mills, Delaware County, Pennsylvania. Therefore, in estimating the greater 
cost of materials we may say that the freight or trucking cost on the sand would be about 
doubled while that for the gravel would be about five times. Therefore, the extra charge per 
cubic yard of concrete would be about $4.00 without allowing for special runs of the mixing, 
and storage of the special sand and stone. This latter is important since it may take as long 
as six weeks to make all the pours. It is not expected however that the cost will more than 


double, that is, not more than $30 per cubic yard. 
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Bulletin 192, Jan.6-10. 1958. 


Com puters 


228. “A Rapid Statistical Data Processing 


System for Radio Propagation Re- 


search,”’ by S. Charp, in “Proc. Sym- 
posium on Statistical Methods in Radio 
Wave Propagation,” U. of Calif., 1959. 


Engineering Psychology 


197. “A Review and Summary of Tracking 
Research Applied to the Description of 
Human Dynamic Response,” by D. T. 
McRuer and E. S. Krendel, in “IRE 
WESCON Convention Record, Part 4, 
1958," p. 254. 
“A Survey of Suggested Mathematical 
Methods for the Study of Human Pilot's 
Responses,”’ by E. S. Krendel, in ‘‘Re- 
port of 2nd Piloted Aircraft Flight Con- 
trol System Symposium, June 1952,” 
BuAer Report AE-61-5-III. 
“Human Dynamics and System Syn- 
thesis,” by E. S. Krendel and D. T. 
McRuer, in ‘Proc. of National Syms 
posium on Human Factors in System- 
Engineering, HFSA, IRE, Phila., Pa. 
Dec., 1957." 
“Human Transfer Functions,” by E. S. 
Krendel, in “Proceedings, 2nd Annual 
Human Factors Discussion Group, 11th 
Annual International Air Safety Semi- 
nar, Nov. 1958,” Flight Safety Founda- 
tion, New York, p. 17. 
“The Human Operator as a Servo Sys- 
tem Element,”’ by D. T. McRuer and 
E. S. Krendel, Jour. FRANKLIN INsT., 
267, 381 and 511 (1959). 
Friction-Lubrication 
11. “What About Friction?: I, Classical 
Laws; II, Friction as the Result of 
Molecular Forces; III, Friction Saw, 
Electrical Theory, Lubrication,"’ by F. 
Palmer, Am. J. Phys., 17, 181, 327, 
336 (1949). 
“Mixed Friction Conditions in Lubrica- 
tion,” by D. D. Fuller, Lubrication 
Eng., Sept-Oct. 1954. 
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155. “AStudy of the Effect of Wear Particles 

and Adhesive Wear at High Contact 

Pressures,”’ by E. B. Sciulli and G. M. 

Robinson, Trans. ASLE, 1, 312 (1958). 

“Determination of Optimum Propor- 

tions for Hydrostatic Bearings,” by A. 

M. Loeb and H. C. Rippel, Trans. 

ASLE, 1, 241 (1958). 

211. “The Influence of the Molecular Mean 
Free Path on the Performance of Hy- 
drodynamic Gas Lubricated Bearings,” 
by A. Burgdorfer, Trans. ASME, J. 
Basic Eng., 1, 94 (1959). 

215. “Axial, Relative Motion of a Circular 
Step Bearing,” by L. Licht, Tvrans. 
ASME, J. Basic Eng., 1, 109 (1959). 


166. 


General Physics 

203. ‘Technical Data on Border-Punched 
Cards,"’ by G. Cohn, Am. Documenta- 
tion, X, 116 (1959). 

“The Theory of Low-Temperature Bo- 
lometer Detectors Applied to the Meas- 
urement of Low-Level RF Power,” by 
D. L. Birx and N. Fuschillo, Trans. 
IRE Instrumentation, 1-7, 310 (1958). 
213. “Early Detection of Fatigue in Metal 
Alloys by Ultrasonics,” by K. Sittel, 
Nondestructive Testing, XVII, 165 (1959). 
“On the Reaction of Fluids and Fluid 
Jets,” by R. Ekesrgian, Jour. FRANK- 
LIN INstT., 237, 385 (1944). 
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Magnetism 


139. ‘Magnetic Properties of Some Iron and 

Iron Oxide Micropowders,”’ by R. B. 

Campbell, H. Amar, A. E. Berkowitz 

and P. J. Flanders, AIEE Spec. Pub. 

T-91, 1956. 

222. ‘Magnetic Properties of Some Ferrite 
Micropowders,” by A. E. Berkowitz 
and W. J. Schuele, J. Appl. Phys., 30, 
134S (1959). 


Mechanical Engineering 
230. “Vibration of Beam 
with Time-Dependent Boundary Con- 
ditions,” by T. C. Yen and S. Kao,* 
ASME Paper No. 59-APM-22. 


Mass Systems 


Metallurgy 
202. “Cellular Substructure in Zn Crystals 
Grown from the Melt,” by V. Damiano 
and M. Herman, Trans. AIME, 215, 
136 (1959). 
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212. “Etch Patternsin Zinc,” by V. Damiano 
and M. Herman, Jour. FRANKLIN INST., 
267, 303 (1959). 
“Precipitation in a Beta-Brass-Fe Al- 
loy,”” by A. E. Berkowitz and P. J. 
Flanders, J. Appl. Phys., 30, 1115S 
(1959). 


Metals Physics 


110. ‘“Self-Diffusion in Single and Polycrys- 
tals of Zinc at Low Temperatures,” by 
F. E. Jaumot, Jr. and R. L. Smith, 
J. Metals, 8, 164 (1956). 

“Preparation of Fine Particles from Bi- 
metal Oxalates,”’ by W. J. Schuele, J. 
Phys. Chem.,63, 83 (1959). 

“Properties of Multidomain Particles,” 
by H. Amar, /. Appl. Phys., 30, 139S 
(1959). 
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Plastic Deformation of Metalst 


201. “A Study of Dislocations in Thin Alu- 
minum Foils Elongated in the Electron 
Microscope,” by H. G. F. Wilsdorf, in 
“Symposium on Advances in Electron 
Metallography,”’ Spec. Tech. Publ. No. 
245, ASTM, 1958, p. 43. 
“A Study of Slip Lines in Alpha-Brass 
as Revealed by the Electron Micro- 
scope,”’ by J. T. Fourie* and H. G. F. 
Wilsdorf, Acta Met., 7, 339 (1959). 
218. “Dislocation Patterns in Deformed 
Alpha-Brass,”’ by H. G. F. Wilsdorf, to 
be published in “Proc. Symposium 
on Internal Stresses and Fatigue in 
Metals,”” Elsevier Publ. Co., 1959, in 
press. 
“Experimental Techniques for the Plas- 
tic Deformation of Metal Foils in the 
Electron Microscope,”” by H. G. F. 
Wilsdorf, L. Cinquina, and C. J. Varker 
to be published in ‘Proc. International 
Congress of Electron Microscopy, 1958.”’ 
“The Effect of Plastic Deformation on 
the Resistivity of Copper-Palladium Al- 
loys,” by E. Klokholm and B. Hyatt, 
to be published in Trans. AIME. 


Rheology 


205. ‘Use of Rheological and Other Data in 
Asphalt Engineering Problems,” by J. 
G. Brodnyan, Highway Research Board 
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217. “The Rheology of Various Solutions of 
Cellulose Derivatives," by J. G. Brod- 
nyan, F. H. Gaskins, W. Philippoff and 
E. G. Lendrat,* to be published in 
Trans. Soc. Rheology, 11, in press 

“The Capillary Experiment in Rhe- 
ology,” by W. Philippoff and F. H. 
Gaskins, to be published in Trans. Soc. 
Rheology, 11, in press. 

“The Experimental Check of Theories 
of the Viscosities of Solutions,”” by W. 
Philippoff and F. H. Gaskins, J. Phys. 
Chem., 63, 985 (1959). 
“Instrumentation for Rheological In- 
vestigation of Viscoelastic Materials,” 
F. H. Gaskins and W. Philippoff, Jnd. 
Eng. Chen:., 51, 871 (1959). 
“Stress-Optical Analysis ef Fluids,”’ by 
W. Philippoff, Ind. Eng. Chem., 51, 883 
(1959). 
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Semiconductors 


196. “Transport Properties of the Pseudo- 
Binary Alloy System Bi,Te;_ySe,,"" by 
N. Fuschillo, J. N. Bierly and F. J. 
Donahoe, J. Phys. Chem. Solids, 8, 430 
(1959). 

“Transport Phenomena in Semiconduc- 
tors for Warped Energy Surfaces,"’ by 
G. E. Tauber and L. Soffer,* J. Phys. 
Chem. Solids, 8, 138 (1959). 

“A Generalized Variational Principle 
for Transport Phenomena,” by G. E. 
Tauber, Canadian Journal of Physics, 
36, 1308 (1958). 


Thermoelectricity 


231. “Theory of the A. F. Joffe Method for 
Rapid Measurement of the Thermal 
Conductivity of Solids,” by W. F. G. 
Swann, Jour. FRANKLIN INstT., 267, 
363 (1959). 
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“Convertible Aircraft,”” by R. S. Barnaby, 
S. M. Berkowitz and W. H. Colcord, Aero. 
Eng. Rev., 8, 4 (1949). 
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Micelles,” by W. Philippoff, Faraday Soc. 
Disc., No. 11, 96 (1951). 

“Reactions Involving Ozone Nitrogen Di- 
oxide, and Organic Compounds at Low Con- 
centrations in Air,”” by P. L. Hanst, E. R. 
Stephens and W. E. Scott, J. Air Poll. Cont. 
Assn., 5, 219 (1956). 

“Recent Developments in the Study of the 
Organic Chemistry of the Atmosphere,” by 
E. R. Stephens, W. E. Scott, P. L. Hanst 
and R. C. Doerr, API Reprint, 1956. 

“Journal-Bearing Operation at Superlaminar 
Speeds,’’ by M. I. Smith* and D. D. Fuller, 
Trans. ASME, 78, 469 (1956). 

Theory and Practice of Lubrication for Engi- 
neers, by Dudley D. Fuller, John Wiley & 
Sons, Inc., 1956. 

“Film-Pressure Distribution in Grease-Lubri- 
cated Journal Bearings,"’ by G. Cohn and 
J. W. Oren,* Trans. ASME, 71, 555 (1949). 

“The Design of Lightweight Trains,’ by R. 
Eksergian, Trans. ASME, 56, 667 (1934). 

“The Dynamical Analysis of Machines,"’ by 
R. Eksergian, International Congress of 
Applied Mechanics, Stockholm, 1930. 


“The Fluid Torque Converter and Coupling,” 
by R. Eksergian, Jour. FRANKLIN INsT., 
235, 441 (1943). 

“General Dynamical Considerations Applied 
to Piezo-Electric Oscillations of a Crystal in 
an Electrical Circuit,”” by W. F. G. Swann, 
Jour. FRANKLIN INsT., 242, 167 (1946). 

“Certain Matters Concerning Scaling in the 
Magnetron with Special Reference to the 
Relative Efficiency of Magnetrons of Differ- 
ent Sizes,” by W. F. G. Swann, Jour. 
FRANKLIN INsT., 246, 149 (1948). 

“Fluctuations in Proportional Counters,"’ by 
W. F. G. Swann, Jour. FRANKLIN INsT., 
249, 133 (1950). 

“Electrets,” by W. F. G. Swann, Phys. Rev., 
78, 811 (1950). 

“On Certain Matters Pertaining to Electrets,”’ 
by W. F. G. Swann, Jour. FRANKLIN INsT., 
250, 219 (1950). 

“Engineering and Pure Science,"’ by W. F. G. 
Swann, Physics Today, 4, 9 (1951). 

“The Engineer and the Scientist," by W. F. G. 
Swann, Jour. Engineering Education, 43, 
241 (1952). 

“Fundamentals in the Behavior of Electrets,"’ 
by W. F. G. Swann, Jour. FRANKLIN INsrT., 
255, 513 (1953). 

“The Case of Charged Density Distribution 
Versus Semi-Permanent Polarization as a 
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Basis for Electret Behavior,” by,W. F. G. 
Swann, Jour. FRANKLIN INsT., 256, 167 
(1953). 

“Engineering in our Epoch,” by W. F. G. 
Swann, University of Toronto Press, 7th 
Wallberg Lecture, 1954. 

“Stress and Strain,” by W. F. G. Swann, J. 
Appl. Phys., 25, 1108 (1954). 

“Creative Intuition in Engineering Science,” 
by W. F. G. Swann, Am. Management 
Assn., Report No. 3, Planning Ahead for 
Profit, 1957, p. 156. 

“The Formation of Mechanical Twins,"’ by 
N. Thompson* and M. Hingley, Acta Met., 
3, 289 (1955). 

“On the Representation and Extrapolation of 
Creep Data of Metals and Technical Al- 
loys,"’ by P. Feltham* and J. D. Meakin, 
Rheologica Acta, 2, 176 (1958). 

“Creep in Copper,”’ by P. Feltham* and J. D. 
Meakin, Acta Met., in press. 

“Quenching Stresses and Strains in Ionic 
Crystals,"’ by B. H. Kear and P. L. Pratt,* 
Acta Met., 6, 457 (1958). 

“Quench-hardening in Sodium Chloride Crys- 
tals," by B. H. Kear and P. L. Pratt,* 
Phil. Mag., 4, 56 (1959). 

“Some Dislocation Interactions in Simple 
lonic Crystals," by B. H. Kear, A. Taylor* 
and P. L. Pratt,* Phil. Mag., 4, 665 (1959). 

“Work-hardening in Sodium Chloride Crys- 
tals," by B. H. Kear and P. L. Pratt,* 
in press. 

“The Half Life of I'* Labeled Normal Human 
Gamma Globulin in Patients with Hepatic 
Cirrhosis," by J. N. Bierly, P. Havens,* 
W. Dickensheets* and T. Eberhard,* J. Jm- 
munology, 73, 256 (1950). 
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“The Mechanism of Exchange of Hydrogen 
between Nitrogen and Oxygen Atoms,” by 
C. G. Swain,* J. T. McKnight,* M. M. 
Labes, and V. P. Kreiter,* J. Am. Chem. 
Soc., 76, 4243 (1954). 

“The Mechanism of Exchange of Hydrogen 
between Ammonium and Hydroxyl Groups. 
1," by C. G. Swain* and M. M. Labes, 
J. Am. Chem. Soc., 79, 1084 (1957). 

“Molecular Compounds. V. Solvent Effects on 
Complexing of 1,3,5-Trinitrobenzene with 
Aniline and with Naphthalene,” by S. D. 
Ross* and M. M. Labes, J. Am .Chem. Soc., 
77, 4916 (1955). 

“Molecular Compounds. VI. A Reinvestiga- 
tion of the Picryl Chloride-Hexamethy]- 
benzene Complex; The Effect of Triethyla- 
mine on the Equilibrium Constant for Com- 
plex Formation,” by S. D. Ross,* M. M. 
Labes and M. Schwarz,* J. Am. Chem. 
Soc., 78, 343 (1956). 

“Molecular Compounds, VII. Interactions 
between: 1,3,5-Trinitrobenzene and Pyri- 
dine-N-Oxides in Chloroform Solution," by 
S. D. Ross,* D. J. Kelley* and M. M. Labes, 
J. Am. Chem. Soc., 78, 3625 (1956). 

“Molecular Compounds. VIII. Evidence for 
the Existence of 2:1 and 1:2-Complexes in 
Chloroform Solutions of 1,3,5-Trinitroben- 
zene and N,N-Dimethylaniline,” by S. D. 
Ross* and M. M. Labes, J. Am. Chem. 
Soc., 79, 76 (1957). 

“The Solvolyses of Some Tertiary Halides in 
Dimethylformamide and N-Methylpropi- 
onamide,” by S. D. Ross* and M. M. Labes, 
J. Am. Chem. Soc., 79, 4155 (1957). 

“The Dimethylamine-1,3,5-Trinitrobenzene 
Complex in Dioxane,”” by M. M. Labes and 
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FRAMES AND ARCHES, by Valerian Leontovich, 
with foreword by S. P. Timoshenko. 472 
pages, diagrams, 6 X 9 in. New York, 
McGraw-Hill Book Co., Inc., 1959. Price, 
$20.00. 


Even though the exact equations for the 
elastic analysis of arches was developed (by 
Weyrauch) as early as 1879, their complexity 
precluded their extensive applications to ac- 
tual design problems, which were often solved 
by more or less elaborate graphical methods. 
This state of affairs was summarized in the 
transactions of the American Society of Civil 
Engineers in 1898: 

“This modern and most exact method, how- 
ever, is not free from criticism. While the 
fundamental principles of the theory are al- 
most axiomatic, their final application to the 
solution of stresses is extremely complicated, 
so much so that few engineers can be credited 
with the patience and earnest endurance to 
master either the method or the solution of a 
problem to which it is applied. . . .” 

Even today, the level of ‘“‘patience and en- 
durance”’ required for carrying out the analy- 
sis of indeterminate structures by such meth- 
ods as slope deflection, virtual work and mo- 
ment distribution is still considerable. Every 
practicing design engineer will therefore wel- 
come the work of the present author, who can 
be credited with a major step towards the 
reduction of the work of computation, as well 
as the chances of error, in the analysis of 
frames and arches. 

Condensed, yet clear solutions for shears 
and moments are given for twenty principal 
types of indeterminate structures. Over 400 
loading conditions are dealt with, including 
concentrated and partial distributed loads, 
sidesway loads, and loads due to displacement 
of supports and changes in temperature. The 
simple algebraic operations developed by the 
author lead to determination of the redun- 
dants, from which the moments and forces for 


any section of the structure may be quickly 
calculated. 


The book is divided into two main parts: 
the first comprising 12 sections, treats frames 
and arches with members of constant cross 
section, with both hinged and fixed supports, 
including portal frames, trapezoidal frames, 
gable frames, parabolic arches, and parabolic 
frames. 

The second part presents what may be con- 
sidered the first convenient and_ practical 
method of analysis for frames and arches of 
variable cross section. Sections 13 to 19 deal 
with frames with straight members. An ex- 
tension of the method of elastic centers (called 
method of elastic parameters), developed by 
the author, not only allows treatment of mem- 
bers with one or two haunches, straight or 
parabolic, but enables the designer to analyze 
members not falling into these classifications 
as well. Portal, trapezoidal and gable frames 
with hinged and fixed supports are dealt with. 
Eight tables giving elastic parameters and 
load constants for the standard members are 
given in the Appendix. 

Sections 20 to 24 treat frames with curved 
members, including parabolic arches and 
frames, with hinged and fixed supports, which 
may be solved with the help of 12 additional 
tables in the Appendix. These allow reduc- 
tion to a standard form of any symmetrical 
arch, on the basis of the ratios of rise to span, 
and thickness at crown and springing. Fac- 
tors for a great variety of loading conditions 
are tabulated, including concentrated, dis- 
tributed and lateral loads, and factors for the 
effect of axial forces are also given for the three 
main types of arches considered by the author : 
arches of constant cross section; ‘prime’ 
arches, in which the cross section varies at a 
more or less constant rate from the crown to 
abutments; and “quadratic”’ arches in which 
the cross section increases rapidly near the 
abutments. 

These details can only give a superficial idea 
of the comprehensiveness of the work. Ample 
cross references are given, and a variety of 
examples are worked out in full to show the 
application of the author’s method. The fig- 
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ures and notations are clear and unambiguous, 
and the attractive type does justice to the sub- 
ject matter. An appealing feature is the in- 
clusion, before each chapter, of photographs of 
large and unusual structures, either completed, 
or in the process of construction, or, in some 
cases, both. The concise captions draw atten- 
tion not only to their salient engineering de- 
tails, but also to their aesthetic aspects, thus 
showing again that something of an artist 
is hidden in every outstanding structural 
engineer. 
K. L. CAapPEL 
The Franklin Institute Laboratories 


Cerenkov Rapration, by J. V. Jelley. 304 
pages, diagrams, plates, 53 X 8} in. New 
York, Pergamon Press, 1958. Price, $10.00. 


The appearance of this excellent book, the 
only monograph devoted exclusively to the 
subject of Cerenkov radiation, is indeed timely, 
since it coincides with the awarding of the 
1958 Nobel Prize in Physics for the discovery 
of this phenomenon 24 years earlier. 

The effect first manifested itself as a feeble 
bluish-white light emitted by transparent sub- 
stances under the influence of the radiations 
from radioactive substances. When an elec- 
trically-charged particle, such as an electron, 
passes through a medium with a velocity ex- 
ceeding the velocity of light in that medium, 
a kind of optical shock-wave is produced. 
This is quite analogous to the shock-wave 
created in air by a missile travelling at super- 
sonic velocity, or the bow-wave from a ship 
moving faster than the velocity of the surface 
waves. The light output can be measured by 
a photomultiplier tube, and information about 
the nature and properties of the individual 
particle producing the pulse can thus be 
deduced. 

The Cerenkov detector has recently become 
an important tool in the hands of the cosmic 
In 1951, Dr. Jelley 
succeeded in detecting single fast charged par- 
ticles for the first time by this method, and 


ray and nuclear physicist. 


since then he has made numerous significant 
contributions to the development of the Ce- 
renkov counter. An eminent authority on the 
subject, he is exceedingly well qualified to dis- 
cuss Cerenkov radiation and its applications. 

This comprehensive book will be a welcome 
addition to the reference library of the experi- 
mentalist who is interested in exploiting the 
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possibilities of this new device. It contains a 
compendium of useful technical information, 
including some, such as tables summarizing 
characteristics of probably all commercially- 
produced photomultiplier tubes, which will re- 
quire frequent revision. Relevant curves and 
data from standard reference sources with 
which workers in the field are thoroughly fa- 
miliar are collected in the Appendix, and there 
is an extensive and apparently complete bibli- 
ography, including the considerable work 
carried out in the USSR. 

From the point of view of fundamental 
physics, the treatment is essentially historical. 
Chapter 1 is a general introduction, followed 
by a presentation of theoretical considerations 
in Chapters 2 and 3. After a survey of later 
experimental work in Chapter 4, there is an 
abrupt switch to a discussion of the technology 
of Cerenkov counters in the next four chapters. 
These deal with photomultipliers, radiators, 
design factors involved in the most efficient 
transfer of light from the radiator to the multi- 
plier, and descriptions of some _ practical 
counters. Although the logical justification 
for next inserting a chapter on Cerenkov radi- 
ation in the atmosphere is not apparent, this 
review article is of interest in its own right to 
at least one segment of the readers of this 
book, cosmic ray investigators. The text then 
again reverts to counters, discussing the little- 
used potentialities of the gas Cerenkov counter 
which should have some application in experi- 
ments conducted with the very large acceler- 
ating machines at present under construction. 
The final chapter is a potpourri of miscellane- 
ous topics, some controversial, others highly 
speculative. Needless to state the author is 
competent to distinguish fact from fancy, and 
indeed does so. 

For whatever it is worth, the reviewer would 
like to state that he has already derived much 
benefit from the use of Dr. Jelley’s book in 
his own work. 

Martin A, POMERANTZ 
Bartol Research Foundation 


SouNDING Rockets, by Homer E. Newell, Jr. 
334 pages, illustrations, 6 X 9 in. New 
York, McGraw-Hill Book Co., Inc., 1959. 
Price, $12.50. 


Man lives at the bottom of a turbulent 
ocean of air. Everything coming through this 
atmosphere is attenuated, scattered or altered. 
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Beyond the atmosphere we find the full range 
_of the electromagnetic spectrum. Thus, on 
the surface of the Earth most of these radia- 
tions have been filtered out. We are told that 
if the entire electromagnetic spectrum was 
considered a keyboard with 60 keys, one single 
key represents that which we can see with our 
eyes and but a few more represent what can 
be detected with instrumentation from the 
Earth's surface. To obtaina definitive picture 
of what energies are present in space and their 
flux, it is necessary to ascend beyond the 
atmosphere. 

Until a few years ago this was impossible. 
The balloon was the only vehicle which could 
ascend to great heights to collect significant 
information. With the advent of the rocket 
this has been changed. 

Today scientists all over the world have 
been probing radiations beyond the atmos- 
phere. With the fantastic improvements in 
instruments we are able to parlay the short 
periods of time during which the sounding 
rocket spent above most of the atmosphere to 
draw a picture of what takes place in this 
region. 

One of the true pioneers in this work is the 
senior editor of this new book. Dr. Newell 
began his work in this field when the U. S. 
Army looked for something to put into the 
payload sections of the V-2's which had been 
brought over from Germany. His work as 
head of the Rocket Sonde Branch of NRL's 
rocket upper-air program gave him the back- 
ground to contribute to this science. 

Sounding Rockets is a basic book covering 
the sounding rockets used in all countries ex- 
cept, and for obvious reasons, the U.S.S.R. 
Together with his colleagues Dr. Newell has 
covered rocket theory, instrumentation tech- 
niques, handling and launching techniques 

One of 
the most significant tables this reviewer has 


and the future of sounding rockets. 
seen is Table 2-2. Here are listed the results 
obtained from sounding rockets. 

Dr. Newell states the object of study and 
then describes the methods used to attack the 
problem. Finally, the results obtained are re- 
vealed. When you consider that this was done 
for the determination of atmospheric pressure, 
density composition and temperature, winds, 
solar radiation, ionosphere, magnetic field of 
the Earth, etc., it proves to be the best synopsis 
and guide to the scientific attack on our lack 
of knowledge of our atmosphere. 
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Dr. Newell and his contributors have put 
together a commendable work. It will prove 
to be a most useful reference. 

I. M. Levitt 
The Franklin Institute Planetarium 


ProGress 1N METAL Puysics, Vol. 7, edited 
by Bruce Chalmersand R. King. 408 pages, 
diagrams, 6 X 93 in. New York and Lon- 
don, Pergamon Press Ltd., 1958. Price, 
$16.00. 


, As in previous years many research workers 


in physical metallurgy or metal physics will 
find at least one article among the five con- 
tained in the 7th volume of the series which is 
closely related to their own work. The timing 
of the articles could not have been better. Al 
the papers describe work in areas where ¢ 
great deal of activity is being displayed and 
where further rapid advances are expected 
during the next decade. 

Hobstetter’s article “. . . Equilibrium, Dif- 
fusion and Imperfections in Semiconductors” 
contains a good deal of solid state physics 
while the following paper by Jaffee, ‘The 
Physical Metallurgy of Titanium Alloys,” is 
on the extreme metallurgical side. However, 
the articles are written in clear language so 
that the former will make enjoyable reading 
also to the metallurgist and the latter to the 
physicist. The remaining reviews deal with 
“The Thermodynamics and Kinetics of Mar- 
tensitic Transformations’ (Kaufman and 
Cohen), “The Stored Energy of Cold Work” 
(Titchener and Bever) and ‘‘The Properties of 
Metals at Low Temperatures” (Rosenberg). 
All three articles have in common that the 
subject has been reviewed very critically, and 
that much thought kas been used to show 
new approaches for further research. 

The new volume is, in the opinion of the 
reviewer, the best balanced book of the series 
and is fully recommended. There is an in- 
crease of 50 to 408 pages as compared with the 
previous volumes. This expansion seems to be 
responsible in part for the increase of price by 
30 per cent, and many readers may find the 
price of $16.00 rather high. It is felt that 
the series may serve its purpose better by 
limiting the size of the single volume to 300 


pages and thereby keeping the book in reach 
of the graduate students who are a substantial 


group of its readers. 
H. G. F. Witsporr 
The Franklin Institute Laboratories 
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On NUMERICAL APPROXIMATION, edited by 
Rudolph E. Langer. 462 pages, 6} * 9} in. 
Madison, University of Wisconsin Press, 


1959. Price, $4.50. 


On Numerical Approximation consists of 
twenty papers delivered at the Symposium on 
Numerical Approximation conducted by The 
Mathematics Research Center, United States 
Army, at the University of Wisconsin, April 21 
to 23, 1958. The general theme was threefold: 
Linear Approximation, Extremal Approxima- 
tion, and Algorithms. An additional paper by 
R. Creighton Buck, Linear Spaces and Ap- 
proximation Theory, was included for the pur- 
pose of clarifying the role of functional analy- 
sis in its relation to present day investigations 
of approximation theory. 

Those contributing to the symposium and 
consequently to this book can be identified as 
some of the most outstanding figures in con- 
temporary approximation theory. It is well, 
then, to list their names, the titles of their 
papers and to add remarks where they may 
be appropriate. 

Alexander Ostrowski, On Trends and Prob- 
lems in Numerical Approximation 

In this brief paper the stage is set for the 
remaining papers of the symposium. In his 
historical account Ostrowski discusses the tre- 
mendous influence of the digital computer on 
the development of approximation art and 
science. He also admonishes those who, in 
their quest for results “‘with a capital R,"’ may 
be prone to abandon rigor. 

Zdenek Kopal, Operational Methods in Nu- 
merical Analysis Based on Rational Ap- 
proximations 

The aim of this paper is to develop new 
methods for numerical integration and differ- 
entiation of continuous functions in terms of 
their tabular differences. Remarkable ac- 
curacy is obtained by rational approximation 
using low order differences. 

Philip J. Davis, On the Numerical Integration 
of Periodic Analytic Functions 

The notion that the trapezoidal rule is the 
“best’”’ rule when integrating a periodic func- 

tion over its period is investigated in this 

paper. A comparison is made between the 
trapizoidal and the Gaussian families to de- 
termine which is better, according to a certain 

It is shown that the Gaussian is better 


norm. 
if the analytic periodic function possesses iso- 
lated singularities of a certain type, sufficiently 
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close to the interval of integration. The trap- 
ezoidal method is better for regular analytic 
periodic functions. 
Herbert E. Salzer, Some New Divided Differ- 
ence Algorithms for Two Variables 
Preston C. Hammer, Numerical Evaluation of 
Multiple Integrals 
An accounting is made of the somewhat 
deficient state of the art in this area of nu- 
merical approximation. 
Michael Golomb and Hans F. Weinberger, 
Optimal Approximation and Error Bounds 
This excellent paper displays the power and 
the “‘all-embracing’”’ characteristic of func- 
tional analysis, whose general purpose is : ‘Show 
to derive from a given amount of information 
about the unknown function wu the total range 
of possible values of the functional F(u).” 
Arthur Sard, The Rationale of Approximation 
More on functional analysis pertaining to 
approximations of the sort 7(f + 5f) = Ff, 
where f, 5f are functions, F is a given operator 
and T is an approximating operator. Several 
sets of hypotheses are formulated relative to 
T, F, f and 6f and the results are described. 
Joseph L. Walsh, On Extremal Approximations 
Eduard L. Stiefel, Numerical Methods of 
Tchebycheff Approximations 
Leslie Fox, Minimax Methods in Table Con- 
struction 
Theodore S. Motzkin, Existence of Essentially 
Nonlinear Families Suitable for Oscilla- 
tory Approximations 
I. J. Shoenberg, Jn Variation Diminishing 
Approximation Methods 
The property ‘variation diminishing" is 
defined and several approximation methods 
are shown to possess it, including the Bern- 
stein polynomials and the de la Vallee Poussin 
means. 
Michael Golomb, Approximation by Function 
of Fewer Variables 
J. C. P. Miller, Extremal A pproximations- 
A Summary 
This is a discussion of polynomial approxi- 
mation from the practical side. Smoothing 
and Numerical Differentiation are discussed ; 
also Quadrature, Tchebycheff approximation 
and rational approximations. 
R. Creighton Buck, Survey of Recent Russian 
Literature on Approximation 
This paper does not read as interestingly as 
the title might suggest to the casual reader of 
mathematical journals, if there are any. In- 
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deed, it reminds this reviewer of a sophisti- 

cated edition of an article in Variety usually 

called “What They're Doing and Where 

They're Playing."’ Nevertheless, it reflects a 

great amount of research on the part of the 

author and an excellent bibliography of writers 

in the U.S.S.R. (and other countries as well) 

in approximation theory has resulted. 

Friedrich L. Bauer, The Quotient-Difference 
and Epsilon Algorithms 

J. Barkley Rosser, Some Sufficient Conditions 
for the Existence of an Asymptotic Formula 
or an Asymptotic Expansion 

John W. Tukey, The Estimation of (Power) 
Spectra and Related Quantities 

Lothar Collatz, Approximations in Partial 
Differential Equations 

John Todd, Special Polynomials in Numerical 
Analysis 
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Because most first (and often only) courses 
in numerical analysis tend toward the practi- 
cal and away from the theoretical, the author 
in this paper proposes an outline of a course of 
a more theoretical nature on approximation. 
The need for this is indicated by the existence 
of many books such as the one under review. 

As a general observation on On Numerical 
Approximation, it must be granted that only 
a monumental presentation can result from 
such a notable lineup of contributors. It is 
written by experts and will be fully appreci- 
ated only by those who are equally expert or 
nearly equally expert. Not for browsing, it is 
unquestionably a warehouse of modern think- 
ing in approximation theory. 


LaurENCE C. McGINN 
The Franklin Institute Laboratories 


PUBLICATIONS RECEIVED 


SERVOMECHANISMS AND REGULATING SysTEM DesiGN, by Harold Chestnut and Robert W. 
Mayer. Volume 1, second edition, 680 pages, diagrams, 6 X 9 in. New York, John 


Wiley & Sons, Inc., 1959. Price, $11.75. 


SAMPLING INSPECTION TABLES, by Harold F. Dodge and Harry G. Romig. Second edition, 
224 pages, 84 K 11 in. New York, John Wiley & Sons, Inc., 1959. Price, $8.00. 


A Way oF LIFE AND OTHER SELECTED WRITINGS, by Sir William Osler. Unabridged and 
unaltered republication of the work first published in 1951 by Oxford University Press 
under the title “Selected Writings of Sir William Osler."" 278 pages, illustrations, 54 X 8 


in. New York, Dover Publications, 1959. 


Price, $1.50 (paper). 


SCIENTIFIC RussIAN, by George E. Condoyannis. 225 pages, 5 X 7 in. New York, John 


Wiley & Sons, Inc., 1959. Price, $3.50. 


MANUAL OF SCIENTIFIC RusstaN, by Thomas F. Magner. 101 pages, 8} K 10} in. New 


York, Prentice-Hall, 1959. Price, $3.95. 


STEREO—How It Works, by Herman Burstein. 224 pages, illustrations, 54 X 84 in. New 
York, Gernsback Library, Inc., 1959. Price: $2.90 (paper) ; $5.00 (hard cover). 


Mars RevisitED, by Donald L. Cyr. 131 pages, illustrations, 43 X 7 in. Philadelphia, 


Dorrance & Co., Inc., 1959. Price, $3.06. 


FUNDAMENTALS OF PrpE DRAFTING, by Charles H. Thompson. 66 pages, illustrations, 
8} X 11 in. New York, John Wiley & Sons, Inc., 1958. Price, $3.50 (paper). 


PROCEEDINGS OF THE EIA CONFERENCE ON MAINTAINABILITY OF ELECTRONIC EQUIPMENT. 
104 pages, illustrations, 8} X 11 in. New York, Engineering Publishers, 1958. Price, 


$5.00 (paper). 


THE PIROTECHNIA OF VANNOCCIO BiRINGUCCIO, translated from the Italian with an intro- 
duction and notes by Cyril Stanley Smith and Martha Teach Gnudi. Reissue of the 
1943 printing, 477 pages, illustrations, 7} X 10} in. New York, Basic Books Inc., 1959. 


Price, $8.50. 
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FUNDAMENTALS OF ELECTRON DEVICES AND 
Circuits, by Herman R. Weed and Wells 
L. Davis. 591 pages, illustrations, 6 X 9 in. 
New York, Prentice-Hall, Inc., 1959. Price, 
$9.50. 


This new text was classroom tested for three 
years before publication. It was successfully 
used with both electrical and non-electrical 
students. Prerequisites are a working knowl- 
edge of calculus and a conventional back- 
ground in electrical physics. Each of the 
twenty-one chapters is a unit covering one 
type of device or circuit. The newer devices 
such as transistors, magnetic amplifiers and 
diodes are interspersed with the older ones 
such as gas tubes and regulators. Each chap- 
ter ends with a set of problems which em- 
phasize fundamental principles rather than 


specific applications. 


Locic or Screntiric Discovery, by 
Karl R. Popper. Translated from the Ger- 
man edition by the author. 480 pages, 
6 X 9} in. New York, Basic Books, Inc., 
1959. Price, $7.50. 


First published in German in 1934, this 
book now appears in English, translated by 
the author and augmented by 150 pages of 
new material (principally footnotes and ap- 
pendices), identified by asterisks. His theme 
originally was that there is logic in scientific 
discovery and that philosophy—in spite of 
many opinions to the contrary—can pose 
genuine problems about the whole of scien- 
tific knowledge. Twenty-five years later, Dr. 
Popper, in his additions to the original work, 
brings up-to-date some of his earlier concep- 
tions, without basically altering his original 
theme. The new material includes twelve ap- 
pendices covering such topics as probability, 
statistics and imaginary experiments. 

The final sentences in the original text sum 
up the author's conclusions: ‘Science never 
pursues the illusory aim of making its answers 
final, or even probable. Its advance is, rather, 
towards the infinite yet attainable aim of ever 
discovering new, deeper, and more general 
problems, and of subjecting its ever tentative 


answers to ever renewed and ever more rigor- 


ous tests.” 
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CALCULATEURS ANALOGIQUES by 
Rajko Tomovic. 186 pages, diagrams, 
6} X 94 in. Paris, Masson et Cie, 1958. 
Price, 3000 fr. 

This is the first of a series of monographs to 
be published under the auspices of |’ Associa- 
tion Internationale pour le Calcul analogique. 
The author, Director of Research at the Boris 
Kidrich Institute of Nuclear Sciences in Bel- 
grade, divides his discussion of repetitive ana- 
log computers into three main parts: theory, 
construction and applications. In Part I (four 
chapters), he covers the principles of operation 
of a repetitive differential analyzer and the 
theory of error analysis. Part II (six chap- 
ters) deals with the componet parts of the 
computer and with installation problems, in- 
cluding sources of errors in the various com- 
ponents, both for linear and non-linear opera- 
tion. Part III (three chapters) is devoted to 
applications for ordinary differential equations 
and for integral equations. A short chapter 
on high speed analyzers contains two interest- 
ing photographs of sets of solutions for a 
differential equation—one with 32 given 
points and the other with 32 values for one 
parameter. 


“ARCINOGENESIS BY ULTRAVIOLET LiGutT, by 
Harold F. Blum. 340 pages, diagrams, 
54 X 84 in. Princeton, Princeton Univer- 
sity Press, 1959. Price, $6.50. 


This book represents over fifteen years of 
experimental and theoretical studies on the 
effect of ultraviolet light in inducing cancer. 
It discusses various effects of ultraviolet light 
on biological systems as well as experimental 
carcinogenesis with this agent. A description 
of the course of carcinogenesis is based on 
experimental data in which cancer is treated 
Having arrived at 


as a growing tissue mass. 
this description, the author does not hesitate 
to speculate boldly—always being careful to 
remain within the limits set by the experi- 
mental evidence and by modern concepts of 
the dynamics of living systems. A new view- 
point with regard to cancer emerges. The re- 
lation of ultraviolet light to skin cancer in man 
is discussed together with some practical 
matters concerning environmental carcino- 
genic agents in general. 


| 
ae 


Sept., 1959.] Book 


PropuctTION ContROL, by Nyles V. Reinfeld. 
339 pages, diagrams, 6 X 9 in. New York, 
Prentice-Hall, Inc., 1959. Price, $8.00. 


Based on an extensive survey, this volume 
presents the results of seminars with more 
than 500 companies throughout the United 
States and Canada. The questions and an- 
swers emerging from round-table discussions 
form the nucleus of this practical and flexible 
guide. Viewing production control as a train- 
ing ground for top management, the author 
shows why this , remise holds true for all types 
and various sizes of organizations. Every 
phase of production control is carefully linked 
to management problems and methods. An 
up-to-date discussion of financial accounta- 
bility and management is provided by the 
author, giving the reader invaluable back- 
ground material. 


OPERATIONS RESEARCH, by Maurice Sasieni, 
Arthur Yaspan and Lawrence Friedman. 
316 pages, diagrams, 6 X 9 in. New York, 
John Wiley & Sons, Inc., 1959. Price, 
$10.25. 


Much of the subject matter of this text 
arose from a problems course in operations 
research, designed as a companion course to 
one offered at Case Institute of Technology. 
This work can serve as a textbook for a one- 
semester introductory techniques course, at 
either graduate or advanced undergraduate 
level. The reader is assumed to have a work- 
ing knowledge of the differential and integral 
calculus. Students majoring in such fields as 
industrial engineering, statistics, economics 
and applied mathematics will find the text 
invaluable. 


PROGRAMMING Business COMPUTERS, by 
Daniel D. McCracken, Harold Weiss and 
Tsai-Hwa Lee. 510 pages, diagrams, 6 K 9 
in. New York, John Wiley & Sons, Inc., 
1959. Price, $10.25. 

This volume is written primarily for those 
involved in the day-to-day application of elec- 
tronic computers to business data processing 
problems. The book begins with a discussion 
of fundamental topics, such as the nature of 
the data processing problem, the central con- 
cept of the file, flow charting and the general 
characteristics of electronic computers. The 
authors employ numerous examples to explain 
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all the standard techniques of coding. These 
examples are written in terms of a hypothetical 
computer called DATAC, which is a compila- 
tion of the features of many machines. The 
reader who understands the principles pre- 
sented in terms of DATAC, will find no diffi- 
culty in applying them to real machines. The 
concluding chapters deal with a summary of 
the steps involved in establishing a computer 
application, and a critical examination of the 
accounting and auditing problems associated 
with electronic data pro.essing. 


RADIOISOTOPES FOR INDUSTRY, by Robert S. 
Rochlin and Warner W. Schultz. 190 
pages, illustrations, 5 X 7} in. New York, 
Reinhold Publishing Corp., 1959. Price, 
$4.75. 


The purpose of this small handbook is to 
show the benefits that industry can derive by 
the use of radioisotopes. The reader will dis- 
cover techniques that can be used directly, or 
be stimulated in devising analogous solutions 
to his own problems. Chapters include intro- 
ductions to the subjects of instrumentation, 
facilities and safety. The book is not organ- 
ized in textbook form, but is intended to be 
useful to engineering students as well as engi- 
neers who wish to keep abreast of our advanc- 
ing technology. Assuming no prior knowledge 
of radioactivity and nuclear science, the author 
aims to introduce the reader to the benefits of 
radioisotopes with a minimum of unfamiliar 
jargon or highly technical discussion. Newer 
applications and techniques have been de- 
scribed at greater length, since details on the 
older ones can readily be found elsewhere. 


AsBeEsTos: Its INDUSTRIAL APPLICATIONS, by 

D. V. Rosato. 213 pages, illustrations, 

5 X 7} in. New York, Reinhold Publishing 

Corp., 1959. Price, $5.75. 

This text represents an accumulation of all 
data available on the properties, manufactur- 
ing methods and uses of asbestos in industry. 
Included also are the types of asbestos ma- 
terial available—asbestos cement, tile, asbes- 


tos heat and electrical insulation, friction 
materials, textiles, plastics, packings, gaskets, 
filters and others. A census of asbestos prod- 
ucts is also given. The volume serves as a con- 
venient guide to various industrial branches, 
research and development groups, manufac- 


turers, and market research groups. 
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Giant Balloon to Study Cosmic 
Rays.—The first in a new series of in- 
strumented balloons was launched on 
June 29 by the Air Research and De- 
velopment Command of the United 
S.ates Air Force near Mineral Wells, 
Texas, in an effort to gather new data 
about cosmic rays near the top of the 
Earth’s atmosphere. 

The 128-ft. unmanned balloon, 
launched by personnel from Holloman 
Air Force Base under the direction of 
Major L. M. Bogard, Chief, Balloon 
Branch, Air Force Missile Develop- 
ment Center, carried a gondola with 
about 100 pounds of scientific appa- 
ratus and was designed to reach an 
altitude of about 100,000 ft. The in- 
struments, developed jointly by RIAS 
(a research division of The Martin 
Company located near Baltimore, 
Md.) and the Air Force Cambridge 
Research Center (Bedford, Mass.), are 
designed to measure only a small por- 
tion of the cosmic particles which strike 
the gondola—the so-called ‘heavy 
nuclei.” 

Since “‘primary’’ cosmic rays disin- 
tegrate when they collide with atoms 
of gas in the Earth’s atmosphere, pro- 
ducing and transferring their energy to 
“secondary’’ cosmic particles, a bal- 
loon at this altitude is expected to en- 
counter both primary and secondary 
radiation. Both will be detected, but 
a unique device developed by RIAS 
will distinguish between the two by 
noting the change in velocity as parti- 
cles pass through a beryllium plate 
and record them separately. Findings 
are radioed back to Earth. 

Most of the particles in cosmic rays 
are the nuclei of hydrogen atoms origi- 
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nating from our sun and other cosmic 
bodies. Some of them, however, are 
the nuclei of heavier atoms, and these 
are the ones on which the current ex- 
periments will concentrate. These 
particular particles are believed to be 
especially sensitive to solar influences, 
and scientists will watch closely for 
variations in their intensity during 
each 10- to 20-hr. flight. 

The heart of the cosmic ray experi- 
ment is a metal cylinder about the size 
of a portable fire extinguisher, filled 
with argon gas under pressure. A wire 
carrying a positive electrical charge 
runs the length of a cylinder and the 
plate of beryllium bisects it. When a 
high-speed cosmic particle enters the 
chamber it “‘ionizes’’ some of the gas 
atoms inside (that is, it knocks an elec- 
tron off the outer shell of each atom). 
These free electrons are attracted to 
the wire at the center of the cham- 
ber and counted. By measuring the 
amount of ionization, scientists can 
determine the precise nature of the 
particle that caused it. 

The balloons will be tracked con- 
tinuously, and eventual recovery will 
be made. 


New Development in Cancer Treat- 
ment.—A radiotherapy unit for the 
treatment of cancer in an atmosphere 
of high pressure oxygen has been 


opened at St. Thomas’ Hospital, 
London. 

The unit, which includes a pressure 
chamber in the form of an 8-ft. cylin- 
der in which the patient lies, is the 
only one of its kind in the world. It 
was installed at a cost of nearly 


$140,000 by a British nuclear engineer- 
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ing firm, after development at St. 
Thomas’ by a team consisting of radio- 
therapists, an anaesthetist, a patholo- 
gist, a physicist and an electronics 
engineer. 

The team has pioneered the treat- 
ment to its present degree and is being 
followed by a team at University Col- 
lege Hospital, London. Teams in the 
United States, Europe and Australia 
are also following this lead. 

Treating cancer by radiotherapy is 
based on the fact that all living cells 
can be killed by the action of ionizing 
radiations and there is a relationship 
between the sensitivity of cells to 
damage by them and the concentra- 
tion of oxygen present. 

During the time of irradiation the 
patient cannot be seen and nobody can 
remain in the treatment room. The 
patient is given an anaesthetic and en- 
closed in the pressure chamber. 

Observations of the patient are made 
by continuous electrical readings of 
the heart beat, respirations, blood flow 
and pressure, and of brain-waves. 

The patient lies on a couch which 
runs into the chamber on rails. Radi- 
ation is applied through the walls of 
the chamber. Twelve large cylinders 
supply oxygen. The walls of the 
treatment room are of dense concrete 
up to 2-ft. thick, the door is lead-lined 
and the observation window is a lead 
glass block 5 in. thick. 

The technique has been used on 
more than 80 patients without mishap. 
The members of the British team will 
say no more than that they are en- 
couraged by the results, but it is hoped 
the apparatus will become a perma- 
nent contribution towards the treat- 
ment of cancer. 


Flat Speaker Developed.— An elec- 
tro-dynamic loudspeaker of revolu- 
tionary design has been developed in 
the Electronics Department of the 
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Weizmann Institute of Science in Re- 
hovot, Israel. The main outward fea- 
ture of the invention—called an Iso- 
phase speaker—is its flatness. The 
paper cone, which is part of any con- 
ventional speaker, has disappeared. 

The new device, which is expected 
to be of great advantage in high fidel- 
ity, stereophonic and television sys- 
tems, utilizes a membrane with a 
printed zig-zag winding to convert the 
current into vibratory motion, in place 
of the conventional mobile coil. 

One of the main advantages of the 
new amplifier is that the phase of the 
acoustic waves is completely true and 
is governed by the phase of the electric 
current. Also, owing to its flatness, 
the new speaker needs much less space 
than the traditional type. 


Multi-Tap Lumped Delay Line.- 
The development of a new lumped- 
constant Delay Line having 10 sepa- 
rate taps of 0.05 usec. each providing a 
maximum total delay of 0.5 usec. has 
been announced by Technitrol Engi- 
neering Company, Philadelphia. 

Designated as CT-18, the Delay 
Line has an impedance of 550 ohms 
and a maximum rise time figure of 0.1 
usec. Attenuation is 1.0 db. max. and 
temperature range is from —25 C. to 
+85 C. Some useful applications for 
a delay line of this type would be as a 
trimmer for a long delay line (5 usec. 
or more) to adjust the tolerance of the 
line ; as a variable pulse generator with 
the taps controlling the width of the 
output pulse ; to generate clock pulses ; 
or asa variable pulse forming network. 

Delay Line CT-18 is 2} in. long, lin. 
wide and approximately 13 in. high in- 
cluding threaded mounting studs on 
the underside of the case. 


Automatic Accounting Machines. 
A complete new line of automatic ac- 
counting machines that will obsolete 
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thousands of conventional bookkeep- 
ing machines in use today has been 
developed by Burroughs Corporation. 
Ray R. Eppert, Burroughs president, 
said the new series of electro-mechani- 
cal machines for business and industry 
represents one of the most important 
advances in conventional bookkeeping 
equipment in the past decade. Ten 
models of the machine were placed on 
the market August 1. 

Designed for all types of businesses 
—from dress shops to public utilities— 
the general purpose machine is the first 
fully automatic mechanical accounting 
device equipped with a ‘‘dual printer,” 
a special feature that increases oper- 
ating speeds by as much as 50 per cent. 

The new device, a numeric data 
processing accounting machine, Series 
F5000, is an outgrowth of Burroughs’ 
famed Sensimatic. Like the Sensi- 


matic, the new line fills the middle 
ground between the company’s small 
portable bookkeeping machines and 


high speed electronic data processing 
equipment, Eppert declared. 

Offering a new standard of economy 
and versatility, the machine is de- 
signed to improve all bookkeeping op- 
erations that require two original rec- 
ords, such as a statement and ledger. 
The machine prints both records si- 
multaneously. 

Accounts receivable, accounts pay- 
able, payroll, inventory, bank book- 
keeping and posting to subsidiary 
ledgers are among the scores of im- 
portant accounting tasks the machine 
will handle with greater speed and 
accuracy. 

One of the machine’s two printing 
heads can be de-activated for account- 
ing functions that only require a single 
printing, the president explained. 

All models will be equipped with a 
unique “program panel”’ that controls 
the automatic functions of the ma- 
chine, telling it to add, subtract, retain 
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figures for later use, print out totals 
and the like. The panel is capable of 
directing the machine to perform up 
to 35,000 functions without assistance 
from the operator. 

The compact device can be “pro- 
grammed” to perform four totally 
different accounting jobs with a simple 
turn of a knob. 

Other features include modern styl- 
ing, conventional keyboard input, con- 
trolled ‘‘gang printing’’ output, two to 
19 memory-storage units with a maxi- 
mum capacity of 228 digits, and auto- 
matic readout of data accumulated in 
storage units. 

The machine not only will reduce 
the time needed to process mounting 
volumes of paper work, but also will 
permit other economies through sim- 
plification of accounting forms. 

Credit unions, department stores, 
banks, insurance brokers, manufac- 
turers, wholesalers, retailers, federal, 
state and local government agencies— 
virtually every business that keeps 
books—will be able to realize substan- 
tial savings, Eppert concluded. 


Microwave Techniques to Improve 
Computers.—New electronic comput- 
ing techniques, approaching the oper- 
ating limits imposed by the speed of 
light, are under development by the 
Radio Coporation of America for use 
in computers that will function up 
to 1,000 times faster than present 
systems. 

Announcement of the advances was 
made recently by E. D. Foster, Vice 
President and General Manager, RCA 
Electronic Data Processing Division, 
and Dr. James Hillier, Vice Presi- 
dent, RCA Laboratories. The two 
RCA divisions are collaborating in 
an ultra-high-speed computer develop- 
ment study program under a contract 
with the U.S. Navy’s Bureau of Ships. 

Dr. Hillier predicted that such fu- 
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ture computers will be not only more 
compact and economical than today’s 
systems, but will be capable of han- 
dling problems toocomplex for solution 
within a useful time by present tech- 
niques. Such problems, he added, 
might relate to ‘‘the frontiers of practi- 
cally every science, as well as many so- 
cial, economic and strategic matters.” 

The two executives stated that the 
new techniques now under develop- 
ment by RCA will result in ‘‘a revolu- 
tionary broadening of electronic com- 
puter functions, with profound impact 
upon science, industry, and defense.”’ 

“The startling advances achieved 
through this and similar projects will 
find their place during the years im- 
mediately ahead in the commercial 
data processing equipment which RCA 
will make available to commerce and 
industry as well as to the military 
services,’’ said Mr. Foster. 

Experimental use of some of the new 
techniques in the laboratory already 
has proven their ability to perform 
computations such as adding two digits 
together at speeds of 500 million times 
a second, according to Dr. Hillier. 
Stating that the techniques are based 
upon microwave principles employed 
in super-high-frequency radar and 
radio systems and upon new develop- 
ments in solid-state electronics, he 
added that continuing development is 
expected to lead soon to ‘‘a considera- 
ble further increase”’ in speed. 

“Since the end of World War II, we 
have progressed from electromechani- 
cal systems using devices capable of 
less than 100 basic switching opera- 
tions per second to electronic com- 
puters with operating speeds more 
than a thousand times faster,”’ said 
Dr. Hillier. 

“This advance already has had 
major effects upon our science and 
technology by permitting solution of 
problems far too complex for the earlier 
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systems. These problems relate to 
aerodynamic design, weather predic- 
tion, missile guidance, and many as- 
pects of research in physics, chem- 
istry, electronics and other sciences. 
At the same time, one of the most 
significant results of this increase in 
computing speeds has been the de- 
velopment and broadening application 
of electronic data-processing systems 
in commerce and industry. 

“The new ultra-high-speed tech- 
niques now under development by 
RCA are aimed at increasing the com- 
putation rate by as much as 1,000 
times beyond today’s electronic com- 
puter speeds. With such an achieve- 
ment, we can expect a new revolution- 
ary effect at least equal to that result- 
ing from the advances of the past 
decade.” 


Scope of Development Program 


The research and development pro- 
gram on ultra-high-speed computer 
techniques involves scientific and en- 
gineering groups of several RCA divi- 
sions. These include not only the 
RCA Electronic Data Processing Di- 
vision and RCA Laboratories, but also 
the RCA Electron Tube and RCA 
Semiconductor and Materials Divi- 
sions. The RCA Laboratories’ con- 
tribution discussed here has consisted 
of a large-scale research effort involv- 
ing specialists in microwave, semicon- 
ductor and computer techniques, di- 
rected by Dr. Jan A. Rajchman, an 
Associate Director of the Technical 
Staff at the David Sarnoff Research 
Center, Princeton, N. J. 

According to Dr. Rajchman, em- 
phasis in the research program has 
been placed upon new methods of per- 
forming all basic computer functions, 
such as switching, storing and ampli- 
fying the coded electrical signals that 
represent the “language”’ of electronic 
computers. 
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He stated that one promising 
method, out of several that are under 
study and development, involves a 
combination of microwave and solid- 
state techniques in a possible system 
based on the use of super-high-fre- 
quency principles. Dr. Rajchman ex- 
plained that a computing system based 
on microwave techniques would handle 
information in the form of frequency 
pulses varying from one to ten billion 
cycles per second, in the following 
manner : 

The language of the computer would 
consist of signals in the form of elec- 
trical oscillations at these frequencies. 
The oscillation in one phase represents 
the digit ‘1,’’ while an oscillation of 
the same frequency in the opposite 
phase—that is, having a directly op- 
posite wave pattern of peaks and val- 
leys—would represent In the 
binary code system used in computers, 
all numbers can be expressed by vari- 
ous combinations of these two digits, 
so that all numerical information could 
be conveyed through the system en- 
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tirely by means of the two opposite 
phases of the high frequency signal. 
With the super-high-frequencies, ex- 
tremely high computing speeds can 
thus be achieved. 

“At the operating speed achieved 
experimentally in the laboratory, light 
travels only two feet in the time re- 
quired for a simple computation such 
as the adding of two digits,” said 
Dr. Rajchman. “Electrical impulses 
travel through wire at slightly less 
than the speed of light. Thus the 
speed of light, which normally might 
appear incredibly fast, is in this case 
so slow that it imposes design limita- 
tions upon our systems. The circuit 
connections must be made as short as 
possible so that the gain in operating 
speed is not wasted by a time lag in 
communicating information from one 
point to another in the system. It 
follows from this that circuits designed 
for computers based on microwave 
techniques must be miniaturized as 


far as possible, leading in turn to 
smaller as well as faster systems.” 


ic 
pind 
sf 
3 


: 
: 


THE FRANKLIN INSTITUTE LABORATORIES 


FOR RESEARCH AND DEVELOPMENT 


In the fields of engineering and the physi- 
cal sciences, a competent and versatile 
staff of several hundred scientists and en- 
gineers, working with modern equipment 
in a creative climate, is trained to bring a 
fresh scientific approach to the solution of 


industrial problems. 


THE FRANKLIN INSTITUTE OF THE STATE OF PENNSYLVANIA 


BENJAMIN FRANKLIN PARKWAY AT 20TH ST. PHILADELPHIA 3, PA. 


Rte 


: 


